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abstract 


In  order  to  determine  the  effects  of  the  continental  shelf  bathy¬ 
metry  on  coastal  geomorphology,  a  series  of  wave-refraction  diagrams 
were  generated  for  the  S.  C.  coast  from  Bull  Island  to  the  Isle  of 
Palms.  Wave  rays,  as  they  approach  the  shore,  converge  or  diverge 
depending  on  the  uneven  offshore  bottom  topography.  Therefore,  zones 
of  magnified  or  reduced  wave  energy  are  created  from  the  interaction 
between  waves  and  the  offshore  topography.  REFRAC,  a  computerized 
wave- refract ion  program  developed  for  this  study,  generates  refrac¬ 
tion  diagrams  which  delineate  the  patterns  of  longshore  variation  in 
wave  energy.  A  variety  of  input  wave  conditions  were  used  in  REFRAC 
to  model  the  various  possible  wave  conditions  existing  in  nature. 

Data  input  included  waves  propagated  in  deep  water  from  the  east, 
southeast  and  south  for  several  different  periods.  To  improve  the 
accuracy  of  the  refraction  diagrams,  bathymetric  charts  of  increasing 
detail  near  the  shore  were  used  in  tracing  the  path  of  a  wave  onto 
the  shore.  Qualitatively,  zones  of  potential  erosion  or  deposition 
can  be  inferred  to  correspond  to  converging  or  diverging  wave  rays 
respectively. 

The  results  indicate  that  the  offshore  bathymetry  does  partially 
control  the  coastal  geomorphology  by  creating  zones  of  potential  ero¬ 
sion  and  deposition  and  by  influencing  the  direction  of  sediment 
transport.  An  analysis  of  the  refraction  diagrams  reveals  the  fol¬ 
lowing  observations.  (1)  Bull  Island  and  Capers  Island,  areas  of 
long-term  erosion,  are  located  in  zones  of  higher  than  average  wave 
energy.  (2)  The  southern  section  of  the  Isle  of  Palms  has  under¬ 
gone  extensive  accretion  due  to  its  location  in  a  zone  of  lower 


than  average  wave  energy  and  to  a  sediment  supply  from  the  north.  (3) 
The-oblique  orientation  of  Dewees  Inlet,  as  compared  to  the  normal  orien¬ 
tation  of  Price  Inlet,  results  from  the  increased  wave  energy  noted  at 
Dewees  Inlet,  especially  for  a  10-second  wave  from  the  east.  (4)  The 
large  downdrift  offset  of  Dewees  Inlet  appears  to  be  related  to  a  sud¬ 
den  reduction  in  the  southward  directed  longshore  drift  at  the  Isle  of 
Palms  as  compared  to  Dewees  and  Capers  Islands. 

Additional  information  on  the  sediment  transport  patterns  at  the 
ebb-tidal  delta  of  Price  Inlet  was  gained  from  the  increase  in  detail 
of  the  refraction  diagrams  for  Price  Inlet.  For  waves  from  the  east 
and  southeast,  the  predominant  direction  of  littoral  drift  is  to  the 
south.  Sediment  transport  reversals  to  the  north  were  seen  on  the 
south  side  of  the  inlet.  These  reversals  in  transport  direction,  re¬ 
sulting  from  refraction  around  the  ebb-tidal  delta,  are  capable  of 
reintroducing  sand  into  the  inlet  and  building  up  the  beach  and 
shoals  south  of  the  inlet.  On  the  other  hand,  shorter  period  waves 
from  the  east,  because  of  their  oblique  angle  of  approach  to  the 
shoreline,  may  enhance  sediment  bypassing  around  the  distal  portion 
of  the  ebb-tidal  delta  at  Price  Inlet. 
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INTRODUCTION 


The  relationship  between  the  nonuniform  wave-energy  distribution 
along  a  coastline  and  the  coastline's  stability,  both  past  and  future, 
is  of  primary  importance  in  coastal  management.  A  graphical  method, 
based  on  the  principles  of  refraction,  set  forth  by  Snell's  Law,  is 
commonly  used  in  displaying  wave-energy  distribution.  These  wave-re¬ 
fraction  diagrams  serve  as  a  useful  tool  in  t^e  interpretation  of 
coastal  processes  when  used  in  conjunction  with  field  data  (Colonell 
et  al. ,  1973).  Goldsmith  and  others  (1970)  concluded  from  refraction 
computations  on  Monomy  Island  that  the  interaction  between  incoming 
waves  and  variations  in  the  offshore  bathymetry  may  be  able  to  pre¬ 
dict  changes  in  the  configuration  of  the  shoreline.  Fico  (1977) 
depicted  areas  on  the  S.  C.  coast  of  varying  susceptibility  to  storm 
damage  by  refracting  waves  from  several  different  directions. 

This  study  exemplifies  the  significance  of  the  offshore  bathy¬ 
metry  in  modifying  the  configuration  of  the  S.  C.  coast  from  Bull  Is¬ 
land  to  the  Isle  of  Palms  (Fig.  1).  Through  the  interaction  of  waves 
and  the  offshore  topography,  which  results  in  an  uneven  longshore 
wave  energy  distribution,  differential  rates  of  deposition  or  erosion 
occur  along  the  coast.  A  series  of  wave-refraction  diagrams  were 
generated  to  show  the  correspondence  between  the  variation  in  the 
wave-energy  distribution  and  the  resulting  coastal  geomorphology. 
Furthermore,  detailed  refraction  diagrams  of  Price  Inlet  were  able 
to  illustrate  the  effects  of  the  offshore  bathymetry  on  the  morpho¬ 
logy  and  sediment  transport  patterns  of  an  ebb-tidal  delta. 

A  computerized  wave  refraction  procedure,  REFRAC,  has  been  im¬ 
plemented  on  the  IBM  System  370/168  at  the  University  of  South 


! 
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Carolina.  Dobson  (1967)  developed  the  initial  program  and  based  the  re¬ 
fraction  of  incoming  waves  on  linear  wave  theory.  Since  he  assumed  bot¬ 
tom  friction  to  be  insignificant,  the  wave  height  at  any  point  of  inter¬ 
est  (usually  the  breaker  zone)  is  a  function  of  the  initial  deep  water 
wave  height  and  the  refraction  and  shoaling  coefficients  at  that  point. 
Modifications  of  this  program  made  by  Senter  (1972)  at  the  Waterways  Ex¬ 
periment  Station  include  the  addition  of  Calcomp  plot  routines  and  a 
window  feature  which  allows  wave  data  generated  from  a  small-scale  map 
to  be  used  as  input  into  a  larger-scale  map  for  a  more  detailed  refrac¬ 
tion  analysis. 

REFRAC,  including  the  above  modifications,  has  been  further  al¬ 
tered  in  the  current  study  to  accept  digitized  bathymetric  input  for 
generation  of  a  depth  grid.  Previously,  areas  of  interest  were  over- 
lain  by  a  grid  and  the  depths  interpolated  from  a  hydrographic  chart 
at  grid  points.  Next,  the  interpolated  depth  data  were  keypunched 
onto  cards.  When  using  large  areas  or  a  series  of  desired  large- 
scale  window  plots,  this  was  a  very  time-consuming  process.  These 
tasks  are  eliminated  by  digitizing  the  shoreline  and  the  bathymetric 
contours  onto  magnetic  tape.  The  digitized  data  are  processed  by  the 
program  which  creates  a  grid  of  depth  values  by  using  the  straight-line 
slope  formula  for  depths  between  the  contours. 


LINEAR  WAVE  THEORY 


General  Description 

Linear  wave  theory,  as  developed  by  Airy  (in  1845)  is  a  first-order 
theoretical  approximation  describing  wave  behavior.  The  theory  assumes 
long  trains  of  uniform  waves  with  long  unbroken  wave  crests.  Obviously, 
this  condition  is  not  met  in  areas  of  locally  generated  seas,  where 
wave  forms  may  consist  of  several  periods  and  directions  and  broken  wave 
crests.  However,  swell  waves,  which  are  outside  the  area  of  generation, 
approximate  the  above  assumption  (Komar,  1976).  Higher-order  theories, 
usually  referred  to  as  finite  amplitude  theories,  increase  the  number  of 
successive  approximations  in  order  to  describe  more  closely  the  actual 
wave  behavior.  These  successive  approximations  serve  as  correction  fac¬ 
tors  for  preceding  terms  (Coastal  Engineering  Research  Center,  1973, 

p.  2-2). 

Linear  wave  theory  is  more  commonly  used  in  the  study  of  wave  be¬ 
havior  because  it  is  reliable  over  a  number  of  wave  conditions  and 
mathematically  easy  to  apply.  Dobson  (1967),  using  linear  wave  theory 
in  the  refraction  of  incoming  waves,  developed  the  original  version 
of  the  computer  program  used  in  this  study.  Therefore,  the  celerity 
of  the  wave  is  dependent  on  the  water  depth.  Refraction  or  the  bending 
of  wave  orthogonals  (rays  or  lines  parallel  to  the  direction  of  wave 
propagation)  are  a  result  of  the  different  speeds  along  the  wave  crest 
as  it  passes  over  an  uneven  topography.  Simply  stated,  a  wave  ray 
will  refract  toward  shallower  water.  Because  of  this  interaction  be¬ 
tween  waves  and  the  offshore  topography,  there  is  a  variation  in  wave 
characteristics  (height,  energy,  etc.)  in  the  nearshore  zone.  This 
spatial  variation  in  wave  characteristics  may  be  a  driving  mechanism 
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for  nearshore  circulation  patterns  (Noda,  1974;  Sonu,  1972).  Qualita¬ 
tively,  zones  of  potential  erosion  or  deposition  can  be  inferred  to 
correspond  to  converging  or  diverging  wave  rays,  respectively. 

Breaker  wave  height  is  a  function  of  the  deep  water  wave  height  af¬ 
ter  it  has  been  modified  by  refraction  and  shoaling  in  shallow  water. 
Although  Dobson  considered  frictional  attenuation  of  wave  height  to  be 
insignificant.  Goldsmith  (1976)  states  that  friction  may  cause  a  sig¬ 
nificant  loss  in  wave  energy  and  height.  Because  of  the  frictional 
loss  of  energy  over  the  wide  shallow  shelf  of  the  Virginian  Sea,  wave 
heights  were  reduced  50  to  75  percent  for  longer  wave  periods.  This 
percentage  will  be  less  for  smaller  period  waves. 

A  friction  routine  similar  to  the  one  used  by  Goldsmith  (personal 
communication)  and  Coleman  and  Wright  (1971)  was  added  to  REFRAC,  the 
wave-refraction  program  used  in  this  study.  However,  the  resulting 
wave  heights  and  energy  are  questionable.  Therefore,  a  quantitative 
evaluation  of  wave  heights  was  not  attempted. 

Assumptions  and  Limitations 

In  relation  to  wave  refraction,  the  main  assumptions  in  linear 
wave  theory  are  (Coastal  Engineering  Research  Center,  1973,  p.  2-65): 

(1)  Wave  energy  between  wave  rays  remains  constant.  This  assump¬ 
tion  is  suspect  when  wave  orthogonals  bend  sharply  as  energy  may  be 
transferred  along  the  wave  crests.  Caustics,  which  are  rays  that 
bend  sharply  enough  to  cross,  no  longer  present  a  problem.  Accord¬ 
ing  to  Chao  (1972),  wave  rays  continue  on  the  same  path  after  they 
pass  through  a  caustic  as  before  the  caustic,  the  only  difference  be¬ 


ing  a  phase  shift. 


(2)  Wave  celerity  is  a  function  of  water  depth.  This  holds  true  in 
linear  theory  but  not  necessarily  in  higher-order  theories.  An  increase 
in  wave  height,  resulting  from  either  shoaling  or  refraction  will  cause 
a  slight  increase  in  the  wave  velocity.  This  effect  is  small  (Komar, 
1976). 

(3)  The  bottom  slope  is  gentle  (less  than  1:10).  Linear  theory  is 
strictly  valid  only  for  constant  depths,  but  it  will  successfully  pre¬ 
dict  wave  velocities  over  a  gently  sloping  bottom  (Dobson,  1967). 

(4)  Waves  are  long-crested,  constant  period,  and  of  small  amplitude. 
This  is  not  true  for  'confused  seas’  in  their  area  of  generation  but 
does  approximate  swell  conditions.  In  refraction  studies,  a  spec¬ 
trum  of  swell  conditions  needs  to  be  analyzed  to  simulate  the  real  world 
(Goldsmith,  1976). 

(5)  Reflection  of  wave  energy  from  a  gently  sloping  bottom  is  neg¬ 
ligible.  Caldwell  (1949)  supports  this  assumption  for  slopes  of  four 
degrees  or  less. 


DATA  INPUT 


South  Carolina  Wave  Climate 

Brown  (1977)  used  SSMO  data  (U.S.  Naval  Weather  Service  Command, 

1970)  to  derive  wind  and  wave-energy- flux  diagrams  for  the  South  Caro¬ 
lina  coast.  As  expected,  seasonal  trends  in  wind  and  wave  conditions 

were  noted  (Fig.  2).  The  average  annual  wave  energy  flux  on  the  South 

6  3  6  3 

Carolina  coast  is  1.7  x  10  g-m/s  ,  with  a  maximum  of  3.1  x  10  g-m/s 

from  the  northeast.  Wave  energy  flux  is  defined  as  the  rate  at  which 

wave  energy  per  unit  surface  area  is  transmitted  in  the  direction  of 

wave  propagation. 

The  dominant  winds  (highest  velocities)  from  the  north  and  east 
during  fall  and  winter  are  reflected  in  the  direction  of  the  deep  water 
wave-energy  flux.  Strong  storm  winds  from  the  northeast  are  generated 
by  northward  passing  extratropical  storms  and  are  considered  by  Finley 
(1976)  to  be  the  "most  important  wave  generators".  Finley  documented 
7.3  meters  of  foredune  erosion  on  Debidue  Island  in  a  two-week  period 
following  an  extratropical  cyclone  in  February  1973.  Kana  (1977)  re¬ 
corded  process  measurements  during  a  minor  northeast  storm  in  September 

3 

1974.  An  average  of  3.8  m  of  sand  per  foot  of  shoreline  at  Debidue 
Island  was  eroded  in  a  6-hour  period.  During  this  short  period  of 
erosion,  breaker  heights  were  approximately  120  cm  (4  ft)  with  an  average 
wave  period  of  6  seconds. 

In  spring  and  summer,  an  increasing  frequency  of  winds  are  observed 
from  the  south  and  southwest.  These  winds  are  generated  by  an  anti- 
cyclonic  circulation  pattern  associated  with  a  high  pressure  zone  set¬ 
tling  over  Bermuda.  According  to  Crutcher  and  Quayle  in  Brown  (1977), 
an  average  of  1.4  hurricanes  and  tropical  storms  affect  South  Carolina's 
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Figure  2.  Seasonal  and  annual  wind  and  wave  energy  flux  (P^s) 
roses  computed  from  the  1970  version  of  SSMO  data.  Bar 
length  in  the  wind  roses  represents  the  percentage  of  time 
the  wind  blows  from  any  given  direction.  Bar  length  in  the 
energy  flux  roses  is  a  relative  measure  of  the  wave  energy 
coming  from  a  given  direction  (from  Brown,  1977). 
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coast  annually.  The  probability  of  a  tropical  storm  of  hurricane  force, 
winds  in  excess  of  74  mph,  striking  the  South  Carolina  coast,  was  de¬ 
termined  by  Nummedal  and  Humphries  (1977)  to  be  about  .2,  correspond¬ 
ing  to  one  every  5  years.  The  highest  breakers  recorded  on  the  coast 
were  approximately  12  ft.  They  were  recorded  at  Myrtle  Beach  during  a 
hurricane  in  1958  (Nummedal,  1977).  Nummedal  and  Humphries  (1977)  also 
noted,  from  numerous  observations  at  Debidue  Island,  an  11  ciu  decrease 
in  wave  height  for  the  spring  and  summer  as  compared  to  fall  and  winter. 

The  dominant  deep  water  wave-energy  flux  in  the  summer,  as  calcu¬ 
lated  by  Brown  (1977),  is  from  the  southeast  and  lower  in  magnitude 
than  the  winter  flux  values  from  the  northeast  and  east.  Finley  (1976) 
computed  the  highest  sea  and  swell  which  could  affect  the  S.  C.  coast 
to  be  from  the  northeast  and  east.  Figure  3  graphically  displays  the 
frequency  of  occurrence  of  sea  and  swell  and  their  approach  direction. 
These  SSMO-based  wave  climate  evaluations  support  the  morphologic  evi¬ 
dence  of  a  net  southward  sediment  transport. 

Based  on  the  above  information  and  1975  SSMO  data,  it  was  decided  to 
use  Goldsmith's  (1976)  "scatter  gun"  approach  with  respect  to  wave  in¬ 
put  conditions.  In  this  way,  a  variety  of  conditions  are  modeled.  Un¬ 
fortunately,  SSMO  data  have  several  inherent  biases  which  make  it  dif¬ 
ficult  to  determine  percentage  of  occurrence  for  a  given  wave  condition 
(Goldsmith,  1976;  Nummadal  and  Stephen,  1976). 

The  initial  deep  water  wave  conditions  consist  of  waves  propagating 
from  the  east  (90  degrees),  southeast  (135  degrees),  and  from  the  south 
(180  degrees)  at  periods  of  10,  8,  6  and  4  seconds  with  a  wave  height 
of  1  foot.  An  initial  wave  height  of  one  foot  was  used  at  the  time  of 
this  analysis  because  frictional  attenuation  of  the  wave  height  in 
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Figure  3.  Frequency  of  occurrence  of  sea  and  swell  and  their 
approach  direction  for  SSMO  observation  square  off  South  Caro 
lina  (from  U.S.  Army  Corps  of  Engineers,  1970).  (from  Finley 
1976). 


shallow  water  was  not  considered.  Therefore,  these  diagrams  are  of  qual¬ 
itative  value.  Calculations  of  energy  and  longshore  transport  can  be 
used  for  relative  comparisons  between  sections  of  the  coast  but  should 
not  be  used  as  indicators  of  absolute  magnitudes. 

Bathymetry 

Shelf  morphology.  -  According  to  Swift  (1976),  the  Atlantic  shelf 
sands  are  predominantly  generated  by  erosional  shoreface  retreat  during 
the  Holocene  transgression  about  11,000  B.P.  Two  constructional  fea¬ 
tures  formed  from  this  sand  sheet  were  shoal  retreat  massifs  which  are 
overlain  by  linear  sand  shoals  (Swift  et  al. ,  1972;  Duane  et  al. ,  1972). 
Refer  to  Figure  4. 

Shoal  retreat  massifs  are  broad  sand  ridges  of  subdued  relief,  trans¬ 
verse  to  the  shelf,  which  mark  the  retreat  of  nearshore  depositional  cen¬ 
ters.  These  depositional  centers  form  off  capes  or  cuspate  forelands 
and  are  the  result  of  littoral  drift  convergence  (Swift,  et  al. ,  1972). 
Because  of  the  closely-spaced  forelands  south  of  Cape  Romain,  the  shoal 
retreat  massifs  tend  to  coalesce. 

Overlying  the  shoal  retreat  massifs  are  northeast  trending  linear 
shoals.  These  shoals  form  an  angle  of  approximately  35  degrees  with 
the  shoreline,  exhibit  up  to  30  feet  of  relief  and  may  extend  for  many 
miles  (Duane  et  al. ,  1972).  The  shoals  may  be  connected  to  the  shore¬ 
line  or  isolated.  Duane  jst  al.  hypothesize  that  the  shoreface  con¬ 
nected  ridges  are  formed  by  storm-generated  currents  interacting  with 
the  shoreface.  These  ridges  become  isolated  as  the  shoreline  retreats 
in  response  to  sea  level  rise.  Because  of  the  similarity  in  orienta¬ 
tion  of  both  shoreface-connected  and  isolated  shoals  with  respect  to 
the  shoreline,  Duane  et  al^  propose  that  the  shoreline  orientation  pro- 


17 


Figure  4.  Cuspate  forelands  and  cape  shoal-retreat  massifs 
(stippled)  of  the  South  Carolina  shelf.  Note  overprinting 
of  ridge  and  swale  topography.  Contours  are  in  fathoms, 
(from  Swift,  1976). 


bably  remained  essentially  the  same  during  the  Holocene  marine  trans¬ 
gression. 

Hydrographic  charts.  -  Depth  data  is  obtained  by  contouring  hydro- 
graphic  charts  and  boat  sheets  printed  by  NOAA  (National  Oceanic  and  At¬ 
mospheric  Administration) .  The  contoured  depths  are  hand  digitized  on  a 
Bendix  Datagrid  Digitizer  and  stored  on  magnetic  tape.  This  tape  is 
read  by  a  routine  in  REFRAC,  which  constructs  a  depth  grid  b>  inter¬ 
polating  depths  between  the  digitized  contours.  (For  more  details,  re¬ 
fer  to  Appendix  II) 

A  first-order  chart  of  the  South  Carolina  shelf  (scaled  at  1:432, 

720  to  1:449,659)  was  used  to  genetate  input  data  onto  a  larger  scaled 
second-order  chart  (1:80,000)  of  the  section  of  coast  between  Cape  Ro- 
main  and  Folly  Island.  This  chart  generated  data  for  the  larger  scaled 
third-order  charts  (1:20,000)  of  Price  Inlet.  Refer  to  Fig.  5  and 
Table  1. 

The  advantage  of  this  technique  lies  in  the  use  of  more  detailed 
bathymetric  data  for  second  and  third  order  charts.  This  allows  in¬ 
creased  detail  and  greater  accuracy  in  the  resulting  wave-refraction 
diagrams . 

In  relation  to  depth  data,  two  considerations  are  necessary. 

First  is  the  accuracy  with  which  the  depths  were  measured.  Accuracy 
criteria  for  the  depths  and  navigational  positioning  has  been  compiled 
by  Sallenger  et  al.  (1975).  Refer  to  Figure  6  and  Table  2  for  a  gra¬ 
phical  summary  and  explanation  of  this  information.  Second,  the 
amount  of  area  distortion  produced  by  the  Mercator  projection  from  a 
sphere  to  a  flat  map  need  be  insignificant  or  corrected  for.  Gold¬ 
smith  £t  al. (1974)  had  a  special  Mercator  projection  constructed  for 
the  Virginia  shelf  to  minimize  distortion.  Robinson  (1969)  and  Green- 
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Figure  5.  Graphical  representation  of  the  1st-,  2nd-, 
and  3rd-order  charts  and  their  relation  to  each  other. 
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Table  1.  Hydrographic  Chart  Information 


Charts 

Area 

Year 

Scale 

1ST 

ORDER 

11520 

S.C. 

1975 

1:432,720 

11480 

coast 

1975 

1:449,659 

2ND 

ORDER 

1238 

Cape  Romain 

1973 

1:80,000 

to  Folly  Is. 

1973 

1:80,000 

3RD 

ORDER 

H-8779 

Price 

1963 

1:20,000 

H-4179 

Inlet 

1921 

1:20,000 

H-4180 

1921 

1:20,000 

Figure  6.  Maximum  accuracy  criteria  used  for  soundings  on 
hydrographic  charts  since  1860.  Modified  from  Sallenger 
(1975).  Refer  to  Table  2  for  detailed  explanation  of 
this  figure. 
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hood  (1964)  state  that  the  amount  of  distortion  decreases  closer  to 
the  equator  and  that  for  small  areas  (as  compared  to  the  whole  globe) 
may  be  negligible.  The  South  Carolina  coast  lies  between  32°N  and  34°N 
and  so  qualifies  as  a  small  area.  Therefore,  the  standard  Mercator  pro¬ 
jections  issued  by  NOAA  were  used  without  modification. 
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DATA  OUTPUT  -  WAVE-REFRACTION  DIAGRAMS 

Description 

A  complete  description  of  the  types  of  output,  including  sample 
listings  of  the  wave  parameters  printed  by  REFRAC,  can  be  found  in  Ap¬ 
pendix  LL.  Refraction  diagrams  referred  to  in  this  section  are  in  Ap¬ 
pendix  l.  (Arrows  refer  to  direction  of  sediment  transport.) 

Wave-refraction  diagrams  are  a  graphical  representation  of  the  bend¬ 
ing  of  wave  rays  as  the  waves  interact  with  the  continental  shelf.  This 
interaction  results  in  an  uneven  distribution  of  wave  characteristics 
along  the  shoreline,  as  inferred  by  zones  of  converging  and  diverging 
wave  orthogonals.  Converging  wave  orthogonals  depict  areas  of  increased 
wave  height  and  energy  commonly  correlative  with  areas  of  erosion.  Con¬ 
versely,  diverging  orthogonals  depict  areas  of  lower  wave  height  and 
energy  commonly  correlative  with  areas  of  deposition.  To  accurately 
delineate  these  areas  of  erosion  and  deposition  for  a  segment  of  the 
S.  C.  coastline,  a  set  of  first-,  second-,  and  third-order  wave-refrac¬ 
tion  diagrams  have  been  generated. 

A  first-order  diagram  shows  the  gross  energy  distribution  over  an 
extended  area.  Figure  1,  for  example,  shows  the  over-all  wave-refrac¬ 
tion  pattern  on  the  S.  C.  continental  shelf  for  waves  with  a  10-second 
period  from  the  east.  The  shelf  break  ranges  from  between  47  km  and  75 
km  offshore  at  the  75  to  100  fathom  contour.  Figure  2  is  a  second-order 
diagram  showing  the  energy  distribution  between  Bull  Island  and  the 
Isle  of  Palms.  This  figure  was  generated  from  the  first-order  diagram 
shown  in  Figure  1.  The  increase  in  detail  is  the  result  of  an  increase 
in  the  density  of  wave  rays  refracted  and  the  greater  bathymetric  de¬ 
tail  of  the  2nd-order  chart.  The  rays  were  started  in  deep  water 
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about  1  km  apart.  A  subset  of  the  rays  refracted  in  Figure  2  generated 
the  third-order  diagrams  in  Figures  6a  and  6b.  The  deep  water  spacing 
between  these  rays  is  the  same  as  in  Figure  2.  These  third-order  dia¬ 
grams  show  the  detailed  wave-energy  distribution  on  Price  Inlet’s  ebb- 
tidal  delta. 

The  bathymetry  of  Price  Inlet  was  surveyed  in  1921  (PI  -  1921)  and 
again  in  1963  (PI  -  1963) .  The  maximum  depth  of  both  surveys  is  between 
25  and  30  feet.  Incoming  waves  with  periods  of  6  seconds  or  greater 
have  already  felt  bottom  and  begun  to  refract  by  the  time  they  reach 
the  maximum  depth  of  the  above  surveys.  Therefore,  what  is  being  ob¬ 
served  on  these  charts' is  the  ’fine  tuning’  of  a  process  that  has  al¬ 
ready  begun  in  water  deeper  than  30  feet.  Since  the  same  second-order 
hydrographic  chart  generated  input  data  at  the  seaward  extent  of  both 
3rd-order  charts,  wave  refraction  diagrams  for  PI-1921  and  PI-1963  are 
similar.  Changes  of  the  inner  shelf  bathymetry  have  been  noted  by 
Goldsmith  et;  al.  (1975)  off  the  southern  Delmarva  Peninsula,  Virginia. 
Unfortunately,  a  smaller  scale  chart  for  1921  was  unobtainable. 

Only  2nd-  and  3rd-order  refraction  diagrams  will  be  described  in 
detail.  Figures  2-5  are  2nd-order  diagrams  of  waves  approaching  from 
the  east.  These  waves,  typical  of  winter  months,  are  generated  by 
extratropicaL  storms  passing  north  of  South  Carolina.  A  study  of 
these  figures  reveals  the  strong  effect  offshore  bathymetry  has  on 
an  easterly  wave  approach.  Waves  with  a  period  of  10  seconds  (Fig.  2) , 
after  interacting  witli  the  continental  shelf,  have  an  exceptionally 
nonuniform  distribution  of  wave  energy.  A  very  strong  zone  of  con¬ 
vergence  exists  at  Capers  and  Dewees  Inlets  and  Capers  Island.  A 
dramatic  decrease  in  wave  energy  is  observed  south  of  Dewees  Inlet 


and  north  of  Price  Inlet.  The  direction  of  sediment  transport  was  cal¬ 
culated  by  REFRAC  to  be  predominantly  to  the  south.  In  analyzing  these 
diagrams,  it  is  important  to  remember  the  variability  in  the  amount  of 
sand  being  transported  as  a  result  of  the  nonuniform  distribution  of 
wave  energy.  For  example,  in  Figure  2,  more  sand  is  capable  of  being 
transported  at  Capers  and  Dewees  Islands  and  their  adjoining  inlets  be¬ 
cause  of  the  increased  wave  energy  than  further  south  on  the  Isle  of 
Palms.  Therefore,  from  Figure  2  alone,  deposition  may  be  expected  on 
the  Isle  of  Palms  because  of  the  lower  wave  energy  and  the  resulting 
decrease  in  competency  of  the  longshore  currents. 

Figure  3  is  a  wave  refraction  pattern  for  8-second  waves  from  the 
east.  The  wave  energy,  as  compared  to  Figure  2,  is  more  evenly  dis¬ 
tributed.  A  strong  zone  of  convergence  exists  in  the  central  section 
of  the  Isle  of  Palms  and  the  northern  tip  of  Bull  Island.  Price  In¬ 
let  and  the  remainder  of  Bull  Island  seem  to  be  in  a  zone  of  reduced 
wave  energy  for  both  10-  and  8-second  waves.  The  predominant  direc¬ 
tion  of  sediment  transport  for  an  8-second  wave  is  to  the  south. 

Refraction  diagrams  for  waves  with  periods  of  6-  and  4-seconds 
(Figs.  4  and  5,  respectively)  have  a  progressively  more  uniform  wave 
energy  distribution  than  the  energy  distribution  for  8-  and  10-second 
waves.  This  results  from  longer  period  waves  beginning  to  refract  in 
deeper  water.  Furthermore,  as  a  result  of  the  different  depths  in 
which  waves  of  varying  periods  begin  to  refract,  the  orthogonals  of 
waves  with  longer  periods  tend  to  be  more  perpendicular  to  the  coast 
(compare  Figs.  2-5).  Figure  4  (period  (T)=6  seconds)  shows  a  zone 
of  convergence  on  the  shoals  north  of  Price  Inlet.  Sediment  trans¬ 
port  is  predominantly  to  the  south  except  on  these  northern  shoals 


wtiere  it  is  to  the  north.  In  Figure  5  (T-4  sec) ,  the  zone  of  conver¬ 
gence  has  shifted  to  the  southern  shoals  of  Price  Inlet  and  the  nor¬ 
thern  end  of  Capers  Island.  The  sediment  transport  is  predominantly  to 
the  south,  including  the  northern  shoals  of  Price  Inlet.  The  difference 
in  transport  direction  on  the  northern  shoals  of  Price  Inlet  is  ex¬ 
plained  by  the  more  perpendicular  orientation  of  the  wave  rays  with 
the  shore,  for  waves  with  a  longer  period.  This  implies  that  shorter 
period  waves  may  actually  enhance  sediment  bypassing  at  Price  Inlet. 

Third-order  refraction  diagrams  show  the  energy  distribution  and 
resulting  sediment  transport  at  Price  Inlet  in  greater  detail.  The 
direction  of  sediment  transport  at  the  northern  shoals  of  Price  Inlet 
is  to  the  south  in  Figures  9a  and  9b  (T=4  sec);  whereas,  for  10-,  8-, 
and  6-second  waves  (Figs.  6-8)  transport  is  to  the  north.  Comparing 
the  direction  of  sediment  transport  in  Figs  6a-9a  and  6b-9b  reinforces 
the  above  observation  that  sediment  bypassing  at  the  distal  margin  of 
the  ebb-tidal  delta  of  Price  Inlet  may  be  enhanced  by  waves  with  a 
shorter  period.  Further  analysis  of  these  diagrams  reveals  the  exis¬ 
tence  of  sediment  transport  reversals  south  of  the  inlet  (Figs. 7a  and 
9b).  Sediment  transport  reversals  are  caused  by  the  refraction  of 
waves  around  an  ebb-tidal  delta.  Finley  (1976)  documented  a  similar 
transport  reversal  at  North  Inlet,  South  Carolina.  Another  observa¬ 
tion,  seen  in  Fig.  8a  (T=6  sec),  is  the  strong  convergence  of  wave 
rays  at  the  seaward  extent  of  the  northern  marginal  flood  channel. 

IJuring  the  summer  months,  storms  generate  waves  from  the  south¬ 
east  and  south.  Second-order  diagrams  of  waves  from  the  southeast 
are  shown  in  Figs.  10-13.  The  influence  of  the  offshore  bathymetry 
on  these  waves  is  evidenced  by  the  alternate  zones  of  converging  and 
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diverging  wave  rays.  Contrary  to  waves  from  the  east,  the  zones  of 
wave  energy  concentration  do  not  shift  with  a  change  in  period.  For 
waves  from  the  southeast  with  periods  of  10- ,  8-,  and  6-seconds  (Figs.  10- 
12) ,  strong  zones  of  convergence  occur  at  the  southern  end  of  the  Isle 
of  Palms  (Breach  Inlet)  and  the  northern  section  of  Capers  and  Bull  Is¬ 
lands.  A  weaker  concentration  of  wave  energy  is  seen  at  Dewees  Inlet. 

The  only  difference  between  periods  is  the  degree  of  concentration; 
smaller  period  waves  are  usually  not  as  focused.  Zones  of  reduced  en¬ 
ergy,  such  as  Price  Inlet,  exist  between  the  above  zones  of  concentra¬ 
tion.  Waves  with  a  4-second  period  show  a  comparatively  even  wave- 
energy  distribution. 

The  direction  of  sediment  transport  was  determined  to  be  the  same 
for  all  4  periods.  Transport  is  predominantly  to  the  south  with  an 
indication  of  transport  to  the  north  at  the  northern  shoals  of  Price 
Inlet  and  Capers  Inlet.  As  seen  from  the  third-order  refraction  dia¬ 
grams  (Figs.  14-17),  the  predominant  sediment  transport  is  also  to 
the  south,  but  with  one  important  difference.  A  sediment  transport 
reversal  to  the  north  exists  at  the  southern  shoals  of  Price  Inlet 
in  all  the  diagrams  except  Figure  17b,  which  is  for  a  wave  with  a 
period  of  4-seconds. 

Waves  from  the  south  are  not  affected  as  much  by  the  offshore 
bathymetry  as  waves  from  the  east  and  southeast.  As  a  result,  wave 
energy  is  more  evenly  distributed  (Figs.  18-21).  Fig.  18  (T=10  sec) 
shows  zones  of  minor  concentration  at  the  southern  and  northern  ends 
of  the  Isle  of  Palms.  Bull  Island  experiences  a  slight  increase  in 
wave  energy  at  the  southern  and  central  sections  of  the  island.  A 
more  even  energy  distribution  is  seen  in  Fig.  19  (T*8  sec),  except 
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for  a  strong  concentration  of  wave  rays  on  the  mid-section  of  the  Isle 
of  Palms.  This  zone  is  surrounded  on  either  side  by  a  reduction  in 
wave  energy.  Figs.  20  (T=6  sec)  and  21  (T=4  sec)  have  a  fairly  uniform 
energy  distribution.  The  only  difference  is  seen  in  the  direction  of 
sediment  transport.  Figs.  18-21  show  a  dominant  transport  to  the  north. 
However,  for  waves  witti  a  10- ,  8-,  and  6-second  period,  the  sediment 
transport  on  the  southern  shoals  of  Price  Inlet  is  to  the  south  (Figs. 
18-20).  Transport  is  to  the  north  for  a  wave  with  a  4-second  period 
(Fig.  21).  By  comparing  the  above  4  figures  with  the  more  detailed 
third-order  diagrams.  Figs.  22-25,  it  is  apparent  that  the  lack  of  de¬ 
tail  in  Fig.  21  accounts  for  the  discrepancy  in  the  transport  direc¬ 
tions.  Fig.  25a  and  b  shows  a  sediment  transport  direction  to  the 
south  for  4-second  waves  on  the  southern  shoals  of  Price  Inlet. 
Discussion 

Wave- refraction  analysis  provides  an  efficient  method  for  inter¬ 
preting  coastal  processes.  When  used  in  conjunction  with  field  data, 
refraction  diagrams  aid  in  the  understanding  of  coastal  geomorphology. 
Stephen  £t  al.  (1975)  measured  the  rates  of  shoreline  change  from  ver¬ 
tical  aerial  photographs  between  the  years  1939  and  1973  for  Charles¬ 
ton  County,  South  Carolina.  The  resulting  classification  scheme  has 
four  categories:  areas  of  long-term  erosion,  long-term  accretion,  un¬ 
stable  areas  and  stable  areas.  Refer  to  Fig.  7  for  an  explanation 
and  the  geographical  location  of  these  categories.  Included  in  this 
figure  are  the  locations  of  the  zones  of  wave  energy  concentration 
and  reduction  described  in  the  previous  section.  A  careful  analysis 
of  Fig.  7  reveals  a  relation  between  the  zones  of  high  and  low  wave 
energy,  as  depicted  by  wave- re fraction  diagrams,  and  the  coastal  geo- 
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Figure  7.  Classification  scheme  of  the  shoreline  changes  for 
the  S.  C.  coast,  based  on  aerial  photographs  from  1939  to  1973. 
Modified  from  Stephen  c^t  al.  (1975). 

long-term  erosion:  areas  which  have  undergone  relatively 
continuous  erosion  over  the  study  interval, 
long-term  accretion:  areas  which  have  undergone  relatively 
continuous  deposition  over  the  study  interval, 
unstable:  areas  with  fluctuations  greater  than  50  ft 
over  the  study  interval. 

stable:  areas  with  fluctuations  of  less  than  50  ft  over 
the  study  period. 

Superimposed  on  these  shoreline  changes  are  the  various  degrees 
of  wave  ray  concentrations  observed  from  the  wave  refraction  dia¬ 
grams. 


exceptionally  high  concentration 
high  concentration 
uniform  distribution  of  wave  rays 
low  concentration  or  diverging  rays 
exceptionally  low  concentration. 
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morphology. 

There  is  a  dynamic  interaction  between  the  tidal  currents  and  long¬ 
shore  currents  generated  by  waves  which  affects  the  orientation  and  mor¬ 
phology  of  ebb-tidal  deltas  (Oertel,  1975;  Oertel  and  Howard,  1972;  Hub¬ 
bard,  1977).  Tide-dominated  inlets  have  well-developed  shoals  forming 
a  high  oblique  angle  with  the  shoreline.  Waves  have  a  tendency  to 
straighten  a  coastline  by  focusing  energy  on  headlands.  As  a  result, 
wave  dominated  ebb-tidal  deltas  have  a  relatively  parallel  orientation 
with  the  shoreline.  Coleman  (1976)  cites  similar  behavior  in  river 
deltas,  only  the  interaction  is  between  sediments  supplied  by  the  river 
and  the  wave  regime. 

Analysis  of  Fig.  7  shows  the  orientation  of  Price  Inlet's  ebb-tid¬ 
al  delta  to  be  more  normal  to  the  shoreline  than  the  orientation  of 
the  Capers-Dewees  shoal  complex  further  to  the  south.  One  explanation 
for  this  is  obtained  by  comparing  the  refraction  patterns  on  these  in¬ 
lets.  Price  Inlet  exists  essentially  in  a  shadow  zone  of  reduced  wave 
energy  for  waves  from  the  east,  which  is  the  direction  of  dominant 
energy  flux.  Therefore,  tidal  currents  play  a  more  active  role  in 
the  development  of  the  ebb-tidal  delta,  resulting  in  an  almost  per¬ 
pendicular  orientation  to  the  shoreline.  The  Capers-Dewees  complex, 
however,  lies  in  a  zone  of  exceptionally  high  wave  energy  concentra¬ 
tion  for  a  10-second  wave  from  the  east.  As  a  result,  waves  tend 
to  push  the  sand  up  against  the  shore  causing  a  more  oblique  or  para¬ 
llel  orientation  wich  the  shoreline. 

Besides  the  oblique  orientation  of  Dewees  Inlet,  a  large  down- 
drift  offset  is  exhibited.  The  ebb-tidal  delta  will  act  as  a  bar¬ 
rier  against  wave  attack,  but  the  downdrift  offset  is  probably  more 
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related  to  a  sudden  reduction  in  southward  directed  longshore  drift  at 


I 


i 
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the  Isle  of  Palms.  This  reduction  in  longshore  drift,  caused  by  the 
nonuniform  energy  distribution,  allows  the  deposition  of  sediment  on 
the  northern  end  of  the  Isle  of  Palms.  The  area  is  unstable  because 
of  its  proximity  to  Dewees  Inlet. 

Figure  7  shows  long-term  erosion  further  north  at  Bull  Island. 
North  Bull  Island  has  no  protecting  shoals  and  is  therefore  open  to 
wave  attack,  especially  from  an  8-  or  6-second  wave  from  the  east. 
Figs.  10-13  (in  Appendix  I)  show  a  concentration  of  wave  energy  on 
Bull  Island  for  all  waves  from  the  southeast.  This  concentration 
of  wave  rays  results  from  the  interaction  of  waves  and  the  section 
of  the  shoal  retreat  massif  off  the  northern  end  of  Bull  Island 
(30  foot  contour  and  shallower).  Also,  there  is  no  apparent  sedi¬ 
ment  source  to  allay  the  erosional  trend  on  Bull  Island. 

Capers  Island  is  also  undergoing  long-term  erosion.  As  seen 
from  Fig.  7,  strong  concentrations  of  wave  energy  are  noted  for 
waves  from  the  east  with  a  wave  period  of  10-seconds  and  waves 
from  the  southeast  with  a  period  of  10-,  8-,  or  6-seconds. 

A  maximum  accretion  of  400  feet  since  1941  has  occurred  on  the 
southern  end  (spit)  of  the  Isle  of  Palms.  This  area  is  one  of  the 
few  areas  on  the  S.  C.  coast  categorized  as  long-term  accretion. 

An  analysis  of  Fig.  7  reveals  this  section  of  the  Isle  of  Palms  to 
be  in  a  zone  of  reduced  wave  energy  for  most  wave  conditions. 
Therefore,  deposition  and  spit  growth  would  be  enhanced. 

A  comparison  of  Figs.  1-5,  10-13,  18-21  (in  Appendix  I)  and 
Fig.  7  suggests  that  longer  period  waves  have  more  effect  on  the 
shape  of  the  coast  than  shorter  period  waves.  Because  of  the  more 


dramatic  concentration  or  divergence  of  wave  rays  for  longer  period 
waves,  sections  of  the  coast  will  accrete  or  erode  at  different  rela¬ 
tive  rates,  therefore  influencing  the  gross  morphology  of  the  coast 
to  a  larger  degree  than  shorter  period  waves.  Shorter  period  waves 
will  have  relatively  even  rates  of  erosion  or  deposition  because  of 
their  more  uniform  wave  energy  distribution.  The  accretion  on  the 
Isle  of  Palms,  the  downdrift  offset  at  Dewees  Inlet,  and  the  orien¬ 
tation  of  Dewees  Inlet  as  compared  to  Price  Inlet  support  this  ob¬ 
servation.  The  possibility  that  longer  period  waves  have  more  effect, 
on  coastal  geomorphology  implies  that  the  deeper  offshore  bathymetry 
(greater  than  approximately  48')  affects  the  coast  more  than  the  near¬ 
shore  bathymetry.  To  determine  the  validity  of  this  observation, 
more  precise  data  on  the  frequency  of  various  wave  conditions  is 
necessary. 

To  gain  an  understanding  of  wave  refraction  and  its  possible  ef¬ 
fects  on  the  morphology  of  an  ebb-tidal  delta,  refraction  diagrams 
with  greater  detail  were  generated  for  Price  Inlet.  Price  Inlet  is 
an  example  of  a  model  developed  by  Hayes  (1975)  on  ebb-tidal  deltas 
(Fig.  8).  Basically,  waves  will  bend  toward  the  inlet.  The  angle 
the  wave  makes  with  the  shoreline  is  of  critical  importance  in  deter¬ 
mining  the  direction  of  sediment  transport.  Waves  from  the  east  and 
southeast  have  a  predominant  sediment  transport  to  the  south.  A 
careful  examination  of  Figs.  6-9,  14-17  (in  Appendix  I)  shows  sediment 
transport  reversals  south  of  the  ebb-tidal  delta.  Transport  rever¬ 
sals  result  from  the  sharp  bending  of  wave  rays  around  the  delta 
and  provide  a  means  for  both  reintroducing  sand  into  the  inlet, 
through  the  marginal  flood  channels,  and  building  the  beach  and 
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Figure  8.  General  ebb-tidal  delta  model 
(from  Hayes,  1975). 
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shoals  south  of  the  Inlet.  This  mechanism  was  proposed  by  Hayes  et 
al.  (1970)  as  being  responsible  for  downdrift  offsets. 

The  angle  at  which  waves  strike  the  ebb-tidal  delta  may  enhance  or 
inhibit  sediment  bypassing  around  the  distal  shoals  of  the  delta.  Sand 
may  be  transported  landward  onto  the  swash  platform  or  along  the  shoal 
margins  depending  on  the  tide  level.  At  low  tide,  the  shoals  are  sub¬ 
aerial.  As  a  result,  sand  transport  will  ocv.ur  along  the  shoal  margins. 
For  waves  from  the  east,  the  net  longshore  transport  is  to  the  south. 
However,  an  analysis  of  Fig.  6-9  (in  Appendix  I)  reveals  that  on  the 
northern  shoals  of  the  delta,  sediment  transport  is  to  the  north  for 
waves  with  a  period  of  10- ,  8-,  and  6-seconds.  Therefore,  sand  trans¬ 
port  to  the  south  around  the  distal  end  of  the  delta  is  inhibited. 

For  a  4-second  wave,  sediment  transport  on  the  northern  shoals  is  to 
the  south  and,  therefore,  may  actually  enhance  sediment  bypassing. 

This  is  explained  by  the  more  oblique  approach  of  a  4-second  wave  to 
the  shoreline  than  the  longer  period  waves. 

Many  of  the  ebb-tidal  deltas  on  the  S.  C.  coast  exhibit  channels 
which  run  parallel  to  the  shore.  During  early  flood  when  the  swash 
bars  are  exposed,  these  marginal  flood  channels  provide  an  avenue 
for  water  flow  and  sediment  transport  into  the  inlet.  In  the  vicinity 
of  the  marginal  flood  channels.  Figures  6-9,  14-17,  22-25  (in  Appen¬ 
dix  I)  show  converging  sediment  transport  directions.  This  conver¬ 
gence  may  create  a  hydraulic  head  at  the  seaward  extent  of  the  mar¬ 
ginal  flood  channel,  thereby  enhancing  water  flow  and  sediment  trans¬ 
port  through  the  channels.  The  larger  size  of  the  northern  marginal 
flood  channel  probably  reflects  the  effects  of  the  direction  of  domi¬ 
nant  wave  energy  flux  from  the  northeast  and  east.  The  wide  funnel 
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shape  of  this  channel  may  result  from  the  strong  convergence  in  wave 
rays  noted  in  Fig.  8  (in  Appendix  I).  This  convergence  of  rays  would 
cause  both  an  increase  in  wave  setup  and  energy  resulting  in  increased 
sediment  transport  and  erosion  at  the  seaward  extent  of  the  flood  chan¬ 
nel. 

In  nature,  a  spectrum  of  wave  conditions  strike  the  shore  at  any 
one  time.  Each  component  of  this  spectrum  (different  period  and  di¬ 
rection)  has  its  own  distribution  of  wave  height  and  energy.  This 
nonuniform  distribution  of  wave  characteristics,  caused  by  wave  re¬ 
fraction,  may  partially  explain  the  intermittent  character  of  the 
flood-directed  stresses  noted  by  Huntley  and  Nummedal  (pers.  comm.) 
in  the  marginal  flood  channels.  These  authors  made  a  series  of  along- 
channel  and  cross-channel  velocity  measurements  in  both  channels.  An 
analysis  of  the  variability  of  these  two  directional  velocities  pro¬ 
vided  a  measure  of  the  radiation  stress  available  to  drive  longshore 
currents  in  the  channels.  These  stresses  were  intermittent  and  re¬ 
lated  to  obliquely  incident  waves. 
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CONCLUSIONS 


The  offshore  bathymetry  creates  zones  of  converging  and  diverging 
wave  rays  by  interacting  with  waves  of  different  periods  and  directions. 
The  importance  of  this  interaction  of  coastal  geomorphology  is  seen  in 
the  following  observations  from  wave  refraction  diagrams  of  the  S.  C. 
coast  from  Bull  Island  to  the  Isle  of  Palms: 

1.  Zones  of  high  and  low  wave  energy,  as  depicted  by  converging 
and  diverging  wave  rays,  tend  to  have  corresponding  zones  of  erosion 
and  deposition.  Both  Bull  Island  and  Capers  Island  which  are  under¬ 
going  long-term  erosion,  are  located  in  areas  of  wave  convergence. 

Bull  Island  lacks  an  apparent  sediment  source  which  may  partially  ac¬ 
count  for  the  erosion.  Further  south,  the  southern  tip  of  the  Isle 
of  Palms  has  undergone  extensive  deposition.  This  section  of  the 
Isle  of  Palms  corresponds  to  a  zone  of  reduced  wave  energy. 

2.  The  large  downdrift  offset  of  Dewees  Inlet  appears  to  be  re¬ 
lated  to  a  sudden  reduction  in  the  southward-directed  longshore  drift 
at  the  Isle  of  Palms  as  compared  to  Dewees  Island  and  other  barriers 
further  north. 

3.  The  orientation  of  Dewees  Inlet  and  Price  Inlet  is  related  to 
the  balance  between  wave  energy  and  tidal  energy.  Price  Inlet,  be¬ 
ing  in  a  zone  of  wave  energy  reduction,  is  relatively  dominated  by 
the  tidal  current  resulting  in  a  normal  orientation  to  the  shoreline. 
Dewees  Inlet  exhibits  an  oblique  orientation  with  the  shoreline.  This 
orientation  is  strongly  related  to  the  exceptionally  high  concentra¬ 
tion  of  wave  energy  received  from  a  10-second  wave  from  the  east. 

4.  This  study  suggests  that  longer  period  waves  and,  therefore, 
the  deeper  offshore  bathymetry  may  have  a  stronger  effect  on  the 
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coastal  geomorphology  than  shorter  period  waves  and  the  nearshore 
bathymetry.  Evidence  for  this  is  seen  by  the  accretion  on  the  Isle  of 
Palms,  the  downdrift  offset  at  Dewees  In let,  and  the  shore-parallel 
orientation  of  Dewees  Inlet. 

5.  Sediment  transport  reversals  are  seen  at  Price  Inlet  for  waves 
from  the  east  and  southeast.  The  calculated  direction  of  transport 
implies  that  these  reversals  provide  a  means  for  reintroducing  sand 
into  the  inlet  and  depositing  sand  on  the  beach,  and  shoals  south  of 
the  inlet. 

6.  Sediment  bypassing  around  the  distal  portion  of  the  ebb-tidal 
delta  at  Price  Inlet  may  be  enhanced  by  shorter  period  waves  from  the 
east. 

7.  In  relation  to  marginal  flood  channels  at  Price  Inlet,  the  fol¬ 
lowing  observations  were  made: 

a)  Hydraulic  heads  resulting  from  the  convergence  of  two  dif¬ 
ferent  directions  of  transport  may  enhance  water  flow  and 
sediment  transport  through  the  flood  channels. 

b)  A  wave  with  a  period  of  6-seconds  from  the  east  may  account 
for  the  funnel  shape  of  the  seaward  extent  of  the  northern 
marginal  flood  channel  because  of  the  strong  zone  of  conver¬ 
gence  in  the  vicinity. 

c)  The  pulsating  flood  currents  recorded  by  Huntley  and  Numme- 
dal  (pers.  comm.)  may  partially  be  explained  by  the  variations 

in  wave  height  produced  at  a  point  due  to  the  simultaneous 
refraction  of  several  different  wave  conditions. 
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APPENDIX  II 
REF RAC — User's  Guide 


PROGRAM  DESCRIPTION 


Module  Description 

REFRAC  is  divided  Into  a  MAIN  routine  and  several  subroutines: 
RAYCON,  REFRAC,  CURVE,  DEPTH,  HEIGHT,  REFDIA,  ENERGY,  BOTTOM  and 
SORT1.  The  subroutine  ENERGY  is  similar  to  the  subroutine  used  by 
Colonell  and  Goldsmith  (personal  communication)  for  calculations  of 
wave  power  and  energy. 

MAIN  Routine.  -  This  routine  reads  the  input  wave  parameters 
such  as  initial  wave  height,  period,  direction  and  starting  grid 
coordinates.  Calculated  refraction  coefficients  need  be  supplied  if 
the  wave  starts  in  shallow  water.  Controlling  program  parameters  such 
as  window  specifications,  number  of  wave  rays  and  so  on  are  also  read. 
Output  includes  a  title  page  listing  the  input  parameters  for  quick  re¬ 
ferral  and  comparisons  between  sets  of  rays.  The  shoreline,  depth 
contours  and  plot  title  information  are  plotted  at  this  time.  Cal¬ 
culations  of  deep  water  wave  celerity  and  length  and  the  depth  at 
which  refraction  will  begin  are  performed.  Further  calculations  in¬ 
clude  printing  and  plotting  intervals  based  on  a  time  step.  Time 
was  chosen  as  the  independent  variable  in  the  governing  equations 
because  in  deep  water  where  effects  of  refraction  are  small,  the 
celerity  and,  therefore  the  distance,  between  two  points  decreases. 

This  coincides  with  the  area  in  which  the  greatest  detail  is  desired. 
The  first  subroutine  called  in  this  program  is  BOTTOM  which  sets  up 
the  depth  grid  for  later  calculations.  Nested  loops  control  the  num¬ 
ber  of  sets  of  rays  and  the  number  of  rays  per  set  to  be  analyzed. 

In  the  inner  loop,  control  is  passed  to  subroutine  RAYCON  which  cal- 
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culates  each  ray's  path  across  the  depth  grid. 

Subroutine  KAYCON.  -  RAYCON  controls  each  individual  ray  as  it  pro¬ 
gresses  across  the  grid,  initially,  the  wave  is  advanced  one  step  and 
the  depth  calculated  (from  DEPTH  subroutine)  at  this  point.  If  the 
depth  is  greater  than  one-half  the  deep  water  wave  length,  deep  water 
conditions  exist  signifying  no  changes  in  the  wave.  Printing  and 
plotting  options  are  then  checked.  The  above  steps  are  repeated  un¬ 
til  transitional  water  depths  are  reached,  that  is,  until  the  wave 
feels  bottom.  Wave  height  and  boundary  positions  are  continuously 
checked  for  breaking  wave  or  out  of  bound  conditions.  As  the  wave 
progresses  through  transitional  water  depths,  the  ray  begins  to  re¬ 
fract.  RAYCON  repeatedly  calls  CURVE,  REFRAC,  and  HEIGHT  to  calculate 
the  curvature,  step  length  and  wave  height  resulting  from  increased 
refraction  and  shoaling.  The  ray  may  be  stopped  for  a  variety  of 
reasons:  there  is  no  convergence  in  the  calculation  for  curvature, 
the  wave  breaks,  the  ray  reaches  one  of  the  boundaries,  the  maximum 
number  of  points  calculated  has  been  exceeded,  or  the  incremental  dis¬ 
tance  between  steps  is  less  than  the  minimum  specified.  As  mentioned 
previously,  RAYCON  controls  the  printing  and  plotting  for  each  ray. 

At  user-determined  intervals,  a  tick  mark  is  plotted  on  the  wave  ray, 
and  corresponding  wave  information,  denoted  by  an  asterisk,  is  printed. 
Subroutine  ENERGY  is  called  before  each  line  is  printed  in  order  to 
calculate  wave  power  and  energy  at  that  point.  The  boundaries  for  a 
more  detailed  window  are  calculated  if  a  window  has  been  specified  in 
the  input  parameters.  As  a  ray  crosses  these  boundaries,  values  are 
generated  and  stored  on  disk  to  be  used  as  input  data  into  a  later, 
more  detailed  refraction  analysis. 


Subroutine  REFRAC.  -  This  routine  and  the  following  subroutines  are 
called  repeatedly  after  the  wave  ray  has  entered  transitional  water 
depths.  REFRAC  calculates  the  step  length  which  is  a  function  of  the 
wave  celerity  and  time.  If  the  step  length  is  greater  than  the  minimum 
value  specified  in  the  input,  the  curvature  (through  CURVE  subroutine) 
and  next  X,Y  coordinates  are  calculated.  The  initial  curvature  is  used 
to  project  the  next  point  on  the  ray.  A  new  curvature  at  this  point  is 
calculated  and  compared  to  the  original  curvature,  thus  giving  a  better 
estimate  of  the  position  for  this  point.  This  process  is  repeated  un¬ 
til  the  difference  between  two  successive  calculations  is  acceptably 
small;  then  the  new  point  is  considered  fixed.  Two  conditions  of  insta¬ 
bility  can  arise.  One  condition  results  in  the  curvature  being  averaged 
between  solutions;  a  message  is  printed  to  this  effect.  The  other  con¬ 
dition  results  from  the  failure  of  the  calculations  to  converge,  in 
which  case  the  ray  is  stopped. 

Subroutine  CURVE.  -  Based  on  the  depth  of  the  water,  one  of  two 
equations  is  used  to  calculate  the  wave  celerity.  Five  one- thousandths 
(.005)  of  the  wavelength  is  the  program's  boundary  between  shallow 
water  and  intermediate  depths.  The  curvature  of  the  ray  is  then  cal¬ 
culated.  The  curvature  is  a  function  of  the  ratio  between  local  speed 
and  deep  water  wave  speed,  the  depth  and  several  depth-related  coeffi¬ 
cients,  plus  the  azimuth  of  the  ray. 

Subroutine  DEPTH.  -  A  second  degree  polynomial  fitted  locally 
(instead  of  on  the  complete  depth  grid)  provides  an  accurate  method 
for  interpolating  depths  at  intermediate  points  (Dobson,  1967).  The 

2  2 

DEPTH  =  e^  +  e2X  +  e^y  +  e^x  +  e^xy  +e^y 


equation  is: 


The  local  grid  system  consists  of  12  points  (Fig.  1)  which  are  filtered 
through  a  matrix  of  weighting  terms  derived  from  the  relative  positions 
of  the  data  points.  This  local  grid  system  traces  the  ray's  path  shore¬ 
ward  across  the  depth  grid.  The  six  coefficients  in  the  equation  result 
from  a  summing  f  the  different  filtered  depths  and  remain  constant  with 
in  the  local  coordinate  system.  Calculation  of  the  curvature  and  refrac 
tion  coefficient  for  a  point  are  based  on  these  six  coefficients. 

Subroutine  HEIGHT.  -  This  subroutine  calculates  the  shoaling  and 
refraction  coefficients,  which  are  used  in  the  computation  of  the  local 
wave  height.  The  shoaling  coefficient  is  equal  to  the  ratio  of  the 
group  velocity  in  deep  water  to  local  group  velocity  and  therefore  de¬ 
pendent  on  the  water  depth.  The  refraction  coefficient  is  determined 
by  the  change  in  separation  of  wave  rays.  Dobson  (1967)  created  a  sep¬ 
aration  factor  which  is  a  function  of  wave  celerity,  azimuth  of  the  ray 
and  the  six  depth  coefficients.  As  a  result,  the  calculation  of  the 
refraction  coefficient  for  the  present  time  step  is  dependent  on  the 
refraction  coefficient  of  the  previous  time  step. 

Subroutine  REFDIA.  -  This  subroutine  plots  the  X,  Y  coordinates  of 
the  ray  as  it  progresses  shoreward.  The  grid  coordinates  are  multiplied 
by  the  length  of  a  grid  side  to  convert  them  into  plot  inches  before 
plotting.  The  data  input  variable  BOUND  is  subtracted  from  the  deep-end 
of  the  y-axis  if  this  part  of  the  plot  is  not  desired. 

Subroutine  ENERGY.  -  The  equations  used  in  this  routine  and  their 
derivations  may  be  found  in  the  Shore  Protection  Manual  (Coastal  En¬ 
gineering  Research  Center,  1973).  As  part  of  each  printed  line  of  out¬ 
put,  the  total  wave  power  and  energy,  the  longshore  power  and  direc¬ 
tion  (left  or  right)  are  calculated  and  printed.  Longshore  transport 
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rates  based  on  wave  energy  flux  and  Calvin's  (1972)  gross  longshore 
transport  rates  based  on  breaker  height  are  calculated  at  the  breaker 
points.  These  transport  rates  are  in  cubic  yards  per  year  and  the  wave 
energy  flux  transport  rate  is  further  converted  to  cubic  meters  per 
year. 

Subroutines  BOTTOM  and  S0RT1.  -  BOTTOM  creates  an  evenly-spaced 
depth  grid  from  a  digitized  bathymetric  map.  A  maximum  of  20,000  grid 
points  may  be  generated  depending  on  the  spacing  interval  between  grid 
lines  which  is  user  determined.  The  digitized  data  includes  the  X,Y 
coordinates  of  the  map  perimeter  and  the  bathymetric  contours  and 
their  associated  depths.  Up  to  100  different  contours  can  be  digitized. 

The  digitized  data  is  read  and  converted  into  a  numeric  or  charac¬ 
ter  format  for  the  program.  As  the  data  is  read,  two  temporary  disk 
files  are  created.  These  files  store  the  X,Y  coordinates  for  the  shore¬ 
line  and  contours  and  will  be  plotted  when  program  control  returns  to  the 
MAIN  routine.  The  points  on  the  perimeter  are  read  first,  and  the  in¬ 
cremental  distances  and  depths  calculated.  The  depths  are  computed 
from  the  straight-line  slope  formula.  These  calculations  occur  be¬ 
tween  two  corners  of  the  map  unles.s  contours  cross  this  section  of  the 
perimeter,  in  which  case  calculations  occur  between  the  corner  and  the 
first  contour;  between  contours  (if  more  than  one  crossed  the  perimeter); 
and  the  last  contour  and  far  corner.  Upon  completion  of  the  perimeter, 
the  coordinates  of  the  contours  are  read  and  converted  into  the  proper 
format.  After  the  digitized  tape  is  read,  the  tape  is  rewound  and  the 
SORT1  routine  called. 

Subroutine  SORT1  written  in  Assember  language,  links  a  system 
utility  sort  to  REFRAC.  The  depths  and  their  coordinates  are  sorted 


1 1—7 


to  build  the  desired  depth  grid.  Note  that,  at  this  point,  the  only 
depths  are  the  incremental  depths  along  the  perimeter  and  the  depths 
located  at  the  contours.  After  the  sort  is  completed,  BOTTOM  is  re¬ 
entered  and  the  remaining  depths  calculated. 

Grid  depths  are  calculated  parallel  to  the  y-axis,  one  grid  line  at  a 
time.  The  distance  between  two  points  on  two  contours  or  a  contour 
and  perimeter  with  the  same  x-coordinate  is  completed.  The  depth  dis¬ 
tance  between  these  two  points  divided  by  their  distance  determines 
the  slope  which  is  used  to  calculate  the  incremental  depths  between 
the  two  points.  Because  of  the  method  of  interpolation  described,  the 
y-axis  should  be  drawn  as  perpendicular  to  the  contours  as  possible. 

Upon  completion  of  the  grid,  program  control  returns  to  the  MAIN 
routine  and  refraction  calculations  as  described  in  the  previous  sub¬ 
routines  begin.  Fig.  2  is  a  generalized  diagram  of  the  flow  of  control 
in  the  program  REFRAC. 

Output  Description 

REFRAC  has  two  types  of  output  to  describe  ray  refraction.  The 
first,  printed  information  includes  a  title  page  per  set  of  rays. 

Title  of  the  area  being  analyzed  and  date  of  analysis  are  printed  on 
top  of  the  page.  An  echo  print  of  all  input  parameters  allows  a  quick 
referral  when  comparing  sets  of  rays.  The  date  of  the  hydrographic 
chart  and  its  units  of  measurement  are  also  printed. 

Information  on  individual  rays  is  printed  next.  At  the  start  of  a 
ray,  a  variable  GKINC  is  calculated  and  printed.  GRINC  is  the  deep 
water  step  length  expressed  as  a  fraction  of  the  grid  square  and  should 
equal  approximately  .5.  Calculation  of  GRINC  is  further  discussed  in 
the  section  on  PROGRAMMING  CONSIDERATIONS  (Determination  of  Input  Param- 
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Figure  2.  Program  flow 
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meters).  Points  along  the  ray  are  printed  at  a  user  specified  interval. 
Asterisks  printed  beside  the  point  number  correspond  with  tick  marks 
plotted  perpendicular  to  the  ray  orthogonal.  Point  information  includes: 
location  of  point  on  the  grid,  azimuth  of  the  ray,  water  depth,  length, 
speed  and  height  of  wave,  refraction  and  shoaling  coefficients,  bottom 
and  mid  fluid  particle  velocities,  total  wave  energy  and  power,  the 
longshore  component  of  power  and  direction  (left  or  right  as  the  obser¬ 
ver  faces  the  ocean).  If  the  wave  breaks,  the  longshore  energy  flux 
factor  and  longshore  transport  rates  are  printed.  For  convenience, 
the  transport  rates  are  converted  into  metric  units.  Also,  the  gross 
longshore  transport  rate  as  derived  by  Galvin  (1972)  is  printed. 

The  second  type  of  output  is  a  plot  of  the  refracting  rays  as  they 
travel  shoreward.  Plot  output  includes  the  digitized  shoreline  and 
bathymetric  contours.  Wave  orthogonals  or  rays  (lines  parallel  to  the 
direction  of  wave  propagation)  are  plotted  with  perpendicular  tick 
marks.  Title  information  on  the  plot  includes:  the  period  and  deep 
water  wave  height  and  azimuth,  date  of  hydrographic  chart,  date  of 
plot,  water  level  height  of  chart  and  relative  water  depth  of  plot 
to  chart. 

Input  Description 

Card  input  formatting.  -  Input  values  are  of  three  types:  integer 
(1)  -  numeric  values  without  a  decimal  point;  floating  point  (F)  - 
real  numeric  values  which  include  a  decimal;  and  alphanumeric  (A)  - 
alphabetic  or  numeric.  A  field  refers  to  the  columns  of  a  punch  card 
where  a  value  is  to  be  found.  Integer  and  real  values  are  to  be  right- 
justified;  that  is,  punched  in  the  columns  furthest  to  the  right  in 
their  particular  field.  Alphanumeric  data  should  be  lef t-justif ied. 
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Blank  data 

fields  are 

interpreted  as 

zeros  by  the  program.  The  number 

in  parenthesis  following  the  variables  description  are  suggested  values 

Card  Sequence: 

Columns 

Variable 

Description 

Card  It 1 

1-40 

ITITLE  (A) 

Title  of  wave  refraction 
analysis 

Card  It 2 

1-5 

FAC  (F) 

Size  of  plot  desired  rela¬ 
tive  to  original  chart. 

Card  HZ 

1-4 

DATE1  (A) 

Year  of  hydrographic  chart 
soundings. 

7-18 

MEAS  (A) 

Units  in  which  depth  of 
chart  is  expressed. 

21-24 

W  (A) 

Water  level  chart  was 
sounded  at  (e.g.  MLLW,  MLW 
LSL,  MHW,  MHHW) . 

Card  HU 

1-5 

MI  (I) 

Number  of  grid  lines  in  X- 
direction. 

6-10 

MJ  (I) 

Number  of  grid  lines  in  Y- 
direction. 

11-15 

limnpt  (I) 

Maximum  number  of  steps  to 
be  calculated  for  a  ray 
(4000) . 

16-20 

NPRINT  (I) 

Printing  interval  (10). 

21-30 

GRID  (F) 

Number  of  feet  (from  hydro- 
graphic  chart)  that  equal 
the  length  of  a  grid  square 
side. 

31-40 

DCON  (F) 

Conversion  factor  for  depth 
values  to  feet. 

41-50 

DELTAS  (F) 

Minimum  step  length,  ex¬ 
pressed  as  a  fraction  of  a 
grid  square  (.002). 

Card  It  5 

1-10 

BOUND  (F) 

Number  of  inches  to  be  sub¬ 
tracted  from  deepwater  end 
of  plot. 

11-12 


Card  //6 


Card  //7 


11-20 

SCX  (F) 

Length  of  a  grid  square  in  inches. 

21-30 

XSG  (F) 

X-coordinate  of  lower  left  corner 
of  window.  Set  equal  to  MI  if  win¬ 
dow  is  not  desired. 

Note: 

If  window  is  not 
blank. 

desired,  remainder  of  card  is  left 

31-40 

YSG  (F) 

Y-coordinate  of  lower  left  corner 
of  window. 

41-50 

SCNV  (F) 

Magnification  of  window. 

51-60 

DGXL  (F) 

Length  of  window  (x-direction) 
expressed  in  plot  inches. 

61-70 

DGYL  (F) 

Height  of  window  (y-direction) 
expressed  in  plot  inches. 

1-5 

N0SETS  (I) 

Number  of  sets  of  rays  to  be 
processed. 

Note:  Repeat  the  following  cards  if  more  than  one  set  of 
rays  is  to  be  processed. 


1-5 

LPL0T  (I) 

Number  of  steps  between  plot 
points  on  a  ray  (10) . 

6-10 

NORAYS  <I) 

Number  of  rays  in  a  set  (Maxi¬ 
mum  400) . 

11-20 

T  (F) 

Wave  period  for  a  set  of  rays 
(seconds) . 

21-30 

HO  (F) 

Deep  water  wave  height  (feet). 

31-40 

SK  (F) 

Shoaling  coefficient  for  first 
time  step. 

41-50 

SKI  (F) 

Deep  water  shoaling  coefficient 
(usually  1). 

41-60 

THI  (F) 

Clockwise  angle  between  north 
on  map  and  y-axis  of  grid. 

61-70 

STAZ  (F) 

Deepwater  azimuth  of  a  set  of 
rays. 

71-80 

UNIT  (F) 

Timestep 

11-13 


1-5  ISP  (I)  Sets  print  option  on  check 

depth. 

*-l,  information  printed 
and  processing  of  ray  con¬ 
tinued. 

-0,  no  check  depth. 

-+1,  information  printed 
and  processing  of  ray 
stopped. 

6-10  LCK  (I)  Starting  position  for  set 

of  rays. 

=0,  rays  start  in  deep  water. 
=+l,  rays  start  in  shallow 
water. 

11-20  WPI  (F)  Number  of  wave  periods  be¬ 

tween  tick  marks  on  a  ray. 

21-30  CKDEP  (F)  Check  depth,  ray  information 

will  be  printed  for  step 
nearest  this  depth. 

31-40  DF  (F)  Factor  to  convert  depth 

values  from  one  water  sur¬ 
face  datum  to  another. 

The  factor  will  be  added 
to  the  depth  values.  If 
not  needed,  leave  blank. 

41-45  IWR  (I)  Specifies  if  starting  co¬ 

ordinates  for  a  set  of 
rays  is  the  same  as  a  pre¬ 
vious  set. 

•0,  new  wave  information 
to  be  read  off  of  cards 
or  disk. 

■+1,  wave  information  same 
as  previous  set. 

46-50  PD  (F)  Depth  at  which  printing  of 

ray  information  will  begin. 

If  left  blank,  PD  will  equal 
first  depth  value  on  grid. 

56-65  YFA  (F)  Factor  added  to  the  y-co- 

ordinate  of  the  starting 
point  of  a  ray  in  order  to 
decrease  the  amount  of  deep 
water  region  the  wave  tra¬ 
vels  over,  thereby  decreasing 
computation  time. 
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Card  it 9 
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IRED  (I)  If  IWR-0  and  the  information 
for  card  deck  #10  resides  on 
disk,  set  IRED  to  25,  other¬ 
wise  leave  blank.  IRED  equals 
25  implies  that  the  ray  infor¬ 
mation  was  generated  by  a  pre¬ 
vious  run  on  a  smaller  scale 
map  through  a  window  and  stored 


The  following  card  must  be  repeated  for  each  ray  in  a  set.  (Note: 
for  any  set  other  than  the  first  set,  if  IWR  “  1,  these  cards  are 
not  needed.  The  information  is  taken  from  the  previous  set.) 


Card  #10  1-10 

X  (F) 

X-coordinate  for  starting 
position  of  ray. 

11-20 

Y  (F) 

Y-coordinate  for  starting 
position  of  ray. 

Note:  The  following  variables  are 
in  shallow  water,  otherwise  leave 

needed  if  the  ray  is  starting 
blank. 

21-30 

AZIMTH  (F) 

Azimuth  of  a  ray. 

31-4C 

RK1  (F) 

Refraction  coefficient  of 
wave  ray  at  the  time  step 
previous  to  starting  loca¬ 
tion. 

41-50 

RK  (F) 

Refraction  coefficient  of 

wave  ray  at  starting  loca¬ 
tion. 


Digitized  depth  input.  -  This  section  describes  in  a  step  by  step 
format  how  to  prepare  a  hydrographic  chart  for  digitizing  and  the  digi¬ 
tizing  process. 

Step  1  -  Outline  the  area  of  interest.  Enclose  the  area  of  in¬ 
terest  in  a  rectangular  outline.  This  outline  is  the  perimeter  of  the 
grid  system  which  will  be  generated  from  the  digitized  contours.  The 
perimeter  should  be  oriented  so  that  the  wave  rays  travel  from  deep  to 
shallower  water.  In  other  words,  the  Y-axis  is  more  or  less  perpendicu¬ 
lar  to  the  shoreline  or  contours.  Figure  3  shows  the  perimeter  orienta- 
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Figure  3.  Hypothetical  hydrographic  chart  to  be  digitized. 
Contours  are  in  feet. 

X-axis  =  8  inches 
Y-axis  =5.5  inches 

If  SCX  -  .1 

MI  =  81 

MJ  =  56 

MI  x  MJ  =  4536  (Total  no.  of  grid  points) . 


tion  on  a  hypothetical  hydrographic  chart. 

Before  drawing  a  perimeter,  two  questions  need  to  be  considered. 
First,  what  are  the  directions  of  wave  approach?  If  the  waves  will  be 
starting  from  the  perimeter  sides,  extra  distance  along  the  X-axis  will 
be  needed  to  give  the  rays  time  to  refract  into  the  area  of  interest. 
Second,  how  many  grid  points  will  be  generated?  The  number  of  grid 
lines  crossing  the  X  and  Y  axis  (MI  and  MJ  respectively)  are  given  by 


Ml 

_  Length  of  respective  axis  in  inches  +  ^ 

^  Spacing  interval  between  grid  lines 

(SCX)  also  in  inches 

As  interpreted  from  the  above  equation,  the  origin  of  the  grid  is  (1,1). 
The  value  of  SCX  is  user  determined  and  will  depend  on  the  complexity 
of  the  bathymetry.  The  number  of  grid  points  generated  is  the  result  of 
multiplying  MI  and  MJ  together  and  must  be  less  than  the  maximum  of 
20,000.  The  values  of  MI,  MJ ,  and  SCX  will  be  entered  into  the  program 
as  card  input. 


STEP  2  -  What  to  digitize.  First  the  perimeter  and  the  points 
where  the  contours  cross  the  perimeter  are  digitized.  Next,  the  con¬ 
tours  are  digitized.  A  problem  arises  in  digitizing  the  shoreline. 

In  Figure  3,  for  example,  if  the  island  is  super-tidal,  then  primary 
breakers  will  be  oceanward  of  the  barrier.  The  application  of  the  li¬ 
near  wave  theory  stops  at  the  line  of  primary  breakers.  So  for  refrac¬ 
tion  purposes,  the  barrier  acts  as  the  shoreline.  Figure  5  shows  the 
resulting  digitization  of  the  shoreline.  Notice  the  inlet  to  the  right 
of  the  barrier  was  truncated. 

The  orientation  of  the  shoreline  with  respect  to  a  breaking  wave 
is  critical  in  calculations  dependent  on  breaker  height  and  breaker  an- 
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gle.  Two  such  calculations  used  in  the  program  are  the  longshore  energy 
flux  factor  and  longshore  sediment  transport  rate.  To  maintain  the  pro¬ 
per  shoreline  orientation,  a  constant  slope  is  assumed  from  the  first 
oceanward  depth  contour  onto  the  beach.  This  is  accomplished  by  digi¬ 
tizing  a  mirror  image  of  the  first  contour  onto  the  landward  side  of  the 
shoreline  (see  Figure  5,  the  -6  foot  contour).  A  second  mirror  image  of 
the  shoreline  may  be  digitized  landward  of  the  first  mirror  image  if 
changes  in  tidal  elevations  are  to  be  considered.  It  is  suggested  these 
lines  be  drawn  in  before  digitizing. 

STEP  3  -  Digitization  preparation.  Instructions  on  setting  up  the 
digitizer  have  been  previously  written  by  Jim  Crabtree  (personal  com¬ 
munication)  and  only  slightly  modified  in  the  following  explanation. 

The  digitizer  is  turned  on  by  a  switch  located  near  the  upper  cen¬ 
ter  of  the  back  panel.  Before  using  the  digitizer,  the  following 
switches  on  the  Operator's  Display  Panel  should  be  checked: 

1.  The  IRC  switch  should  be  in  the  on  position. 

2.  The  switch  labeled  DIST  POINT  TIME  MODE  should  be  in  the  DIST 
mode. 

3.  The  switch  labeled  AX  -  AX  +  AY  should  be  in  the  AX  +  AY  position. 

4.  The  magnetic  tape  unit  on-off  switch  should  be  turned  off. 

5.  The  DISTANCE  thumbwheels  should  be  set  to  .050  -  that  is,  a 
coordinate  will  be  recorded  every  5  hundredths  of  an  inch. 

To  mount  a  magnetic  tape  reel  on  the  tape  unit,  the  following  pro¬ 
cedure  should  be  followed: 

1.  The  magnetic  tape  unit  must  be  turned  off  during  the  loading 
of  a  tape  reel.  If  the  unit  is  on  during  the  loading  process, 
it  is  possible  to  blow  a  fuse. 

2.  Open  the  front  plastic  cover  to  the  tape  unit. 


3.  Mount  the  tape  reel  on  the  right-hand  hub  by  depressing  the  button 
in  the  center  of  the  hub  and  sliding  the  tape  reel  completely  onto 


the  hub  before  releasing. 

4.  Then,  using  one  hand  to  rotate  the  tape  reel,  follow  the  arrows  in¬ 
dicated  in  Figure  4  to  thread  the  tape  to  the  receiving  reel.  Make 
sure  that  the  tape  passes  between  the  metal  plate  and  the  write 
head  at  point  A  in  Figure  4. 

5.  Wrap  the  tape  around  the  receiving  reel  twice  and  make  sure  that 
all  slack  is  taken  up.  The  easiest  way  to  wrap  the  tape  around  the 
receiving  reel  is  to  rotate  the  receiving  reel  with  the  left  hand  un¬ 
til  the  index  finger  can  be  inserted  through  the  finger  hole  and 
press  the  tape  against  the  center  of  the  reel.  Then  rotate  both 
reels  clockwise  with  the  right  hand  giving  slack  and  the  left  hand 
taking  it  up. 

6.  Close  the  front  plastic  cover  and  turn  the  tape  unit  on. 

7.  Press  the  LOAD  button  once  and  release. 

8.  Press  the  LOAD  button  again  and  hold  until  the  READY  light  comes  on. 

9.  Press  the  FILE  GAP  button  and  release. 

10.  Press  the  FILE  GAP  button  again  and  release. 

11.  Press  the  REWIND  button  and  release. 

12.  Press  the  LOAD  button  and  hold  until  the  READY  light  comes  on. 

The  tape  unit  is  now  ready  to  record  data. 

Step  4  -  Digitizing  depth  data.  Tape  the  chart  to  the  top  of  the 
digitizer  table  so  that  the  x-axis  increases  to  the  right.  Located  on 
top  of  the  table  will  be  a  cursor  and  keyboard.  The  crosshairs  on  the 
window  of  the  cursor  are  used  for  guiding  the  cursor  as  the  data  is 
traced.  As  the  cursor  is  moved,  X  and  Y  coordinates  are  recorded  onto 
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Figure  4 


Magnetic  tape  unit. 
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Receiving  reel 

Sending  reel 
— Finger  hole 


Button 


magnetic  tape  every  .05  inch.  The  depth  corresponding  to  the  contour 
being  digitized  is  entered  on  the  keyboard. 

The  cursor  is  initially  set  at  the  lower  left  corner  of  the  perim¬ 
eter.  Pressing  the  red  button  on  the  cursor  resets  the  digitizer's  X 
and  Y  coordinates  to  zero.  Hold  the  cursor  steady  and  press  the  blue 
(DATA)  button.  This  places  the  digitizer  into  the  data-recording  mode. 
When  the  digitizer  is  in  this  mode,  the  red-orange  DATA  light  on  the 
Operator’s  Display  Panel  is  on,  and  data  can  be  recorded  from  the  cur¬ 
sor  or  the  keyboard.  Next,  press  the  yellow  (HOLD)  button  on  the  cur¬ 
sor.  The  digitizer  is  now  in  the  HOLD  mode  and  locks  in  on  the  last 
digitized  coordinate.  Only  data  from  the  keyboard  can  be  entered 
when  the  HOLD  light  on  the  Operator's  Display  Panel  is  on.  The  HOLD 
mode  prevents  the  recording  of  extraneous  data  if  the  cursor  is  ac¬ 
cidentally  moved.  The  digitizer  is  now  ready  to  record  the  first 
depth  entry. 

Depths  are  entered  as  a  5-digit  number  which  the  program  converts 
to  a  number  having  one  place  to  the  right  of  the  decimal.  A  depth 
entered  as  Df00525/s  will  be  converted  to  52.5  units  below  water 
level.  Depths  below  water  level  are  positive  (+)  and  elevations 
above  water  level  are  negative  (-).  The  few  remaining  codes  will  be 
discussed  as  they  are  used  in  the  digitizing  scheme  presented  below. 

The  bracketed  numbers  correspond  to  the  numbers  in  Figures  5,  6  and 
7. 

(1)  This  is  the  starting  point  on  the  chart.  With  the  crosshairs 
of  the  cursor  centered  on  the  lower  left  corner,  press  the  red  but¬ 
ton,  then  the  blue  and  yellow  buttons  on  the  cursor.  This  point  is 
recorded  as  the  zero  reference  point.  The  digitizer  is  in  the  HOLD 
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Figure  5.  Digitization  of  perimeter  on 
hydrographic  chart. 


hypothetical 
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mode  which  records  keyboard  data  entry  only.  Corner  depths  are  dif¬ 
ferentiated  from  depths  on  the  perimeter  by  the  notation  'PD'.  All 
other  depths  begin  with  'D'.  Using  Figure  5  as  an  example,  enter 
PD+00400/s  on  the  keyboard.  In  other  words,  the  depth  at  this  point 
is  40  feet.  Press  the  yellow  button  to  take  the  digitizer  out  of  the 
HOLD  mode.  Cursor  movement  is  now  able  to  be  recorded.  With  a  steady 
motion,  move  the  cursor  to  (2). 

(2)  Press  the  HOLD  button  to  put  the  digitizer  in  the  HOLD  mode. 
Enter  Df 00300/s  on  the  keyboard.  Press  the  HOLD  button  to  take  the 
digitizer  out  of  the  HOLD  mode  and  move  the  cursor  to  (3). 

(3)  Put  the  digitizer  in  the  HOLD  mode  and  enter  PD4-00250/s  on  the 
keyboard.  Press  the  HOLD  button  to  take  the  digitizer  out  of  the  HOLD 
mode  and  move  the  cursor  to  the  first  depth  that  crosses  this  section 
of  the  perimeter  (4) . 

(4)  Press  the  HOLD  button.  Enter  DH)0180/s  on  the  keyboard.  Press 
the  HOLD  button  and  move  the  cursor  to  the  next  point. 

Continue  digitizing  the  perimeter  in  the  same  fashion  until  point 

(21). 

(21)  The  last  point  on  the  perimeter  is  the  starting  point.  Press 
the  HOLD  button  to  enter  the  HOLD  mode  and  enter  'LS'  on  the  keyboard. 
This  signifies  the  end  of  the  perimeter.  The  contours  are  digitized 
next.  Contour  depths  are  denoted  by  a  'D'  except  for  the  shoreline 
which  is  symbolized  by  an  'H'.  Refer  to  Figure  6  for  the  following 
discussion  on  digitizing  the  contours. 

(21)  Press  the  HOLD  button  to  take  the  digitizer  out  of  the  HOLD 
mode.  Steadily  move  the  cursor  to  the  first  contour  to  be  digitized, 
in  this  example  (22). 


Figure  6.  Digitization  of  contours  and  shoreline  on 
hypothetical  hydrographic  chart. 

The  shoals  outlined  by  a  dashed  rectangle  are  enlarged 
in  Figure  7. 
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(22)  Press  the  HOLD  button  to  put  the  digitizer  in  the  HOLD  mode. 
Enter  IH-00300/  on  the  keyboard.  A  contour  is  signified  when  an  'S' 
does  not  immediately  follow  '/'•  Take  the  digitizer  out  of  the  HOLD 
mode  and  carefully  trace  the  contour.  Remember,  contour  traces  will 
be  plotted  on  the  final  output. 

(23)  Enter  the  HOLD  mode  and  enter  an  S  on  the  keyboard.  This 
completes  the  trace  of  the  30  foot  contour.  Press  the  HOLD  button  to 
take  the  digitizer  out  of  the  HOLD  mode  and  move  the  cursor  to  the  next 
contour. 

(24)  Press  HOLD.  Enter  D+-00180/  on  the  keyboard.  Trace  the  18 
foot  contour. 

(25)  Press  HOLD  and  enter  the  letter  S  on  the  keyboard.  Press 
HOLD  and  move  the  cursor  to  the  next  contour. 

The  remaining  contours  are  digitized  in  the  same  manner,  except 
the  shoreline.  Replace  D  and  an  H  for  the  shoreline  and  enter  H+00000/ 
on  the  keyboard.  As  in  the  contour  depths,  the  end  of  a  segment  of 
shoreline  is  denoted  by  an  S.  Refer  to  Figure  7  for  the  following 
bracketed  numbers. 

(37)  After  the  S  is  entered  on  the  keyboard  signifying  the  end  of 
the  -6  foot  trace,  press  the  HOLD  to  take  the  digitizer  out  of  the  HOLD 
mode.  Move  the  cursor  to  a  point  of  one  of  the  shoals  (38). 

(38)  Press  HOLD  and  enter  the  dpeth  of  the  shoal,  D+00010/  ,  on  the 
keyboard.  Press  HOLD  and  with  a  steady  motion,  trace  the  shoal. 

(39)  Upon  returning  to  the  starting  point,  press  HOLD  and  enter  an 
S  on  the  keyboard.  Press  HOLD  again  and  move  the  cursor  to  the  next 
shoal. 

(40)  Press  HOLD.  Enter  D+00010/  on  the  keyboard.  Press  HOLD,  tak¬ 
ing  the  digitizer  out  of  the  HOLD  mode  and  trace  the  contour. 
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Figure  7.  Enlarged  view  of  shoals  digitized  in 
Figure  6. 


(41)  Again  upon  returning  to  the  starting  point,  press  HOLD  and 
enter  an  S  on  the  keyboard.  This  being  the  last  contour  to  trace,  en¬ 
ter  an  F  on  the  keyboard  (after  the  S)  to  signify  the  end  of  the  digi¬ 
tized  data. 

(To  Finish)  The  trace  being  complete,  press  HOLD  to  take  the  digi¬ 
tizer  out  of  the  HOLD  mode  and  move  the  cursor  in  any  direction  until 
a  beep  is  heard  on  the  digitizer.  This  will  create  a  complete  record 
block  of  4096  characters  on  the  magnetic  tape  which  is  necessary  for 
proper  operation  of  the  program.  On  the  Operator  Display  Panel,  press 
the  FILE  GAP  button  twice.  This  will  indicate  the  end  of  this  data 
file. 

Step  5  -  Dismounting  the  magnetic  tape.  The  following  steps  are 
to  be  followed  when  unloading  the  magnetic  tape  unit: 

1.  Press  the  REWIND  button  and  release. 

2.  After  the  magnetic  tape  has  finished  rewinding,  turn  the  mag¬ 
netic  tape  unit  off. 

3.  Open  the  front  plastic  cover  to  the  tape  unit. 

4.  Use  the  left  hand  to  rotate  the  left-hand  reel  counterclock¬ 
wise  to  give  tape  slack.  At  the  same  time,  use  the  right  hand 
to  rotate  the  right-hand  reel  counterclockwise  to  take  up  the 
slack  tape.  Continue  this  process  until  the  tape  has  been 
wound  onto  the  right-hand  reel. 

5.  Remove  the  tape  reel  from  the  right-hand  hub  by  depressing  the 
button  in  the  center  of  the  hub  and  sliding  the  reel  off  before 
releasing. 

6.  If  no  more  tapes  are  to  be  mounted,  make  sure  the  front  plastic 
cover  is  closed.  This  will  help  to  keep  dust  off  of  the  write- 
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head  assembly. 


Step  6  -  Final  Step.  Only  one  bathymetric  map  can  be  digitized  per 
magnetic  tape.  To  avoid  the  problem  of  running  out  of  tapes,  a  small 
program  written  by  J.  Crabtree  will  copy  the  digitized  data  off  the 
magnetic  tape  and  store  it  onto  disk,  thereby  freeing  the  tape  for 
later  use. 

The  tape  first  has  to  be  brought  to  the  user's  service  window  in 
Computer  Services  Division  (CSD) .  At  the  time  of  this  writing  CSD 
is  located  at  the  corner  of  Wheat  and  Main  Streets.  Remember  to  copy 
the  number  of  the  tape  as  this  information  is  needed  by  the  program. 

A  program  listing  and  information  on  how  to  implement  this  program  is 
provided  in  PROGRAMS:  TTD  (Tape  to  Disk). 
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PROGRAMMING  CONSIDERATIONS 
Determination  of  Input  Parameters  • 

MI,  MJ,  SCX  and  GRID.  Refer  to  the  Digitized  Depth  Input  (Step  1) 
for  the  calculation  of  MI,  MJ,  and  SCX.  GRID  is  the  dimension  of  the 
side  of  a  grid  square  expressed  in  map  feet  and  is  a  function  of  SCX 
and  the  scale  of  the  map.  GRID  is  given  by 
GRID  (feet)  =  (SCX  x  map  scale)  /12 
For  example  if  Figure  3  had  a  1:80,000  scale,  GRID  would  equal  666.67 
feet  (.1  x  80000/12),  if  SCX  -  .1. 

XSG ,  YSG.  SCNV.  DGXL  and  DGYL  -  Window  parameters.  A  window  is  an 
area  inside  the  present  grid  being  processed.  This  window  enables  the 
user  to  study  a  particular  area  in  the  grid  in  greater  detail  by  using 
a  larger  scale  grid.  REFRAC  will  trace  the  rays  through  the  window  and 
generate  the  necessary  input  parameters  to  be  used  in  a  refraction  run 

on  the  larger  scale  grid.  To  create  the  window  on  the  small  scale  grid, 

V 

the  X  (XSG)  and  Y  (YSG)  coordinates  of  the  lower  left  corner  of  the  win¬ 
dow  position,  the  ratio  of  the  two  scales  (SCNV)  and  the  length  of  the 
axes  of  the  window  (DGXL  and  DGYL)  expressed  in  inches  of  the  small 
scale  grid  are  input  in  the  preliminary  run.  XSG  and  YSG  are  given  by 


distance  (inches)  from  respective  X  or  Y  axis  of  the 
small  scale  grid  to  the  lower  left  corner  of  the  window 

SCX 


SCNV  may  be  calculated  by 

SCNV  =  small  map  scale/larger  map  scale 
SCNV  is  considered  a  magnification  factor  and  will  always  be  greater 
than  1. 
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The  following  equation  yields  DGXL  and  DGYL: 

DGXL  or  DGYL  =  respective  axes  length  (inches)  /  SCNV 
(inches)  of  larger  scale  map 

This  initial  run  on  the  small  scale  grid  will  generate  the  starting 
coordinates  (X  and  Y) ,  the  refraction  coefficients  of  that  time  step 
(RK)  and  of  the  previous  time  step  (RK1)  and  the  azimuth  (AZIMTH)  for  a 
ry  to  be  used  as  input  in  a  subsequent  run  on  the  larger  scale  grid. 
These  wave  characteristics,  generated  for  every  ray  that  crosses  the 
window  boundary,  can  be  either  printed  onto  paper,  punched  onto  cards 
or  stored  on  disk.  The  program  is  presently  set  up  to  store  these 
values  on  disk  under  a  user  determined  data  set  name,  so  they  can  be 
easily  accessed.  See  JCL  considerations  under  REFRAC  in  the  section  on 
PROGRAMS  for  further  information.  If  RK,  RK1,  and  AZIMTH  are  not  the 
results  of  a  previous  run,  they  will  have  to  be  hand  calculated  for 
each  ray  and  keypunched  onto  cards  for  input  into  the  program.  (This 
only  applies  to  rays  starting  in  shallow  water.) 

NORAYS.  A  maximum  of  400  rays  can  be  refracted  per  set.  If  the 
number  of  rays  is  not  known,  such  as  in  the  case  of  window-generated 
input,  set  NORAYS  equal  to  400. 

SK  and  SKI.  SK  is  the  shoaling  coefficient  at  water  depth  for 
the  present  time  step.  If  a  set  of  rays  start  in  shallow  water  (SK 
is  greater  than  1),  set  SK  to  1  on  input  because  the  program  will 
calculate  the  proper  value.  SKI  is  the  shoaling  coefficient  at  water 
depth  for  starting  location  of  a  ray.  Most  refraction  diagrams  are 
of  swell  conditions;  therefore,  SKI  usually  equals  1. 

THI  and  STAZ.  THI  is  the  clockwise  angle  between  north  on  the 
map  and  the  positive  y-axis  of  the  grid.  See  Figure  8.  STAZ  is  the 
azimuth  or  the  clockwise  angle  between  north  on  the  map  and  the  deep- 
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Figure  8.  THI  is  the  clockwise  angle  between 
north  on  the  map  and  the  positive  y-axis  of  the 
grid. 


water  direction  of  wave  approach.  See  Figure  9. 


GRINC.  UNIT.  WPI.  GRID  and  KPLOT.  UNIT  is  the  time  step  for  a  set 
of  rays.  WPI  is  the  number  of  wave  periods  between  tick  marks  on  a  ray. 
The  variable  KPLOT  is  calculated  in  the  program  and  determines  the  points 
on  a  ray  where  tick  marks  will  be  plotted.  KPLOT  is  dependent  on  the 
values  of  UNIT,  WPI  and  GRID.  GRINC  is  the  deepwater  step  length,  ex¬ 
pressed  as  a  fraction  of  a  grid  square  and  is  given  by: 

GRINC  =  UNIT  x  CO/GRID. 

CO,  the  deepwater  wave  celerity,  equals  5.12  times  the  period  (T) .  In 
the  above  equation,  UNIT  is  the  only  variable  user  determined  and 
should  be  chosen  so  the  resulting  value  of  GRINC  is  approximately  .5. 

KPLOT  is  given  by: 

KPLOT  =  (T  x  WPI/UNIT)  +  .1 

WPI  should  be  chosen  so  the  value  of  KPLOT  is  greater  than  or  equal  to 
1.  Optimal  results  are  achieved  if  KPLOT  is  approximately  10. 
Miscellaneous 

Depending  on  the  needs  of  the  user,  there  is  the  choice  of  printing 
or  not  printing  the  depth  grid  values  or  wave  information.  Refer  to 
JCL  Considerations  if  either  of  these  options  is  desired.  The  amount 
of  printed  output  may  be  decreased  by  both  of  the  variables  NPRINT  and 
PD.  NPRINT  is  the  user  determined  printing  interval  for  information 
along  a  ray.  A  suggested  value  for  NPRINT  is  10.  PD  is  the  depth  at 
which  printing  will  begin.  The  advantage  of  PD  is  that  it  enables  the 
user  to  limit  printing  to  where  the  rays  are  more  strongly  refracted  - 
that  is,  in  shallow  depths.  Over  10,000  lines  of  printout  can  easily 
be  generated  when  refracting  a  set  of  more  than  50  rays.  Turnaround 
time  (time  it  takes  to  get  printed  results)  and  CPU  time  (execution 
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Figure  9.  Determination  of  deepwater  azimuth  for  a  ray  ex 
pressed  in  degrees  (STAZ) .  Wave  approach  depicted  by  di 
rection  of  arrow. 


time)  are  decreased  when  the  number  of  lines  printed  is  kept  Co  a  mini¬ 


mum. 

Keypunching  and  computation  time  can  be  saved  by  the  two  input  vari¬ 
ables  IWR  and  YFA.  The  card  input  falls  into  two  categories,  the  infor¬ 
mation  needed  for  a  set  of  rays  and  the  information  on  individual  rays 
such  as  starting  coordinates.  If  more  than  one  set  of  rays  is  being  pro¬ 
cessed  and  the  individual  wave  characteristics  will  be  the  same,  set  IWR 
equal  to  1  for  the  sets  after  the  first  set.  The  program  will  store  the 
first  set  of  ray  characteristics  and  use  the  same  information  for  the  sub¬ 
sequent  sets  until  all  the  sets  of  rays  are  processed  or  until  IWR 
equals  zero.  If  IWR  equals  zero,  the  individual  wave  characteristics 
have  to  be  read  off  cards  or  disk.  YFA  adds  a  constant  to  the  y-co- 
ordinates  of  the  starting  locations  for  a  set  of  rays.  This  enables 
the  same  wave  characteristics  cards  to  be  used  for  a  variety  of  periods 
by  moving  the  starting  positions  of  rays  relative  to  the  shoreline. 

The  user  has  a  choice  between  two  plotting  devices,  the  Gould  elec¬ 
trostatic  plotter  or  the  Calcomp  drum  plotter.  The  Gould  plotter  is 
faster  -  minutes  compared  to  hours  -  but  the  Calcomp  plotter  has  better 
resolution.  The  increased  resolution  may  be  desired  if  the  refraction 
diagrams  are  to  be  photographed. 

The  Gould  plotter  has  a  maximum  length  (x-axis)  of  327  inches  and 
height  (y-axis)  of  63  inches.  The  width  of  the  plotter  is  only  10.5 
inches;  therefore,  any  width  greater  than  10.5  inches  (but  less  than 
63)  will  cause  the  plot  to  be  'stripped'.  The  input  variables  FAC  or 
BOUND  can  reduce  the  plot  size  to  avoid  stripping.  Because  of  the  fast 
plotting  time,  the  Gould  plotter  is  excellent  for  debugging  or  obtain¬ 
ing  work  copies  of  the  refraction  diagram. 

The  maximum  length  of  the  x-axis  for  the  Calcomp  plotter  is  limited 
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to  the  length  of  paper  on  the  roll.  The  maximum  length  of  the  y-axis 
is  29.5  inches.  For  purposes  of  photocopying  the  diagram,  liquid  ink 
can  be  requested  instead  of  the  standard  black  ballpoint  pen. 
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REFRAC 


JCL  Considerations 

The  JCL  (Job  Control  Language)  presented  in  this  section  is  coded  for 
an  IBM  370/168  System.  Job  control  cards  begin  with  a  //  (slash-slash). 
REFRAC  is  compiled,  loaded  and  then  stored  on  disk  (mass  storage)  in  a 
program  library  called  WORK. 

The  following  discussion  explains  how  to  implement  REFRAC  and  the 
various  options  available  by  small  changes  in  the  JCL.  (Refer  to  JCL 
listing.)  The  first  3  cards  are  the  JOB  card.  Included  on  the  JOB 
card  are  the  account  number  (N3100138) ,  maximum  number  of  lines  to  be 
printed  (50,000)  and  maximum  number  of  plot  records  (99990).  TIME  has 
the  format  (minutes,  seconds);  therefore,  the  time  requested  in  the 
above  JCL  listing  is  2  minutes  and  30  seconds.  The  time  can  be  de¬ 
creased  depending  on  the  number  and  size  of  diagrams  to  be  plotted. 

USER  and  PASSWORD  information  is  obtained  when  an  account  number  and 
disk  space  are  requested  from  Computer  Services  Division. 

The  remainder  of  the  JCL  between  the  fourth  card,  PR0C,  and  the 
last  card,  PEND,  is  called  an  in-stream  procedure.  This  procedure  ac¬ 
cesses  the  necessary  libraries  and  data  sets  whenever  it  is  called. 

A  call  to  the  procedure  is  made  by  an  EXEC  statement  following  the  PEND 
statement.  More  detail  is  given  in  the  several  examples  following  the 
listing  of  the  in-stream  procedure. 

The  REFRAC  procedure  is  set  up  to  plot  the  results  on  the  Gould 
plotter  at  the  Social  and  Behavior  Sciences  Lab  located  in  Gambrell 
Hall  (University  of  South  Carolina).  The  Calcomp  drum  plotter  is  lo¬ 
cated  in  Computer  Services  Division  (CSD)  and  the  results  can  be 
p Lotted  there  with  the  following  minor  changes  to  the  JCL.  Change: 
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JCL  LISTING  OF  THE  IN-STREAM  PROCEDURE  REFRAC 


//N3100138  JOB  (N3100138.50.9999) ,'WAVE  REFRACT  ION • .MSGLEVEL= (1 . 1 ) ♦ 
//  USER=N310013.PaSSWOPD=CRD.REGION=600K. 

//  T IME= (2.30) 

//REFRAC  PROC  DEHGR ID= ‘DUMMY « • , WAVDA T A=NULLF ILE  »NEWW AV=NULLF 1LE 

//STP1  EXEC  PGM=REFRAC 

//STEPLIB  DO  DSN=N31 00 1 3. WORK ,OISP=SHR 

//  OD  DSN=SM1.LINKLIB,DISP=SHR 

//SYSPRINT  DO  SYSOUT  =  A 

//SYSOUT  DO  SYSOUT  =A 

//S0RTLI8  DD  DSN=SM1 .SORTL Ih.DISP=SHR 

//SORT  WKO 1  Of)  UNIT=SYSDA.SPACE=(TRK.200..CONTIG) 

//SORTWK02  DD  UNIT=SYSDA. SPACE* (TRK.200. .CONT IG) 

//SORTWK03  DD  UN IT=SYSDA .SPACE= ( TRk ♦ 200. .CUNT IG) 

//FT 06F001  DD  SYSOUT  =  A 

//FT07F00 1  DD  DSN=LNEWWAV,UNIT=3330V,MSVGP=USCP» 

//  DISP=(NEW«CATLG. DELETE) .SPACE=(TRK. (5.5) .RLSE ) ♦ 

//  DCe=(LPECL=50.BLKSIZF,=2000.RECFM=FB) 

//FT08F001  DD  &DEPGR ID • SYSOUT  =  A 

//FT  1 1F001  DD  DSN=iV31 001 38.004  ,DISP  = (NEW.DELtTE. DELETE) . 

//  UNIT=SYSDA.SPACE=(CYL. 10) ♦ 

//  DC8=(8LKSIZE=12000.LRECL=16.BUFNO=1.RECFM=FB> 

//FT12F001  DD  DSN=N3 1 003 1 8 . DO  1 . DI SP= ( NEW . DELETE .DELETE ) ♦ 

//  UNIT  =SYSDA» SPACE = (CYL , 1 0 ) ♦ 

//  OCB=(8LKSIZE=12000.LRECL=12.8UFNO=1,RECFM=FR) 

/ /SOH  TIN  DD  OSN=*. FT  1 2F00 1 , D ISP* < OLD, DELETE, DELETE) . 

//  UNIT=SYSDA,VOL=REF=*.FT12F001, 

//  DCB* (BLKS I ZE= 12000 «LRECL=12»BUFNO*l ,RECFM=FB) 

//SORTOUT  DO  DSN=N3100138.D02»DISP= (NEW, DELETE .DELETE) ♦ 

//  UNIT  =SYS0A*SPACE=(CYL,10) , 

//  DCB=(BLKSIZE=12000,LRECL=12,BUFNO=1,RECFM=FB) 

//FT13F001  DD  OSN=*.SORTOUT,DISP= (OLD, DELETE, DELETE) . 

//  UNIT=SYSUA,VOL=REF=«. SORTOUT, 

//  DCB=(BLKSIZE=120U0,LRECL=12,BUFNO=l.RtCFM=FB) 

//FT  1  OF  00 1  DD  OSN=N310013B. D03.DI5P* (NEW, DELETE, DELETE) , 

//  UNIT=SYSDA,SPACE*(CYL,2> , 

//  DCB=(SLKSIZE=12000»LRECL=12,BUFNO=1»RECFM=FB> 

//FT17F001  OD  DSN=&DEPOATA,DISP=SHR,LABEL=(»»»IN> 

//FT25F00 1  DD  OSNs&WAVDAT A ,DI5P=SHR ,L ABEL* ( , * , IN) 

//FT05F001  DD  DDNAME=SYS IN 

//SYSVECTR  DD  OSN=&&VECTORS»UNIT=SYSDA,DISP=(NEW,PASS) ,  ■ 

//  SPACE* ( 1320, (500,200) ) ,nCB=BLKSIZE=1320 

//SYSPOUT  DD  SYSOUT  =  A 

//STP2  EXEC  PGM=BTSTN50 

//STP2.SYSUDUMP  DO  SYSOUT=A 

//STEPLIB  DO  DSN=ACAD.SUBLIB.G5005»DISP=SHR 

//SYSVECTR  DD  DSN=^4VFCT0R5,DISP= (OLD .DELETE) 
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//SYSUT1  OD  OSN=&*SYSUTl»UNIT=SYSDA*SPACEs<ltt?00t200> 
//SYSPLOT  DO  SYSOUT=(b,.PLOT) ,DESTsRMT8 
//SYSPOUT  DO  SYSOUTsA 
//  PEND 


CALLS  TO  ABOVE  PPOC 


1.  card  #5  (EXEC  card)  to 


//ST PI  EXEC  PGM=CALCOMP 

2.  replace  the  10  cards  between  (excluding) 

//FT05001  DD  DDNAME=SYSIN 
and 

//  PEND 

with  one  card 

//  GO.PLOTTAPE  DD  DSN=DRUM. PLOT,UNIT=TAPED, 

DISP=  (NEW, KEEP)  ,DCB=DF.N=2 

If  the  Calcomp  plotter  is  used,  a  DP  number  is  returned  on  the  front 
page  of  the  printout.  This  number  is  the  location  of  the  plot  on 
tape  and  needs  to  be  called  into  CSD  in  order  to  obtain  the  plotted  re¬ 
sults  . 

Three  types  of  data  sets  are  of  primary  concern  to  the  user.  The 
first  type  is  the  digitized  depth  data,  DEPDATA,  which  resides  on  disk 
having  been  previously  copied  from  magnetic  tape  by  the  program  TTD. 
(User  must  have  previously  requested  disk  space  from  CSD.)  The  second 
type,  NEWWAV ,  is  created  during  a  run  in  which  window  boundaries  had 
been  specified.  The  output  wave  data,  NEWWAV  (X,  Y,  AZIMTH,  RK,  RK1) , 
is  stored  on  disk  to  be  read  in  a  subsequent  run.  As  input  data  NEW¬ 
WAV  is  called  WAVDATA,  which  is  the  third  type  of  data  set. 

Data  set  names  (DSN)  may  consist  of  up  to  44  alphanumeric  charac¬ 
ters  (including  periods  used  for  separation).  A  period  must  separate 

f 

each  8  characters  -ot*  less  and  must  be  followed  by  ar.  alphabetic  charac- 

i 

ter.  The  first  7  characters  of  the  DSN  are  the  first  7  characters  of 
the  8  character  account  number.  Data  set  names  must  be  unique  or  an 
error  will  occur  upon  the  creating  the  data  set.  A  data  set  name  such 
as  N3100L3.  I’KICK. INLET  could  be  used  for  the  account  number  N3100138  and 
represent  a  data  set  (of  any  of  tile  three  types)  for  Price  Inlet,  South 
Carol i nn , 
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The  general  JCL  format  for  the  program  is: 

//Job  Card(s) 

//REFRAC  PROC  •  •  •  • 

.  in-stream  procedure 

//PEND 

/* 

//stepname  EXEC  REFRAC,  data  set  type='data  set  name', 
data  set  type=*data  set  name* 

//SYSIN  DD 

.  card  input  as  described  in  PROGRAM  DESCRIPTION 

//stepname  EXEC  REFRAC,  data  set  type='data  set  name' 

//SYSIN  DD 

.  card  input 

// 

As  mentioned  earlier,  a  call  to  the  in-stream  procedure  (REFRAC) 
sets  up  the  program  libraries  and  data  sets  necessary  for  the  refrac¬ 
tion  program.  A  call  to  the  procedure  is  made  when  an  EXEC  REFRAC  card 
is  encountered  in  the  JCL.  a  job  may  consist  of  more  than  one  call  by 
simply  adding  another  EXEC  REFRAC  card  and  its  corresponding  data  cards 
to  the  last  data  input  card  from  the  previous  call.  The  above  example 
is  termed  a  multistep  job  because  it  calls  the  in-stream  procedure  twice. 

The  program  is  terminated  by  a  //  in  columns  1  and  2  of  the  last  card. 

The  stepname  is  from  1  to  8  alphanumeric  characters,  the  first  of 
which  must  be  alphabetic.  The  purpose  of  the  EXEC  REFRAC  card  is  to 
specify  the  data  sets  that  are  to  be  read  or  created  during  the  execu¬ 
tion  of  the  refraction  program.  Examples  of  data  set  types  are  explained 

below.  The  various  options  are  also  shown. 

1.  A  simple  refraction  program  reads  the  depth  data  off  disk, 
doesn't  specify  a  window  and  reads  all  wave  parameters  off  cards.  In 
other  words,  of  the  3  data  types,  only  DEPDATA  needs  to  be  specified. 

//STP  EXEC  REFRAC, DEPDATA='N 310013. SOUTH. CAROLINA* 

2.  The  above  example  results  in  a  refraction  diagram  of  the  South 
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Carolina  coast.  In  the  present  example,  a  window  is  specified  in  the 
card  input;  therefore,  a  data  set  name  will  have  to  be  given  for  NEWWAV 
to  store  the  wave  information  generated.  If  the  window  is  of  Price  In¬ 
let,  then 

/ / STP  EXEC  REFRAC,DEPDATA='N310013. SOUTH. CAROLINA' , 

/ /  NEWWAV= ' N 310013 . PRICE . INLET . WAVE ' 

This  data  set  may  now  be  read  in  a  subsequent  refraction  analysis  on  a 
larger  scale  chart  of  Price  Inlet. 

3.  Using  the  above-created  data  set  and  the  digitized  depth  data 
for  Price  Inlet,  a  more  detailed  analysis  can  be  obtained. 

//STP  EXEC  REFRAC, DEPDATA='N310013. PRICE. INLET' , 

/ /  WAVDATA= ' N  310013 . PRICE . INLET . WAVE ' 

Note:  In  this  example,  the  card  input  variables  IWR  and  IRED  EQUAL 
0  and  25  respectively  (see  PROGRAM  DESCRIPTION  (Card  Input)). 

4.  The  printing  of  the  depth  grid  is  suppressed  unless  explicitly 
requested.  This  is  accomplished  by  punching  DEPGRID=  in  the  data  set 
type  location. 

//STP  EXEC  REFRAC, DEPDATA='N310013. DEPTH. DATA' , 

//  DEPGRID  = 

If  NEWWAV  is  also  specified  then 

//  DEPGRID= ,NEWWAV= 'N310013. WAVE ' 

5.  The  printing  of  the  title  page  and  all  wave  information  may  be 
suppressed  if  only  the  plot  is  desired.  Refer  to  card  input  listing 
for  proper  positioning  of  cards  described  in  examples  5  and  6.  Add 
the  following  JCL  card  after  the  'EXEC  REFRAC'  card  and  before  the 
'SYS IN  l)U  *'  card 

//STP1 . FT06F001  DD  DUMMY 

6.  The  wave  parameters  generated  for  a  window  may  be  punched  on 


cards  instead  of  stored  on  disk  by  adding  the  following  JCL  card 
between  the  'EXEC  REFRAC'  card(s)  and  'SYSIN  DD  *'  card 
/ /STP1 . FT07F001  DD  SYSOUT-B 
Program  Listing,  Card  Input  and  Example  of  Output 

The  following  listing  is  a  deck  setup  including  both  JCL  and  card 
input.  A  refraction  analysis  of  South  Carolina  coast  results  from  the 
first  call  to  the  in-stream  procedure.  The  wave  information  generated 
at  the  window  boundaries  is  used  in  the  second  call  of  the  procedure, 
thereby  producing  a  detailed  analysis  of  the  window  area.  Refer  speci¬ 
fically  to  examples  2,  3,  and  5  for  further  explanations. 

A  sample  of  the  output  includes  the  title  page,  information  on  one 
ray  and  the  refraction  diagram.  The  information  is  from  the  window  area 
analyzed  above.  Following  the  sample  output  is  a  complete  listing  of 
the  program  REFRAC  and  the  necessary  JCL  to  compile,  load  and  store  it 
on  disk  (if  it  has  not  already  been  stored). 


PROGRAMS 


COASTAL  RESEARCH  CENTliR  TECHNICAL  REPORTS 
(Geology  Department,  Univ.  of  Massachusetts) 


1.  Offset  coastal  inlets. 

M.  0.  Hayes,  V.  Goldsmith,  and  C.  It.  Hobbs;  Tech.  Rept.  No.  1-CRC,  1972. 

2.  forms  of  sediment  accumulation  in  the  beach  zone. 

M.  0.  Hayes;  Tech.  Rept.  No.  2-CRC,  1972 

3.  Hydraulic  equivalent  sediment  analyzer  (IIESA) . 

F.S.  Aiian;  Tech.  Rept.  No.  3-CRC,  1972 

4.  Glacial  environment  and  processes  in  southeastern  Alaska. 

J.H.  Hartshorn  and  G.M.  Ashley;  Tech.  Rept.  No.  4-CRC,  1972. 

5.  Sedimentation  and  physical  limnology  in  proglacial  Malaspina  Lake,  Alaska. 

T.C.  Gustavson;  Tech.  Rept.  No.  S-CRC,  1972. 

COASTAL  RESEARCH  DIVISION  TECHNICAL  REPORTS 
(l)cpt.  of  Geology,  Univ.  of  South  Carolina) 

6.  Coarse-grained  sedimentation  on  a  braided  outwash  fan,  northeast  Gulf  of  Alaska. 
J.C.  Boothroyd;  Tech.  Rept.  No.  6-CRD,  1972. 

7.  Bank  margin  morphology  and  sedimentation,  Lucaya,  Grand  Bahama  Island. 

I). k .  Hubbard,  L.G.  Ward,  U.  M.  FitzGerald,  and  A.  C.  Hine;  Tech.  Rept.  No.  7-CRD, 
1976. 

S.  Selected  environmental  criteria  for  the  design  of  artificial  structures  on  the 
southeast  shore  of  Lake  Erie. 

R.ll.  Clemens;  Tech.  Rept.  No.  8-CRD,  1976. 

9.  Coastal  dynamics  and  sediment  transportation,  northeast  Gulf  of  Alaska. 

D.  Nummedal  and  M.F.  Stephen;  Technical  Rept.  No.  9-CRD,  1976. 

10.  Morphologic  and  hydrodynamic  characteristics  of  terrestrial  fan  environments. 

D.  Nummedal  and  J.C.  Boothroyd;  Tech.  Rept.  No.  10-CRD,  1976. 

11.  Terrigenous  clastic  depositional  environments. 

M.  0.  Hayes  and  T.  W.  Kana,  eds.,;  Tech.  Rept.  No.  11-CRU,  1976. 

12.  Coastal  morphology  and  sedimentation  -  Lower  Cook  Inlet,  Alaska. 

M.  0.  Hayes,  I’.  J .  Brown,  and  J.  Michel;  Tech.  Rept.  No.  12-CRD,  1976. 

13.  The  Carolina  Bays:  A  comparison  with  modern  oriented  lakes. 

R.  T.  kaezorowski;  Tech.  Rept.  No.  13-CRD,  1977. 

14.  Variations  in  tidal  inlet  processes  and  morphology  in  the  Georgia  Embayment. 

D. k .  Hubbard;  Tech.  Rept.  No.  14-CRD,  1977. 

15.  Coastal  morphology,  sedimentation  and  oil  spill  vulnerability  -  Northern  Gulf  of 
Alaska. 

C.  II.  Ruby,  Tech.  Rept.  No.  15-CRD,  1977. 

16.  Two  years  after  the  Hetula  oil  spill.  Strait  of  Magellan,  Chile: 
oil  interaction  with  coastal  environments. 

Anne  E.  Blount,  Tech.  Rept.  No.  16-CRD,  1978. 

17.  Influence  of  wave  refraction  on  coastal  geoinorphology  -  Bull  Island 
to  Isle  of  Palms,  South  Carolina. 

C.  Fico,  Tech.  Rept.  No.  17-CRD,  1978. 
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LISTING  OF  THE  PROGRAM  PEFRAC  ANO  THE  NF.CtSSARY  JCL  TO 
COMPILE*  LINK  AMO  STOWE  IT  ON  DISK 


r 


//SETUP  JOH 
//  T1ME=U1B> 
//  EXEC  FTGICf 
//SYSIN  OD  # 

C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c. 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
r. 
c 


JlOftUH.R)  «SF  TUP.MSGLEVEL= (  1  *1)  »PEGI0N=2bQK. 

.USER=N310013.PASSwORD=C»D 

PAWH.FOWTs  •LOAn,NOiJECK*MOSOUPCt  • 


AZIMTH  - 

HOUND  - 


CKQEP  - 

DATE!  - 
DCOM 

ofltas  - 

OF  P 
OF 


DGXl 

DGYL 

FAC 


GW  1 0 

GP INC  - 

HO 

ISP 


WAVE  refraction  diagrams 
PWOoWAM  TO  COMPUTE. PRINT. AND  PLOT  WAVE 

UI A6RAMS 


REFRACTION 


ITITLF  - 


IWR 


DFF INITION  OF  INPUT  PAWAMFTFPS 
AZIMUTH  OF  ,  AY. 

DIMENSION  D  DIRECTION  OF  PORTION  OF  MAP  THAT  WILL  NOT 
HE  SHOWN  ON  PL  EXPRESSED  In  INCHES  AND  ALWAYS  IS 
SUBTRACTED  OFF  octPwATFR  END. 

CHFCK  DEPTH.  PROGRAM  WILL  PRINT  RAY  INFORMATION  FOP  THE 
TIME  ST F P  NEAREST  THIS  DEPTH. 

UATE  OF  hYOROGHAPHIC  CHART 

CONVERSION  FACTOR  FOR  DEPTH  VALUES  TO  FEET. 

MINIMUM  STc.P  LENGTH.  EXPRESSED  AS  FRACTION  OF  6RID 
SQUARE. 

WATER  DFPThS  AT  EACH  GRID  POINT. 

FACTOR  TO  CONVERT  DEPTH  VALUES  FROM  ONE  WATER  SURFACE 
DATUM  TO  ANOTHER.  THE  FACTOR  WILL  BE  ADDED  TO  THE  DEPTH 
VALUE.  IF  NOT  NEEDED  LEAVE  HLANK. 

LENGTH  OF  WINDOW  (X  DIMENSION)  EXPRESSED  IN  INCHES. 
HEIGHT  OF  WINDOW  <Y  DIMENSION)  EXPRESSED  IN  INCHES. 
SCALING  FACTOR  HY  WHICH  ALL  POINTS  ARE  MULTIPLIED 
BEFORE  PEING  PLOTTED. 

DIMENSION  OF  A  SIDE  OK  A  GRID  SQUARE.  EXPRESSED  IN 
MAP  FEET. 

DFFPWATER  STEP  LENGTH.  EXPRESSED  AS  FRACTION  OF  GRID 
SQUARE. 

DEEPWATER  WAVE  HEIGHT 

SETS  PRINT  OPTION  ON  CHECK  DEPTH.  IF  =-1. PROGRAM  WILL 
PROVIDE  WAVE  RAY  INFORMATION  AT  TIME  STEP  NEAREST 
CHECK  DEPTH  AND  CONTINUE  PROCESSING  OF  RAY I  IF  =  0. 
PROGRAM  DOES  NOT  LOOK  FOR  CHECK  OEPTh;  ANO  IF  s  1. 
PROGRAM  PROVIDES  DESIRED  INFORMATION  AND  STOPS 
PROCESSING  OK  RAY. 

ALPHANUMERIC  TITLF  TO  APPEAR  ON  PRINTER  OUTPUT  AND  PLOT 
FOR  JOH  IDENTIFICATION. 

IF  =  0.  NEW  WAVE  INFO  TO  BF  READ 

IF  =  1.  WAVE  INFO  IS  RRt-VIOUSLY  READ  INPUT 


1 1-55 


Jt-’EN 


m 


LCK 


L IMNPT 
LPLOT 
PE  AS 
*  I 
MJ 

NOWAYS 

NOSETS 

NPRINT 

PD 


RK 
RK  1 
SCNV 

sex 

SK 
SKI 
ST  AZ 
T 

ThI 

UNIT 

X 

XSG 

XSLINE 

Y 

YF  A 


YSG 

YSLINE 

W 

WPI 


DEEP  WATER* 
INITIAL 
=  1  *  RAYS 


FOR  A  RAY. 
A  RAY. 


PEN  POSITION  WHEN  MOVING  TO  A  NEW  PLOT  POINT  ON 
SHORELINE . 

TELLS  PROGRAM  IF  SET  OF  RAYS  IS  STARTING  IN 
IF  =  0«  RAYS  ARE  STARTING  IN  DEEP  WATER  AND 
ANGLE  OF  RAYS  WILL  RE  DEEPWATER  AZIMUTH?  IF 
ARE  MOT  STARTING  IN  DEEP  WATER  AND  STARTING 
AZIMUTH  FOR  EACH  RAY  MUST  BE  INPUT. 

MAXIMUM  NUMBER  OF  TIME  STEPS  TO  bE  COMPUTED 
NUMBER  OF  TIME  STEPS  BETWEEN  PLOT  POINTS  ON 
UNITS  WHICH  DEPTHS  ARE  EXPRESSED  IN 
NUMBER  OF  GRID  POINTS  IN  X-DIRECTION. 

NUMhER  OF  GRID  POINTS  IN  Y-DIHFCT ION. 

NUMBER  OF  WAYS  IN  A  SET. 

NUMBER  OF  SETS  OF  RAYS  TO  RF  PROCESSED. 

PRINTING  INTERVAL  FOR  TIME  STEP  INFORMATION. 

MAXIMUM  DEPTH  AT  WHICH  THE  PRINTOUT  WILL  BEGIN. 

THIS  IS  A  MEANS  TO  DECREASE  THE  NUMBER  OF  LINES 
PRINTED  BY  PRINTING  ONLY  THE  INFORMATION  FROM 
SHALLOW  WATER.  (FEET) 

REFRACTION  COEFFICIENT  OF  WAVE  RAY  AT 
LOCATION. 

REFRACTION  COEFFICIENT  OF  WAVE  PAY  AT 
PREVIOUS  TO  STARTING  LOCATION. 

MAGNIFICATION  OF  WINDOW 
SCALE  FACTOR  OR  LENGTH  OF  A  SIDE  OF 
EXPRESSED  IN  INCHES. 

SHOALING  COEFFICIENT  FOR  FIRST  TIME 
DEEPWATER  SHOALING  COEFFICIENT 
DEEPWATER  AZIMUTH  OF  A  SET  OF  PAYS. 

WAVE  PERIOD  FOR  A  SET  OF  RAYS 

CLOCKWISE  ANGLE  BETWEEN  NORTH  ON  MAP  AND  Y-AXIS  OF  GPID. 
TIME  STEP 

FOP  A  PAY. 

OF  LOWER  LEFT  CORNFR  OF  A  WINDOW. 

MI  IF  WINDOW  IS  NOT  DESIRFD. 

FOR  DEFINING  A  POINT  ON  SHORELINE. 

FOR  A  RAY. 

TO  THE  Y-COORDINATE  OF  THF  STARTING  POINT  OF 
ORDER  TO  DECREASE  THt  AMOUNT  OF  DEEPWATER  REGION 
TRAVELS  OVFRf  THEREBY  DECREASING  THE  COMPUTATION 


STARTING 
THE  TIME  STEP 


A  GRID  SOUARE 


STEP 


X-COOROINATF. 
X-COORDINATE 
SET  EQUAL  TO 
X-CUORDINATE 
Y-COORDINATE 
FACTOR  AODEO 
A  RAY  IN 
THE  WAVE 


TIME. 

Y-COOHDINATE  OF  LOWER  LEFT  CORNER  OF  A  WINDOW. 
LEAVE  PLANK  IF  WINDOW  IS  NOT  DESIRED. 

Y-COORDINATE  FOR  DEFINING  A  POINT  ON  SHORELINE. 
WATER  LEVEL  AT  ZERO  CONTOUR 

NUMbFR  OF  WAVFJ  PERIODS  BETWEEN  TIC  MARKS  ON  A  PAY. 


n  o  o  o  n  r»  r> 


THIS  W I NDOW  WOULD  HE  USED  FOR  CASES  WHERE.  IT  IS 

oEsrwEo  to  study  a  particular  area  in  greater  detail 

8Y  USING  A  LARGER  SCALE  GRID.  THE  PROGRAM  WILL  TRACE 
THE  RAY  THROUGH  THE  WINDOW  ANU  PRINT  OUT  STARTING 
INPUT  PARAMETER  KOR  THE  RAY  ON  ThE  LARGER  SCALE  GPIO. 


COMMON  0£P (20000) »Ml«MJ,SCX.SCALE ♦GRIDtPD  : 

COMMON  PI*  B2,  HOUND*  CKDfcP.  CO*  CXY*  0(12) *  DCDH»  OCON. 

1DFLTAS*  DRC*  OTGR*  DXY*  F(t>),  GRINCt 
2H«  HH.  HO*  ICN*  IGO.  1  RE T »  ISP*  JGO*  KEIRST*  KPLOT * 

3LIMNPT*  LPLOT*  NRRINT.  NPT,  Phx*  PHY*  R90.  RADIUS*  PAYNO* 

4HCCO.  KHS.  R K*  SCNV.  SIG*  SK *  SKI.  T.  THI,  TOP. 

SV*  WL «  WLO.  XP.  XSG.  XSGL*  YP»  YSG*  YSuL 
COMMON  HPRFV*SKPRf V  .  RKPRF  V  . UBOT  » UM IP ♦  F. T  . C6R* POW.PL .  OS 
DIMENSION  IOATF(?) . ITITLF (10) *JBUF(S00O>. LEGEND  (15) 

DIMENSION  LEG!  (11).  LEG2  (6) *MFAS  (3) 

DIMENSION  XC (400) . YC (400) *AA2 (400) .PRKI (400) *RRK(400) 

DATA  HLE  G/0 . 1 1 2/ 

DATA  LEG1/»WAVE  •  ♦  •  PER  • «  '  I  OD=  *  *  *  »*•  SFC*.*S  •  *  'OEF'P'  *  'WATF*. 

1  *M  AZ»* ' IMUT». *H=  •/. 

2  LEG2/ *  OEG» . 'REES' » 'DATE* . »  PLO ' . » TTED • . *  :  •/ 

C  FOPMaT  STATEMENTS  IN  MAIN  PROGRAM 

80  FORMAT  (•  * . IS*  *  RAYS  EXCEEDS  ThE  MAX  LIMIT  OF  175') 

90  FORMAT  <•  «♦«  THE  Value  of  KPLOT  IS  TOO  Small*) 

100  FORMAT  ( 1 0 A4 ) 

110  FORMAT  (4I5.3F10.4) 

120  FOHMAT  (8b 10.5) 

130  FORMAT  (15) 

140  FORMAT  (FS.O) 

160  FORMAT  (2I5*bF10.2,Fl0.6) 

170  FORMAT  <215* JF10. 4*15. 2F10.3) 

171  FORMAT  (is) 

178  FORMAT  (1 A4*2X,3A4,?X*A4) 

180  FORMAT  (1H-.2A4) 

190  FORMAT  <•  ».///*42X,*W0VE  REFRACTION  ANALYSIS  *.10A4> 

200  FORMAT  Ml  •♦///////, 39X. 'WAVE  REFRACTION  PROGRAM — COASTAL** 

1  *  RESEARCH  DIVISION* ) 

210  FORMAT  (//,*  mi=» ,I4*13x*»MJsi,l4,13X* 'LIMNPT*'* 

1  J5.8X, *NPRINT=*. I4.9X*  *GRID=* ,F 10.4* bX. *OCON=« ,F 10.4/*  DELTAS** 

2  FS.3) 

220  FORMAT  (8H0SET  NO..I3.10H,  RAY  NO.*I3*10H*  PERIOO=*F7 .3, 

1  6H  StCS.) 

230  FORMAT  (8H  GRINC  =,F10.6*l4H*  TIME  STEP  =.F 8.3.7H  SECS.** 

1  24H  WAVE  FRONT  INCREMENT  =.F7.2.6H  SECS.) 

240  FORMA T (/ '  POINT  X*»SX»*Y**2X. 'AZIMUTH  DEPTH  LENGTH  SPEED  RK'*| 
1  SX» *SK* ,4X* 'HEIGHT  DMID  UHOT *♦ 4X .* TO TAL *. 6X *• GROUP  * «  i 

?  4X. ‘TOTAL  LONGSHORE*.'  DIR'./**  '» 79X* 'ENERGY  * *4X.  | 

3  *  VELOC I T Y ' *  3X  * »P0WFP*  *SX* 'ROWER •/♦ 

4  I6*F7.1 *2F6. 1 .GX.F7.1 *F6. 1. 14X.F7.2) 


28  0  FORMAT  ('1','MI  *  MJ  EXCEEDS  20000') 

27  0  FORMAT  (//,34rt  ALL  SETS  COMPLETED*  NO.  OE  SETS  =*I4) 

280  FORMAT  (////,'  ROUND=',K5.1,9X,'SCXs',F10.S,6X,'XSGs,,F6.2,10X, 

1  «YSG=' ,F6.2» 10X, 'SCNV=»,F6.3,9X,  *D6XL=«  »FS.2/»  DGYL=' »FS.2» 10X 

2  'SCALES' *F 10. 3*4X. 'NOSE TS=» ♦ 13) 

282  FORMAT  <////,*  LPLOT=' « I5*9X, »NORAYS=* , IS*«X» • T=' »F5.2» 13X, 'HOs* 

1  FS.2»12X»'SK='*F5.2»12X»*SK1='»FS.2»/'  THI=' ,F7.2,9X, 'STAZ=* ,F7 

2  *8X. 'UNIT  =  ' «F8. 2/'  I SP= * , 1 3 , 1 3X * • LCK= ' , 1 3 *  1 3X » • WP I  =  • ♦ F7 . 2 . 9X , 

3  'CKDEPs • *F6.2,6X, *DF  =  *  »F"5.2,  1?X»  » IWR= » ♦ 13*/'  PO=» .F8.2»9X, 

4  'YEAS', E5.1) 

284  FORMAT  (•  ' , / , 72X » 2 A4 ) 

286  FORMAT  (///////////////,»  INPUT  PARAMETERS  :•> 

288  FORMAT  (/////,'  MAP  OATES  », A4 , 1  OX »• DEPTH  IN  *,3A4) 

9980  FORMAT  (4A4) 

SINH(DUM)=.5*(EXP(0UM)-U./EXP(()UM)  )  ) 

C  READ  TITLE  INFORMATION 

READ  (S,100>  ITITLE 
RF  AD  ( S , 140)  FAC 
READ  (S.178)  DATE  1 »  MF  AS  »  W 
CALL  DATE  (IDATE) 

C  READ  GRID  INFORMATION  AND  OTHER  BASIC  DATA 

READ (S, 1 1 0)  MI,  MJ,  LIMNPT,  NPRINT,  GRID,  OCON,  DELTAS 
Umax  =  mi  *  mj 
IF  (  UMAX. GT. 20000)  GO  TO  290 
GO  TO  300 
290  WRITE ( 6 » 2S0 > 

STOP 

C  READ  PLOT  INFORMATION 

300  READ ( S  * 120)  HOUND,  SCX.  XSG,  YS6,  SCNV*  DGXL ,  DGYL 

C  READ  ALL  DEPTH  VALUES 

SCALE  =GR I D/SCX 

C  SET  ORIGIN  FOR  PLOTS 

CALL  PLOTS  ( JEUF,5000,9) 

IF  (FAC  .EQ.  0.)  F  AC  =  1 • 

CALL  FACTOR  (FAC) 

CALL  PLOT  (3. ,2., -3) 

CALL  BOTTOM 
WP I TF  (6,200) 

C  SET  SOME  INITIAL  VALUES 

RPS=MT 
RHS=RHS-2.S 
T0P=MJ 
TORsTOP-2.5 
XLIMIT=(MI-1)*SCX 
R90=1.S7 07 96327 

C  PRINT  TITLE  AMD  GRID  DATA 

WPITF (6, 190)  ITITLE 
WRITE  (6,264)  IDATE 
WRITE  (6,286) 

WRITE  (6,210)  MJ ,MJ, LIMNPT, NPRINT, GRID, DCON, DELTAS 


11-58 


C  READ  NUMHLH  OF  SETS  OF  RAYS  ANO  SHORELINE  DATA 

READ  130.NOSE7S 

WRITE  (6 » 280 )  ROUNt)»SCX*XSG»YSG,SCNV,DGXL»DGYL,SCALE, NOSETS 
360  CONTINUE 

C  PROCESS  EACH  SET 

XSLINF  =  0.0 
YSLINE  =  0.0 
DFFsO. 

3M0  00  S10  N0SET=1 .NOSE TS 

READ  ( 10,9990»END  =  4l5>  XSS,  YSS,  JPFN 
9990  FOkMAT ( 3A4 ) 

XSLINF  =XSS/GR 1 0 

yslinf;=yss/grid 

C  PLOT  LOWER  LEFT  CORNER  F OR  A  SET 

CALL  PLOT  (0.0,  1.0,  3) 

CALL  PLOMO.,  1.,  ?_) 

CALL  PLOT (0. •  0.,  2) 

CALL  PLOT ( l • *  0.,  2) 

CALL  pLOT (0. ,  0.,  2) 

CALL  PLOT  ( 0  «  «  0 . « 3 ) 

XSL=XSLINE*SCX 
YSL=YSLlNfc*SCX-HOUNO 
CALL  PLOT  <XSt»YSL»3> 

GO  TO  411 

C  PLOT  ShORLF.  INF  FOR  A  SET 

410  CONTINUE 

READ  ( 10,9990,END=415)  XSS,  YSS*  JPFN 

411  LPEMSJP EN 
XHOLD  =  XSLINF 
YHOLO  =  YSLINE 
XSLINF=XSS/GR10 
YSL1NE=YSS/GR10 

IF  (JPFN  -  1)  400,390,400 
3R0  LPEN=? 

CALL  RfcFOIA  (XSLTNE,  YSLINE,  2) 

CALL  R E F  0 1  A  { X HOLD, YHOLO, 2 > 

4 0 U  CALL  PFF01A  (XSLTNE, YSLINE, LPFN) 

GO  TO  410 
41S  CONTINUE 
RE  w I NO  10 

C  PLOT  LOWER  LEFT  CORNER  OF  WINDOW 

420  IF  (XSG-RHS)  430,440,440 
430  CALL  PEFD1 A (aSG,YSG, 3) 

CALL  REED  I A ( XS6, YSG* 1 • ,2) 

CALL  PFFDIAO.SG,YSG*?> 

CALL  REFOIA (XSG*1.,YSG»2) 

C  plot  lower  right  CORNER  OF  WINDOW 

CALL  REFDTA(XSG,YSG,2) 

XSbL  =  XSG*DGXL/SCX 
YSGL=YSG*DGYL/$CX 


i 


CALL  PEFD1A(XSGL,Y5G,3) 

CALL  WEFDIA(XSGL.YSG*1.,2> 

CALL  PEEDIA(XSGL»Y5G,2) 

CALL  PEF  0 1 A  <  XSGL- 1..YSG.2) 

CALL  PEFD1 A ( XSGL  » YSG, 2) 

READ  BASIC  WAVE  BATA  FOR  A  SET 
440  WEAD (5, 160)  LPLOT,  NOWAYS*  T*  HO*  SK*  SKI*  THI *  STAZ. 

1  UNIT 

IF  (NOWAYS  .GE.  400)  GO  TO  800 
SKPH=SK 

READ  (5.170)  ISP.LCK*WPI,CKDEP*DF, IWW.P0*YFA 
IF  (IWW  .FQ.  0)  READ  (5.171)  I  WED 
IF  ( 1 WED  .NE.  25)  IREO=5 
IF  <P()  .EO.  0.)  PO  =  (JFP  ( 1 )  *DCON*DF 

IF  (NOSET  .EU.  1)  WRITE (6.28?)  LPLOT .NOW AYS. T .HO ,SK . SK 1 , THI , 
*  STAZ.UNIT.ISP.LCK.WPI.CKDEP.DF . IWP.PU.YFA 
IF  (NOSET  .EO.  1)  WRITE ( 6 . 288 )  DATE ) .MtAS 
IF  (NOSET  .EO.  1)  GO  TO  445 
WRITE  (6,200) 

WRITE (6. 190)  ITITLE 
WRITE  (6.284)  IDATE 


WRITE (6. 190) 
WRITE  (6,284) 
WRITE  (6,286) 
WRITE  (6,210) 
WRITE  (6,280) 


WRITE  (6,210)  MI. MJ.LIMNPT.NPRINT.GWID.DCON. DELTAS 

WRITE  (6,280)  HOUND, SCX. XSG,YSG,SCNV,DGXL»DGYL, SCALE. NOSETS 

WRITE (6.282)  LPLOT . NOR AYS. T , HO. SK ,SK 1 . THI . 

*  STAZ.UNIT. ISP.LCK.WPI ,CKOEP,DF ♦ IWP.PO, YFA 
wP I TE ( 6 , 288 )  DATE1.MEAS 
445  CONTINUE 

SET  INITIAL  VALUES  FOR  THE  SET 
IF  ((NOSET  .GE.  2)  .AND.  (DF  .EO.  D FT))  GO  TO  340 
IF  (NOSET  .Gt.  2)  DF=DF-OFF 
IF  (OF)  310,340.310 
310  00  330  J  =  l.IJMAX 

DEP  (J)  =  DEP  (J)  ♦  OF 
330  CONTINUE 
340  OFF=DFF+DF 
KSP=ISP 

SIG=6.28318531/T 

C0=5.120406*T 

WLO=CO*T 

DWC=WLO*0.5 

DT  GW=UN1 T /GRID 

GRINC=DTGP*C0 

WFI=T*WPI 

KPL0T=WF I/UNIT*0.1 
IF  (KF’LOT  .EO.  0)  GO  TO  810 
PROCESS  EACH  RAY  IN  SET 
DO  500  NORAY=l .NOWAYS 
HPWE  V=HO 
nS=GRINC 
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j 
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KFIRST=0 

WRITE ( 6  »  1 HO ) IOATF 
ISP=KSP 

WRITE  (6.220)  NOSET  *  NORAY*  T 
IF  (LCK)  460.450*460 
C  READ  RAY  DATA 

450  IF ( IwR  .EQ.  0)  RFAO ( I RED* 120*EN0=501 )  XC (NORAY) *YC (NORAY) 
IF  (XC(NORAY)  .LF.  0.)  00  TO  501 
X=XC (NORAY) 

Y=YC (MORAY) ♦  YEA 
A7IRTH=STAZ 
RK=  1,0 
SK  =  1  , 

81*1.0 
82  =  1.0 
RKPPEV=1 . 

5KPPE V=1 . 

GO  TO  470 

460  IF ( I WR  .FO.  0)  HEAO<TRI:0.l20.END=501>  XC (NORAY ) .YC (NORAY ) * 
*  AAZ (MORAY) »RRK (NORAY) «  RRK 1 (NOWAY) 

IF  ( XC ( NOWAY )  .LF.  0.)  GO  TO  501 
X=XC (NORAY) 

Y  =  YC  <  MORAY ) ♦YFA 
a7InTh=aaz (Moray ) 

RK  1  =RRK  1  (NORAY) 

RK=RHK (MORAY) 

H1=1./RK1/RK1 
H2=l ./RK/RK 
SK=S*PR 
SKPHfcV=l. 

HKPRK  V= 1 • 

470  A  =  THW70.0-AZIMTFi 
C  SET  INITIAL  VALUES  FOR  RAY 

NPT  =  1 
CXY=CO 
WL=WLO 

C  PRINT  STEP  AND  T I WE  INFORMATION 

WRITE (6.230)  GRINC.  UNIT.  WFI 
WAYNU=N0RAY 
XNO=X*SCX 
YNU=Y*SCX-80UN0 
IF  ( YNO  >  4RO.4MO.4flO 
C  IDENTIFY  RAY  ON  PLOT 

480  KFIRST=l 

CALL  NUM8ER  (XNO.YNO-O. 15.0.1 12. WAYNO.O..-1) 

CALL  REFOIA  (X.Y.3) 

C  PRINT  DEEPWATER  WAVE  INFORMATION 

4R0  WHITE (6.240)  NPT.  X.  Y.  STAZ*  WLO.  CO.  HO 
C  CALL  RAYCON  TO  COMPUTE .PRINT. AND  PLOT  THIS  RAY 

CALL  RAYCON  (X.Y.A) 


I 


500  CONTINUE 

501  CONTINUE 

C  PRINT  TITLE  AND  WAVE  DATA  ON  PLOT 

XLEG=XLlMIT-b.72 

CALL  SYMBOL  < XLEG «-0 .5. 0 . lbft, 1  TITLE . 0. .40 ) 

CALL  SYMBOL  (XLEG*  - 1 . *HLEG.LEG1 ( 1 >  *  0.*14) 

CALL  NUMREW  (999. *999 . . 0. 1 1 2, T* 0. « 3) 

CALL  SYMBOL  (999. * 999 . .HLEG.LEG1 (5)  .0..5) 

CALL  SYMBOL  (XLEG ♦ -1 . 25. HLEG.LEG1 ( 7 >  *0.*20> 

CALL  NUMBER  ( 999. « 999 . . 0 . 1 1 2.STAZ » 0. . 1 ) 

CALL  SYMBOL  (999. *999. .HLEG.LEG2 ( 1 )  ,0..B) 

CALL  SYMBOL  (XLEG.- 1 . 50.HLEG.LEG? ( 3)  .0..15) 

CALL  SYMBOL  ( 99®. . 999 ...112*1  DATE . 0 • « B) 

CALL  SYMBOL  <XLEG*4. ,-l . .HLEG. *  0  CONTOUR  AT  *.0..13) 

CALL  SYMBOL  ( 999. ♦ 999 . . HLEG* W . 0 • * 4 ) 

CALL  SYMBOL  ( XLEG*4 . « - 1 . 25. HLE G. • PEL AT IVE  WATER  LEVFL=  *,0.,??> 
CALL  NUMBER  ( 999. « 999 . .HLEG . DEE . 0 . , 2) 

CALL  SYMBOL  < XLEG* 7 . . - 1 . 25. HLEG. ME AS. 0. « 12) 

CALL  SYMBOL  (XLEG*4. .-1 .50. HLEG. 'MAP  DATE:  *.0,.10) 

CALL  SYMBOL  (999.«999.»HLEG«DATE1*0.*4) 

CALL  SYMBOL (XLEG. -1. 75. HLEG. ‘DEEPWATER  WAVE  HEIGHT  =  *.u..23> 
CALL  NUMBER (999.. 999.. HLEG. HO. 0.. 2) 

CALL  SYMBOL (999.. 999.. HLEG.*  FEET*. 0.. 5) 

C  OPAW  LOWER  RIGHT  CORNER  OF  PLOT  FOR  A  SET 

CALL  PLOT  (XLIMIT.O. .3) 

CALL  PLOT  (XLIMIT.1,.2) 

CALL  PLOT  (XLIMIT.O. .2) 


511 

509 


600 


(XL1MIT  -1..0..2) 
(XLIMIT.O. .2) 
(XLIMIT.O. .3) 


XS.YS. DEPTH, I PN 


509 


CALL  PLOT 
CALL  PLOT 
CALL  PLOT 
1ST  =1 
I  ST  1 =0 
REMIND  11 
CONTINUE 

READ  ( 11 ,9980,FND  =  790) 

XX=XS*SCX/G«ID 
YY=YS*SCX/GRID 
IF  (YY  ,LT.  dOUNO)  GO  TO 
YY=YY— ROUND 

hdepth^depth 
IPEN=3 
GO  TO  730 

READ  (11.9980. ENO=770)  XS. YS. DEPTH, IPN 
IF  (DEPTH  .LE.  0.)  GO  TO  bOO 
IF  ( I PH  .FQ.  IPN)  GO  TO  650 
IPH  =  2 

HOEPTH=DEPTH 

CALL  PLOT  (XH.YK.IPFN) 

IPEN=3 
IST  =  1 
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1ST  1  =  0 

CALL  PLOT  <XH.YH,IP£N> 

XX=XS*SCX/GP1H 

YY=YS*SCX/GRID 

IF  (YY  .LT.  GOUNO )  GO  TO  665 

YY  =  YY-BOUNl) 

IF  { 1ST  1  .tO.  0)  GO  TO  670 


IF  (XX  .LT.  XXH)  GO  TO  660 
TF  (XX  .GT.  XXH*SIZE>  GO  TO  660 
IF  (YY  .GT.  YYH* .06)  GO  TO  660 
IF  (YY  .LT.  YYH-.06)  GO  TO  660 
IPEN=3 
GO  TO  600 
CONTINUE 
1ST  =0 
GO  TO  660 
IPH  =  3 
GO  TO  600 

IF  (1ST  ,FQ.  0)  GO  TO  660 
1ST  1  =  1 
IST  =  0 

aopth=-oept  H 

DPTH=ABS (DEPTH) 

SIZE=3. 

IF  <l,PTH  .GF.  1.)  S I ZE=S I Zt ♦ 1 • 

I F ( DPTH  .GE.  10.)  S I/E=S I ZE ♦ 1 • 

I F ( OPT H  .GE •  100.)  SIZE=SIZE*1. 

IF (OPTH  .GE.  1000.)  S I ZE=SI ZE ♦ 1 . 

TF  ( DEPTH  .GT.  0.)  STZE=SIZE»1. 

SIZE =S IZE *  .09 

IF  (XX  .GE.  FLOAT (MI-3)*SCX)  GO  TO  675 

CALL  NUMBER  (XX^.035»YY-.035».10*ADPTH«O.*l) 

XXH  =  X  X 

Y  Y  H  =  Y  Y 

GO  TO  600 

CONTINUE 

XXX=XX-S I ZE 

CALL  NUMBER  (XXX.YY-.035».10.ADPTH»0.*1> 
XXH=XXX 

yyh=yy 

GO  TO  600 
CONTINUE 

CALL  PLOT  ( X X * YY » IPEN) 

TPFN=? 

XH  =  XX 
YH  =  YY 
GO  TU  600 

CALL  PLOT  ( X X ♦ Y Y ♦ 3 ) 

CONTINUE 

SET  ORIGIN  FOR  A  NEW  PLOT 


1.)  S I ZE=S I Zt ♦ 1 • 
10.)  S I/E=S I ZE ♦ 1 • 
100.)  SIZE=SIZE*1. 
1000.)  S I ZE=SI ZE ♦ 1 • 
.  0.)  STZE=SIZE*1. 


( XXX. YY-. 035 ♦.lO.AOPTHtO.*!) 


( X  X  « YY  « I PEN) 
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CALL  PLOT  <XLIMIT*3..0.»-3) 
510  CONTINUE 

WRITE(h,270>  NOSFTS 
GO  TO  a?o 

800  WHITE (6,80)  NORAYS 
GO  TO  820 
810  WHITE  (8. HO) 

820  CALL  PLOT  ( 0 . . 0 . « 9HH) 

STOP 

END 
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o  n  n  do 


RAYCON  SUBROUTINE  DECK 

controls  computations, pointing, and  plotting  of  each  pay 

SUBROUTINE  PAYCON  (X.Y,  A) 

COMMON  DEP ( 20000 ) ,MI .mJ.SCX.SCALF.oRID.PD 
COMMON  H 1 «  H2.  ROUND.  CKOEP,  CO.  CXY.  0(12) »  DCDH.  DCON. 
10ELTAS,  DPC.  UTGR.  DXY.  F (6) ♦  Gu INC  » 

2  H,  HH,  HO.  ICN.  1GO.  IPfc'T .  ISP.  JGO.  NF 1  PST .  KPLOT * 

3  LIMNPT.  LPLOT.  NPPINT.  NPT .  PHX.  PHY.  PRO.  RADIUS.  PAYNO, 
4PCC0.  PHS.  H K.  SCNV.  SIG.  SK.  SKI.  T.  THI*  TOP. 

S  V.  Wt  .  WLO.  XP.  XSG,  XSGL.  YP.  YSG,  YSGL 
COMMON  HPPE V.SKPRfcV.PKPPFV.UrOT.UM I D. E  T. CGH.POW.PL. OS 
COMMON  ] DIP 
0 1 MENS  T ON  p  t 1 6 ) 

100  FORMAT  ( •  • .IS. .4F6.1.F7. 1 ,F6.  1  ,2F7.  3. F7.2.2F6.2.4F  i  0.  ' 
110  FORMAT  ( *  •  .  ls.F7.U3F6. 1  .f  7.1  ,F6.  1 .2F7.3.F 7 .2.2F6.2, 4F  1 0 .2. 
120  FORMAT  (2PH  CURVATURE  AVE’RAGE  U  AT  POINT.  I A ) 

130  FORMAT  (34H  RAY  INSIDE  DETAILED  GRID  AT  POINT. 14, HH  OFT  A ] 
1  29H  GRID  STARTING  LOCATION  IS  X  =  .F7.?.4h  Y=.F7.2) 

140  FORMAT  (/.42H  PAY  STOPPED.  NO  CONVEPGFnCE  FOR  CURVATURE) 

ISO  FORMAT  ( / , 31n  RAY  STOPPED,  WAVE  BREAKS  AT  X=,F7.2.4H  Ys.Fl 


3X.AI) 

.3X.AI) 


140  FORMAT 

iso  format 

160  FORMAT 
170  FORMAT 
180  FORMAT 

l y o  format 


</.33H  RAY  STOPPED. 
(/,32m  RAY  STOPPED* 
</,34h  RAY  STOPPED. 


REACHED  BOUNDARY  X  =  .F7.2.4H  Y  = 
LIMNPT  EXCEtPFO  X=,F7.2.4H  Y=. 
OELT A »S  TOO  Small  X  =  .F7.?,4H  y 


7.2) 

.F7.2) 

F7.2) 

=  ,F7.2) 


( 9H  RAY  NEAR.  F6.2.14H  FOOT  CONTOUR  ) 


IF  THF  RAY  STARTS  IN  w ATEW  LESS  THAN  HREAKER  DEPTH  GO  TO  THE  N£XT 


ANtX=A*.0174s3?925 
XNEX=COS ( A) *GPINC*X 
YNFX=SIN (A)*GRINC*Y 
CALL  DEPT h ( XNF  X ♦ YNE X ) 

IF  UH*DXY-HO)  492.492,494 
492  WRITE  (6,289) 

289  FORMAT  (‘O', 'RAY  STARTS  TOO  CLOSE  T< 
RETURN 
494  CONTINUE 
END=0. 

IF  ( X  Sti-kHS )  200.210,210 
200  LFIRST=0 

SGX=XSG*1 .S/SCNV 
SGY=YSG. 1 .S/SCNV 
SGXXsXSGL-1 .S/SCNV 
SGYYsYSGL-l. S/SCNV 
GO  TO  220 
210  LFIRST=1 

SET  INITIAL  VALUES  USED  TO  START 
220  LPh InT=0 


SnORE*) 
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FLAG 1=0.0 
ANG=THI*270.0-A 
A=A*U.O 174532925 
C0SA=C0S  (A) 

SINA=STNCA) 

H  =  h0 
IG0=1 

C  SAVE  VALVES  OF  X  AND  Y 

230  PX=X 
PY=Y 

C  ADVANCE  X  AND  Y  ONE  STEP  AND  FIND  NEW  DEPTH 

X=COSA*bHlNC*X 
Y=SINA*GwINC*Y 
CALL  DEPTH ( X  »  Y ) 

NWR I T  E  =  1 

C  CHFCK  TO  SEE  IE  WAVE  HEIGHT  IS  GPFATFR  Than  .«  OF  WATEP  nEPTH 

If  <  .E*I)XY-H>  380,380*240 

C  NO,  CHECK  TO  SEE  IF  RAY  IS  IN  DFFP  WATER 

240  IF(DXY-OPC)  320,320, 2b0 

C  YES,  ADVANCE  STEP  COUNT  AN[)  CHECK  TO  SEE  IF  LIMNPT  IS  EXCEEDED 

250  NPT  =  NPT ♦  I 

CALL  FRICTN 

IF (LIMNPT-NPT)  400,400,260 

C  MO,  CHECK  TO  SFE  IF  PAY  IS  TOO  CLOSE  TO  HOUNDAWY 

260  IF(RHS-X)  270,270,2*0 

C  YES.  SET  TO  PPlNT  OUT  PAY  TOO  CLOSE  TO  POUNDAPY 

270  I  PET  =5 

GO  TO  350 

280  IF (X-l .5)  270,270,200 
200  IF(TOP-Y)  270,270.300 
300  IE(Y-l.S)  270.270,310 

C  PAY  WITHIN  riOUNDARY  AND  IN  DEEP  WATtR,  GO  TO 

C  CHECK  PRINTING  AND  PLOTTING  OPTIONS 

310  GO  TO  480 

C  RAY  WITHIN  BOUND ARY  BUT  NOT  IN  DEEP  WATER 

320  X=PX 
Y  =  P  Y 

C  COMPUTE  CURVATURE  OE  RAY  IN  FIRST  STFP  AFTER  DEEP  WATEP 

CALL  CURVt (X.Y.A.FK) 

C  ADVANCE  STEP  COUNTER  AND  CHECK  TO  SEE  IF  LIF'NPT  IS  EXCEEDED 

330  NPT  =NPT ♦ 1 

IF (NHT-LIMNPT)  340,400,400 
340  NWR I TF  =  1 

C  COMPUTE  X  ANO  Y  COORDINATES  AND  ANGLF  OF  RAY  FOR  NEW  STEP 

FLAG1  =  1 .0 

CALL  RFFRAC (X»Y. A,FK) 

C  GO  TO  CONDITION  CONTROLLED  BY  VALUE  OF  IRET 

350  GO  TO  (410, 3»0, 370, 360»3R0, 420) , IRET 

C  SET  NwPITt  FOP  STATUS  OF  RAY 

360  NWRI TF=2 


\ 


GO  TO  42 0 
370  N*H I TF  =  3 
GO  TO  420 
380  NWPITE=4 

IF  (Fl.AGl  .EU.  0.0)  GO  TO  492 
GO  TO  420 
390  NWPTTF=5 
GO  TO  730 
400  NWR I T F  =6 
GO  TO  730 
410  NWPITF=7 
GO  TO  730 

C  COMPUTE  REFRACTION  AND  SHOALING  COEFFICIENTS  AND  WAVF  HEIGHT 

420  CALL  Hf IGHT (A) 

IF (DXY-CKOEP)  430.430,460 
430  IE(ISP)  450.470,440 
440  I  GO  =  3 
450  ISP=0 

N  W  P  I  T  F  =  9 
GO  TO  470 

460  IF(ISP)  470  .<+70.740 
470  IF(NwPITF-l)  460,480.490 

C  CHFCK  PPInT  OPTION  AM)  SF  T 

480  IF (MOO(NPT.NPPFNT) )  500.490.500 
490  LPPIM  =  1 

C  HflS  NUMBER  HEFN  PLACED  ON  TnIS  PAY 

500  IFIKFIPST)  530,510,530 

C  IS  THIS  STEP  WITHIN  STARTING  PLOT  PUUNOAPY) 

510  YON=Y*SCX-HOuND 

IF (YON)  700,520,520 

C  YES,  PLOT  THE  PAY  NUMPF.P,  SET  TO  NO  AND  PEN  UP,  AND  GO  TO  PRINT 

520  CALL  MUMmEP<X*SCX,YON-0. 15,0. 112. PAYNO,0.,-1) 

KF I  PST  =  1 

CALL  PFFD1A  ( X , Y , 3 ) 

C  CHFCK  TO  SEE  IF  PAY  PFaCHF.D  DETAILED  GP ID  ON  ANY  PREVIOUS  STEP 

530  IF  (LFIPST)  590,540,590 

C  NO,  IS  »AY  WITHIN  DETAILED  GRID  ON  THIS  STEP 

540  IF(X-SGX)  590,550,550 
550  IF  <  X-SGXX)  560,580,550 
660  IF(Y-SGY)  590,570,570 
570  IF (Y-SGYY)  5«0,580,S90 

C  YES,  SET  nWPITE  FOP  THIS  CONDITION 

560  nwkITf=P 
LFIPST=1 
L  P  P  I N  T  =  1 

C  COMPUTE  DETAILED  GRID  STARTING  POSITION 

l)GX=  (X-XSG)  *SCN V 
OGY  =  (Y-YSG)  *SCNV 

STEP  WITHIN  MOUNDAPY,  IS  TIC  MARK  Tu  PF  PLOTTED  ON  PAY 
590  I F  < MOD ( NPT ,F  PLOT ) )  600.620,600 


C 


n  n 


C 

c 


c 


c 


c 


c 

c 

c 


c 


c 


c 


NO.  IS  STEM  TO  RE  PLOTTED 
*00  IE (MOO (NPT.LPLOT) )  700. *10. 700 

YES,  PLOT  STEP  AND  GO  TO  PRINT 
610  CALL  PFFOIAtX, Y.2) 

GO  TO  700 

TIC  MARK  TO  RE  PLOTTED  PERPENDICULAR  TO  WAY, 

COMPUTE  ANGLE  AND  AZIMUTH  IN  OEGREES 
620  AA=A*57. 29576 

AN’G  =  THI  ♦2  70.0-AA 

SET  TO  PLACE  ASTERISK  ON  PRINT  OUT  dY  THIS  STEP 
LPP INT=-1 

CALL  PEF01A(X,Y,2) 

FIND  OUADRANT  THAT  RAY  IS  IN  AND  COMPUTE  TIC  MARK  COORDINATES 


IE 

<  A A -VO. ) 

630,630.640 

630 

YT 

10=0.05* 

SIN (R90-A) 

XT 

I C=0 • 05* 

S  I N  <  A ) 

GO 

TO  690 

640 

IE 

(AA-180. 

)  o50,650,660 

650 

YT 

IC=-0 • 05 

*S IN ( A-R90 ) 

XT 

I C=0 . 05* 

SlN(3. 141592684- 

A) 

GO 

TO  690 

660 

IF 

( AA-270. 

)  o7  0 ♦ 67  0 , 680 

670 

YT 

IC=0 . 05* 

SlN(A-3. 141592654) 

XT 

1C=0.05* 

SIN (4. 7123889*1- 

A) 

GO 

TO  690 

68  0 

YT 

1 C=-0 . 05 

*SIM(A-4. 712388961) 

XT 

1 C=0 . 05* 

SlN(6. 283185308- 

A) 

PLOT  T I 

C  MARKS 

690 

CA 

LI  REED  I 

A  ( X ♦XT  1 C  »  Y-YT I C 

,2) 

C6 

LL  REFOI 

A(X-XTIC,Y*YT JC. 

2) 

CALL  PEEDIA (X.Y.?) 

CHECK  PRINT  OPTION 
700  I E ( LPP I NT )  710,740,720 

PRINT  STEP  INFORMATION  WITH  ASTERISK  BY  STEP  NUMBER  TO  DENOTE 

THAT  A  TIC  MARK  WAS  MADE  PERPENDICULAR  TO  WAY  ON  PLOT 
710  CONTINUE 

CALL  ENERGY  (A, END) 

IF  (DXY  .LE.  PO) 

F.  WRITE  (6,100)  NPT,X,Y,  ANG,DXY,WL,CXY,RK,SK,H,UPin»UROT.ET» 

X  CGR,POW,PL» 1DIR 

RESET  PRINT  OPTIUN  TO  NO  PRINT  AND  GO  TO  CHECK  STATUS  OF  RAY 
LPR I nT=0 
GO  TO  740 

COMPUTE  PAY  AZIMUTH  IN  DEGREES  AND  PRINT  STEP  INFORMATION 
720  ANG=THW70.O-A*57.?RS7P 
CALL  ENERGY  (A, END) 

IF  (DXY  .LE •  PD) 

N  WR I TE (6,110)  NPT,X,Y,ANG,UXY»wL,CXY.WK.SK»H»UMTD.UPOT,ET. 

K  CGR.POW.RL. IOIW 

RESET  PRINT  OPTION  TO  NO  PRINT  AND  00  TO  CHECK  STATUS  OF  RAY 
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LPP INT=0 
GO  TO  7M) 

LAST  bTKP  OF  A  way 

COMPUTE  WAY  A7 1 MIJTH  IN  DEGREES  AM)  PRINT  STEP  INFORM  A  T I  ON 
ANG=TH  I  ♦  i>70  .  U-A>»5  /  .  ^9b/8 
WP1TE (6,110)  NPT  •  >  ,  Y  »  ANG,UX Y 
PLOT  LAST  STFP  OF  PAY 
CALL  PEFI)IA(A,Y,?) 

GO  TO  COM)  I T I  ON  CONTROLLED  BY  VALUE  OF  NWRITF 
GO  TO  (830,750*760,  770,780,790, 800,820,850)  * 

NwPITE 

WWINT  PROGRAM  MESSAGE  AND  CONTTNUE  WAY  COMPUTATIONS 


PRINT  PROGRAM  MESSAGE  ANO  CONTTNUE  WAY  COMPUTATIONS 
WRITE  <6, 120)  OPT 
GO  TO  830 

PRINT  PROGRAM  ERROR  MESSAGE  AND  STOP  PAY  COMPUTATIONS 
WRITE.  <b.  1*0) 

GO  TO  810 

WRITE <b, 150)  X  «  Y 

CALL  PEFUIA  (  X  «  Y  «  2) 

END=1 . 

CALL  ENERGY (A, END) 

FNU  =  (J. 

GO  TO  810 

wP 1 TE  (6. 160)  X  «Y 

GO  TO  610 

WRITE ( 6 1 1 70 )  X,Y 

GO  TO  610 

WRITE  ( 6  » loO)  X  »  Y 

I  GO  =3 

GO  TO  830 

WRITE (6,130)  NPT , OGX ♦ UGY 
I  «»R  n  =7 

WRITE  (IWRlT,13b)  OGX  «06Y*ANG«RK,  RKPPF  \l 


NPT , OGX ♦ UGY 


WKlfE  (IWPI  T,13b)  OGX  ,0GY,  ANG«RK»RKPPF.V 
135  FORMAT  <3F1Q»2»2E'10,4) 

CfE"CK  status  of  pay  AND  WPANCH  to  PKOPFW  CONDITION 
030  GO  to  (230*330,840) , IGO 
Rao  RF  TUREJ 

850  WRITE (6«I90)  CKDFP 
GO  TO  h30 
END 
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RF.FRAC  SUBROUTINE  OECK 

COMPUTES  CURVATURE  OF  RAY  IN  A  NEW  STEP 
CHECKS  TO  SEE  IF  RAY  HAS  REATCHEO  ROUNDAPY 
SUBROUTINE  REFRAC (X,Y, A,FK> 

COMMON  0£P (20000) ,MI ,MJ,SCX,SCALE*GRID,PD 
COMMON  Hl»  B2,  BOUNOt  CKOER*  COt  CXYt  0(12) •  OCDHt  OCON, 
lOELTASt  ORCt  OT 6R »  DXY «  E(b>.  GRINC, 

2  Ht  HH,  HOt  ICNt  I  GO*  IRETt  ISP.  J60.  KFIRST.  KPLOT. 

3  LIMNPT.  LPLOT.  NPRINT.  NPT .  PHX,  PHY.  R90,  RADIUS.  RAYNOt 
4RCC0.  RHSt  RK.  SCNV.  SIG,  SK.  SKI,  T.  THI.  TOP, 

5  V,  WL,  wLO,  XP,  XSG.  XSGL,  YP,  YSG,  YSGL 

COMMON  HPREV.SKPREV,RKPREV,UROT.UMin,ET,CGR,POW,PL.DS 
DIMENSION  RF  <  9 ) 

C  SET  INITIAL  VALUE  FOR  NO  AVERAGE  OF  CURVATURE 

NCUR= 1 

C  GO  TO  CONDITION  CONTROLLED  BY  VALUE  OF  160 

GO  TO  (100,1 I0.3S0) , I GO 

C  FIPST  TIME  THROUGH  SUBROUTINE,  SET  FOR  THIS  CONDITION 

100  FKM  =  F  K 
IG0=2 

C  COMPUTE  STEP  LENGHT 

110  DS=CXY<*DTGR 

C  IS  STEP  LENGTH  SMALLER  THAN  PRESCRIBED  MINIMUM  VALUE 

IF (DS-DELTAS)  120,130,130 

C  YES,  SFT  IRET  FOR  THIS  CONDITION  ANU  RETURN  TO  RAYCON 

120  1PET=1 
RETURN 

C  NO.  COMPUTE  TEST  FOR  CONVERGENCE  OF  CURVATURF 

130  RESMAX=0.00005/DS 

C  DO  A  maXIMUN  OF  20  ITERATIONS  TO  COMPUTE  NEW  STEP 

1 A 0  DO  210  1  =  1,20 

C  COMPUTE  DELTA  ANGLE  FOR  NEW  STEP  TRIAL 

OELA=FKM*OS 

C  COMPUTE  TRIAL  ANGLE  ANO  COORDINATES 

AA=A*DELA 
AM=OELA*O.S,A 
XX=COS (AM) *OS*X 
YY=SIN(AM)*DS*Y 

C  COMPUTE  CURVATURE  OF  RAY  FOR  TRIAL 

CALL  CURVE (XX, YY,AA,FKK) 

C  WILL  WAVE  BREAK  IN  TRIAL  STEP 

if(.b*dxy-h>  isonsomso 

C  YES,  SET  [RET  FOR  THIS  CONDITION  AND  RETURN  TO  RAYCON 

150  IRET=? 

RETURN 

C  NO,  WAS  CURVATURE  AVERAGED  IN  LAST  SOLUTION 

160  GO  TO  (170,240),  NCUR 


i 


a 


'00,,  ■  pt . ' 
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no  o  no  on 


c 


NO «  COMPUTE  AVERAGE  CURVATURE  OF  THIS  TRIAL  STEP 
170  EKM  = <FK*FKK) *0.5 
C  IS  THIS  FIRST  ITERATION 

IF  ( 1 -  1 )  210*210* 1P0 

C  NO*  IS  CONVERGENCE  TEST  SATISFIED 

1«0  IF  (RESMAX-AfaS(FKP-FKM)  >190,  ISO, 24  0 
C  NO,  IS  THIS  THE  18TH  ITERATION 

190  IF (  1-18)  210*200,210 

C  YES,  SET  OUMMY  VARIABLE  TO  USE  IN  AVERAGING  PAY  CURVATURE 

200  FKCK=FKM 

CONVERGENCE  TEST  NOT  SATISFIED 

SAVE  CURVATURE  FOR  THIS  TRIAL  AND  IF  NOT  20TH  START  NEW  TRIAL 
210  FKP=FKM 

CONVERGENCE  TEST  NOT  SATISFIED  AFTER  20  ITERATIONS 
IS  CONVERGENCE  TEST  SATISFIED  BETWEEN  18TH  AND  20TH  ITERATIONS 
IF (RESMAX-AHS (FKCK-FKM) )  220,220,230 

NO,  SET  I RFT  FOR  NO  CONVFRGENCE  AND  RETURN  TO  RAYCON 
220  I RF  T  =3 
RETURN 

YES,  SET  CURVATURE  F.UOAL  TO  AVERAGE  OF  18Th  AND  20TH 
ITERATION  AND  SET 
230  FKM=(FKM*FKCK)*0.5 
NCUR=? 

C  CURVATURE  IS  COMPUTED,  SET  VALUES  OF  NEW  STEP  TO  RETURN 

GO  TO  140 
240  X=XX 


Y  =  YY 
A  =  A  A 
FKsFKK 


IF ( NCUR-2 ) 

280,250,260 

250 

IRFT=4 

RETURN 

260 

IF (RHS-Xl 

270,270,280 

270 

IRET=5 

RF  TURN 

280 

IF (X-l.S) 

290,290,300 

290 

IRE  T=5 
RETURN 

300 

IF (TOP-Y) 

310,310,320 

310 

I RF T=5 

RF  TURN 

320 

IF(Y-l.S) 

330,330,340 

330 

IPET=5 

RETURN 

34  0 

IRET=6 

3S0 

RETURN 

END 

■ 
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C  CURVE  SUBROUTINE  DECK 

SUBROUTINE  CURVE(X.Y.A.FK) 

COMMON  DEP (20000) ,MI *MJ*SCX »SCALE *G«ID*PP 
COMMON  B1 »B2* BOUND*  CKDEP*  CO*  CXY*  0(i?)*  PCDH.  OCON* 

10ELTAS*  ORC*  OTGR «  UXY.  E<6>.  6RINC* 

2  H*  HH»  hO*  ICN*  1  GO*  IRET.  ISP*  JGO*  KF1RST*  KPLOT * 

3  LIMNPT.  LPLOT*  NPRJMT*  NPT  *  PHX*  PHY*  R90.  RAOIUS*  BAYNO. 

4RCCO*  RHS*  PK*  SCNV.  SIG.  SK*  SKI*  T*  THI*  TOP* 

S  V.  XL.  WLO*  XP«  XSG*  XSGL*  YP*  YSG*  YSGL 
COMMON  HPREV«SKPRtV*PKPREV»U80T,UMI0»ET*CGP.P0W*PL»DS 
GO  TO  <  140.100) *  IGO 
100  CALL  DEPTH  (X.Y) 

IF  (DXY*200.0-*L)  110*140*140 
110  IF  (DXY )  120.120.130 
120  RETURN 
130  JG0  =  2 

APG=32. 1725* (DXY) 

CXY=SORT ( ARG) 


DCDH=16.0fi623/CXY 
GO  TO  170 
140  CI=CXY 
JGO  =  l 

DO  150  1=1*50 
ARG=< <DXY)*SIG)/C1 
ARG=TANH(ARG) 

CXY=CO*ARG 

IF  (AHS(CXY-CI)-O.OOOl)  160.150.150 
150  CI=(CXY*CI)*0.5 
160  RCCO=CXY/CO 

SCMC= ( 1 ,0-RCCO*RCCO) *SIG 
V=SCMC* (DXY) ♦RCCO*CXY 
DCUH=CXY*SCMC/V 

170  PHX=E(4)*2.0*XP*E<b)*YP*E(2) 

PHY=E  <6>*2.0*YP*F  (5)*XP»E  (3) 

FK= (STN ( A ) *PHX-COS ( A) «PHY) »DCDH*OCON/CXY 

RETURN 

END 


I 


m 1 ' S f  tQjvr-*  1 -  ***** 


I- 


l! 

1: 


I 
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C  DEPTH  SUBROUTINE  DECK 

SUBROUTINE  DEPTH  (X.Y) 

COMMON  OEP ( 20000 ) ,MI .MJ.SCX.SCALE .GB ID.PD 

COMMON  HI.  R2.  BOUND.  CKDEP.  CO.  CXY.  DUE).  OCOH.  DCON. 

1DELTAS.  UPC.  DTGR.  OXY.  E  (6)  •  •GPINC. 

2  H,  HH.  HO.  ICN.  IGO,  IHET.  ISP.  JGO.  KFIRST,  KPLOT. 

3  LIMNPT.  LPLOT,  NPR I  NT .  NPT.  PHX.  PHY.  R90.  RADIUS.  PAYNO. 

4RCCO.  PHS.  PK.  SCNV.  SIG.  SK.  SKI.  T.  THI.  TOP. 

5  V.  wL.  WLO.  XP,  XSG.  XSG1. .  YP.  YSG.  YSGL 
COMMON  HPPEV.SKPPEV.PKPPEV.UHOT.UM If), FT. CGR.POW.PL. OS 
0IMFNS10N  S  < 12. M 
S ( 1* 1) =0.30861241 
S <2, 1 >=0.23664207 
S (3. 1 ) =0.21770331 
S (4, 1 ) =0.236b4207 
S<5.1)=-0. 08492823 
S  <6. 1 )  =-0.05143541 
S  <7. 1) =-0.05143541 
S (6. 1 ) =-0.08492823 
S (9. 1 ) =0 . 005y8086 
S<10. 11=0.13038277 
S(ll. 11=0.13038277 
S(12, 11=0. 00598086 
S ( 1 .2) =0.05322964 
S<2. 21=0. 19677030 
S(3» 21=0. 14413872 
5(4,21=0.10586122 
S<5. 21=0. 09031100 
S (6,21 =-0.0675*374 
S ( 7.2) =-0.03349283 
S ( 8 . 2 1 =0.03349282 
S(9,21=-0. 18241626 
S ( 10,21 =-0.34031099 
S(11.2>=-0. 12440190 
S(12, 21=0. 12440190 
S ( 1 ,31 =0.05322964 
S (2.31 =0. 10586122 
S<3. 31=0. 14413872 
S (4,3) =0. 19677030 
S  <5.31 =0.03349282 
S<6. 3) =-0.03349282 
S (7.3) =-0.06758374 
S(8» 31=0. 09031099 
S (9.3) =0.12440190 
S(10,3)=-0. 12440191 
S(11.3)=-0. 34031099 
S ( 12.3) =-0. 18241625 
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5(1. 4>=-0. 12499998 
5 (2.4) =-0. 12499998 
S<3.4)=-0. 124^9998  J 

S  (4.4) =-0.12499998 
S<5.4)=0.125 
S<6.4)=0.125 
S ( 7 • 4 ) =0. 0 
5  (8.4) =0.0 
S<9.4)=0. 12499999 
5 ( 10*4) =0. 1249999 
S ( 1 1 « 4 ) =-G  .  0 
S  ( 12*4) =-0.0 
S  < 1 .5) =0.05263157 
S (2.5) =-0.05263157 
5  (3.5) =0.05283158 
S (4.5) =-0.05263157 
S  <5.5)=-0. 15789473 
S<6, 51=0.15789474 
S(7«5>=0. 15789474 
S (8.5) =-0. 15769473 
S (9*5) =-0.15789473 
5(10. 5)=0. 1578^473 
S  C 1 1  *S) =0 . 15789473 
S < 12.5) =-0. 15789473 
S<1.6)=-0. 12499998 
S (2.6) =-0.12499998 
5(3. 6)=-0. 12499998 
5 (4.6) =-0. 12499998 
5 (5.6) =0.0 
S (6.6) =0.0 
S(7«6)=0.1 25 
5 (8.6) =0. 125 
S (9.6) =-0.0 
S ( 10.6) =-0.0 
S(ll.6)=0. 12499999 
5 (12. 6)=0. 12499999 
I =X ♦ 1  • 

J  =  Y  ♦  1  • 

XP= AMOD ( X . 1 • ) 

YP=AMOD(Y.l.) 

IF(NPT-l)  120.120.100 
100  IF ( IP— 1 )  120.110.120 
110  IF(JP-J)  120.150.120 
120  IP=I 
JP=J 
I=JP 
J=IP 

I  I  AX  =  (I  -  1)  *  MI 


H. 

ri  '• 


D  ( 3) 

S 

DEP 

(II  AX  ♦ 

MI  4  J  ♦ 

1) 

n<4) 

= 

OEP 

(IIAX  ♦ 

MI  4  J) 

0(b) 

= 

OFP 

( I I  AX  4 

J  ♦  2) 

0(6) 

= 

OEP 

( I I  AX  4 

MI  4  J  ♦ 

2) 

0(7) 

= 

OEP 

( 1 1  AX  4 

2  *  MI  4 

J  ♦ 

* 

0(8) 

s 

OEP 

( 1 1  AX  ♦ 

2  *  MI  4 

J) 

0(9) 

= 

OEP 

( 1 1  AX  4 

MI  ♦  J  - 

1  ) 

0(10) 

=  OEP 

(  I  I  AX 

♦  J  -  1) 

0(11) 

=  OEP 

<  1 1  AX 

-  MI  ♦  J) 

0(12) 

=  OEP 

(  1 1  AX 

-  MI  ♦  J 

♦  1) 

1) 


DO  140  K  =  l«6 
SUM=0. 

no  130  L  =  1  4 12 

SUM=SUM4D(L)*S(L«K) 

130  CONTINUE 
F (K) =SUH 
140  CONTINUE 

lbO  nXY=(F(l)*E(2)*XP*E(3)*VP*£ (4)*XF  #XP4 
1  E (b) *XP*YP4£ (6) *YP*YP> *DCON 
RETURN 
END 


c  refraction  coeff  subroutine  oeck 

SUBROUTINE  HEIGHT(A) 

COMMON  UEP (20000) .MI .MJ.SCX.SCALE.GRID.PD 
COMMON  B1  *  02*  BOUNO*  CKDEP *  CO*  CXY.  U  (12)  *  DCf)H«  DCON, 
1DELTAS.  ORC*  OTGR «  DXY*  E<6>.  GRINC* 

2  H*  HH,  HO.  ICN,  IGO.  TRET,  ISP*  JGO.  KFIRST •  KPLOT * 

3  LIMNPT,  LPLOT,  NPRINT *  NPT ♦  PHX*  PHY*  RRO.  RADIUS.  RAYNO. 
4QCCO.  RHS .  RK.  SCNV,  SIG.  SK*  SKI.  T*  THI,  TOP. 

S  V,  WL.  WLO*  XP.  XSG,  XSGL.  YP.  YSG.  YSGL 
COMMON  HPREV.SKPRFY.RKPREV.UBOT.UMIO*ET*CGR*POV*.PL*DS 
SKPREV=SK 
RKPREV=RK 
WL=WLO*RCCO 

GN= 12. 5663706 14*DXY/WL 
HS1=EXP(GN) 

HS?  = 1 . /HS 1 

CG=(1.-»GN/(HS1-HS2)*2.)*CXY 
IF  (CG)  100.110*  110 
100  PE  TURN 
110  CONTINUE 

RK  =A6S ( 1 • 0/82) 

SK=SGRT (CO/CG) 

RK  =SQRT ( RK ) 

H=(H0*SK*RK)/SK1 
GO  TO  (120.130) .  JGO 
120  U=-2.0*SIG*RCC0*CXY/ (V*V) 

GO  TO  140 
130  U=— 0 . 5/DX Y 
140  U=U*OCON 

DC»)w=OCDH*OCON 
COSA=COS (A) 

STNA=SIN(A) 

P=-  (COSA‘>PHX*SINA*PHY)  *DCDH*DTGR*2.0 
Q=<  <E  (4) *2 . 0*U*PHX#PHX ) •SINA*SINA-(f (5)« 

1  U*PHX*PHY)*2.0*SINA*COSA* (E (6)*2.0*U* 

2  PHY*PHY) *COSA*COSA) *OCDH*CXY*OTGR*OTGK* 

3  2.0 

63= ( (P-2.0 )*B1* (4.0-0 )*02)/ (P.2.0) 

61=02 

B2=B3 

CALL  FRICTN 

RETURN 

END 
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SUBROUTINE.  TO  PLOT  WAVE.  RATS 

SUBROUTINE  PEF 0 1  A ( XR . YR , I  PEN) 

COMMON  OEP ( 20000 ) »M1 *MJ.SCX.SCALE*GRIU.PO 
COMMON  HI*  fi*?«  ROUND*  CKDEP*  CO*  CXY*  U<1?>*  OCOH,  DCON, 
IDELTAS.  O&C*  PTGR*  PXY*  E(fo>*  0R1NC* 

?  B •  HH,  HO*  ICN*  IGO.  IRFT*  TSP*  JGO*  KFIPST.  kPLOT. 

3  L1MNPT*  LPLOT.  NPPINT*  NPT.  PMX',  PHY*  R90 •  RADIUS*  RAYNO* 
APCCO*  PHS.  0 K*  SCNV.  SIG*  S*  *  SKI*  T«  THl*  TOP* 

5  V,  WL »  WLO*  XP«  XSG*  XSGL*  YP.  YSG.  YSGL 
COMMON  HPPev*SKPREV»RKPRfcV.USOT.UMJD.ET*CGR*POW.PL*OS 
RX=XR*SCX 
RY=YR*SCX-ROUNO 
IF  (PY.LT.O.)  GO  TO  100 
CALL  PLOT ( WX «  WY*  IPfN) 

100  RETURN 
END 
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SUBROUTINE  FRICTN 

COMMON  DfP( 20000 ) ,MI .MJ.SCX .SCALE .GRID. PD 
COMMON  61.  82.  BOUND,  CKDEP,  CO,  CXY,  0(12).  OCDH.  OCON, 
10FLTAS,  DRC,  PTGR,  OXY.  E(b),  GRINC, 

2H,  HH,  HO.  ICN,  IGO,  IRET.  ISP.  JGO,  KFIRST,  KPLOT , 
3LIMNPT.  LPLOT.  NPRINT,  NPT,  PHX,  PHY,  H90.  RADIUS.  RAYNO. 
4RCC0.  RHS.  RK.  SCNV,  SIG.  SK,  SKI,  T,  THI,  TOP, 

SV.  WL.  wLO,  XP.  XSG.  XSGL.  YP.  YSG.  YSGL 
COMMON  HPREV.SKPREV.RKPweV.UBOT.UMIU.E  r.CGR.POW.PL.DS 
OX=DS*GR ID 

Oh=HPREV* (SK/SKPREV) *  <RK/PKP»EV) 

PI=3.141SV26S 

G=32.17 

ARG=2.*PI»DXY/WL 

PHl  =  .6391M ( SK/S I NH ( ARG ) )**3) 

HNE  W=OH/  (  (  •  02*0H*PH  1*1)X)  /  ( SK*T**4 )  *1  .  ) 

HPRFV=H,NEw 

H=HNEW 

THE  ABOVE  CARO  SHOULD  BE  USED  IF  THE  FRICTIONAL  AFFECT 
ON  WAVE  HEIGHT  IS  DESIRED 

UP0T=H*G/(2.*CXY«C0SH(AhG) ) 

UVID=U6U1*C0SH(ARG*.5) 

RETURN 

END 


♦ 
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SUBROUTINE  ENFR6Y (A. END) 

COMMON  OEP<20000> .MI .MJ.SCX. SCALE. GRID. PO 

COMMON  HI.  R2.  ROUND*  C.K.OEP.  CO*  CXY*  0(12).  DCDH.  OCON, 

1  DF.LT AS .  DHC.  OTGR.  DXY.  E(fc>),  G»INC. 

2H.  HH,  HO.  1CN.  IGO.  1RFT.  ISP*  JGO.  KFIRST.  KPLOT  * 

3LIMMPT,  LPLGT  *  NPRINT.  NPT «  PHX.  PHY.  HQO «  RADIUS.  RAYNO. 

4RCC0.  RHS.  RK.  SCNV.  SlG,  SK  .  SKI ,  T.  ThI.  TOP. 

5V »  *L.  WLO*  AP.  XSG.  XSGL.  YP,  YSG.  YSGL 
COMMON  HPREV  .SKPRF. V .RKPRFV  .IJHOT  »UM  10 *E  T .CGR. POw  .PL  *0S 
COMMON  10  I R 

DATA  Lf FT/'L*/.IRIGHT/«W*/ 

PI=3.1415S26S 

IF  (PHY  .FG.  0.)  GO  TO  20 

A  A  =  A 

IF  (AA  .GT.  8.2831RS2)  AA=A-to.2R31fc<S2 
THE TA= AT AM (PHX /PHY) .PI//. 

THET  A=PI-THF  T  A 

I  dtheta=thfta-aa 

IF  CDTHfc'T A  .LT.  0.)  TDIP=IPIGHT 
IF  (PTMFTA  .GT.  0.)  IOlR=LtFT 
OX  =DS*GR I D 
HT=8.*H*H*wL 
ATHFTA=AHS (OTHETA) 

OTHF TA=ATHFTA 
PLP=4.*PI*0XY/WL 
SDLR=SINh (DLR) 

IF  (SDLR  .LT.  l.E-10)  SPLR*1.E-10 
CGR=.S*CXY* ( 1 . ♦OLR/SOLW) 

P0W=8.*H*H*CGH 

PL=POW*SIN  (  A1  MET  A)  *C05  ( ATHF. T A ) 

IF  (END  .NE.  1.)  RETURN 
HP=H 

ALPHA=OTHETA 

PLS  IS  AN  ESTIMATE  OF  ThF  LONGSHORE  COMPONENT  OF  ENEPGY  FLUX  fc* 
SURF  ZONE.  BECAUSE  IN  ThE  CALCULATION  OF  PLS  THE  SIGNIFICANT  Vt 
HEIGHT  IS  USER  INSTEAD  OF  THE  ROOT-MF AN-SuUARE  THE  VALUE  OF  PL? 
APPROXIMATELY  TWICE  THE  VALUE  OF  THE  EXACT  ENERGY  FLUX.  THEREFr 
REAL  YUMY ( 20000 ) 

PLS  IS  MORE  APPROPRI ATtLY  REFERRED  TO  AS  THF  LONGSHORE  ENERGY 
FLUX  FACTOR. 


PLS=32.1*(HH**2.S) *S IN (2* ALPHA) 


PLS=32. I* (HH**2.S) *S IN (£« ALPHA) 
HHM=PH/3.2H1 

PI SM= (2.60*10**10) * (HHM*#2.S)*SIN(2#ALPHA) 
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WRITE  <b,?20)  PLS 

220  FORMAT  CO'.'  LONGSHORE.  ENERGY  FLUX  FACTOR:  PLS='.F14,7. 

*  »  FT-LRS/S2C/FT  OF  HE ACM  FRONT  * ) 

I 

APPROXIMATION  OF  L0N(j5H0WE  TRANSPORT  RATE:  (1)  HASED  on  PLS 
(2)  HASED  ON  GAL V IN < 1 972) -LST I MATE  OF  THE  GROSS  LONGSHORE  TRAN* 

pate. 

OI=7SOO*PLS 
02=20000 0*HH*HH 
QM=.000128*PLSM 
WRITE  (6.230)  Q1.10IR.Q2 

230  FORMAT  <•  LONGSHORE  TRANSPOPT  RATES:*./. •  *.11X» 

1  'FUNCTION  OF  PLS :  0=*.E1*.7,*  CURIC  YARDS/YP * . 5X . *F AC ING  OCEAN: 

2  TRANSPORT  TO  *«Al./«*  «.12X. 

3  'FUNCTION  OF  HR:  G=*«El4.7.'  CURIC  Y ARUS/YR • . SX» • GAL V  IN < 1 R72)  - 
4GP0SS  TRANSPORT  RATE') 

WRITE  (6.2S0)  PLSM,  OM 

2S0  FORMAT  CO'.'  METRIC  CONVEPSI UNS : • ./ . •  *.11X,*  PLS='.El4,7. 

*  '  EPGS/METER*SEC'./.  •  *.13X.»  0='.El4.7.*  CURIC  METERS/YE AQ  • ) 
RETURN 

C 

20  THE  T  A  =  P I 
GO  TO  10 
END 
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SUBROUTINE  BOTTOM 


THIS  PROGRAM  CREATES  A  MAP  WITh  DEPTHS  AT  EQUALLY  SPACED 
DISTANCES  FROM  A  DIBIT  I  Z  If  I.)  H ATHYMtTR I C  MAP  US  I  MB  THE 
CONTOURS 

THIS  PART  OF  THE  PROGRAM  heads  THt  TAPE  AND  CONVERTS  The 
MENU  TO  CORRECT  FORK  FOR  PROCESSING. 

COMMON  DEPTH(20Pn(>>  ,  I X  AX  ,  I  Y  AX  ,  SCX  »  SC  ALt  ♦  GR 1 0  ♦  PD 
COMMON  R 1 .  HP.  ROUND 
DIMENSION  A  <  <*0R6) 

DIMENSION  IPA<5),  IPY(S) 

DIMENSION  IOEPTh(IOO) 

DIMENSION  PPU(IOOO) 

DIMENSION  IA IPOOOO) . TY (POOOO) 

DATA  AO»A]»AP*A3.A4.AS.An,A7»Ab*A9/*0*»*l*.*2*»*3***4***5*,«f,l, 

I  *7***8*»*9»/ 

DATA  ah. an/ *h * t *N * / 

DATA  AST/***/ 

DATA  AL  / *L  •/ 

OAT  A  ASLSH  «  API.  US  ♦  AM  I NUS/  •  /  •  ,*♦*,*-*/ 

data  ad*  aej  «  af  *  ap.  as/  *d  *  «  *e  *«*f*.*p***s»/ 

IPPS4=0 
IPEN=0 
ILL  =  0 
IXSW  =  0 

IP  =  0 
IX ( 1 )  =  0.0 
I Y ( 1 )  =  0.0 
on iv  =  io.o 
I RSW  =  0 
ISHR  =  0 
I END  =  0 
IIISw  =  0 
ICT  =  0 
ISTOP  =  1 
8000  FORMAT  (3A4) 

IS  =  0 
I X  Y  =  1 

II  =  1 
ID  =  0 
IRT  =  -1 
IN  a  B 

I OUT  =  b 
1 T  APE  =  17 
I  I D I SK  =  IP 
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ISLSH  =  0 

inumh  =  0 
ISIGN  =  0 
1  END  =  0 
III  =  0 

PFAD  INPUT  TAPE  F»OM  ulGITIZFP 

200  »FAO  ( I T  APE  1 9000 )  ( A ( I ) *  1  =  1 *4096) 

9000  FOKMAT  (64 ( 64 A  1 )  ) 

00  2000  I  =  t ♦ 4096 
IK  (A(I).EO.AO)  GO  TO  bOO 
IF  ( A ( I) .FO. Al )  GO  TO  S20 
IF  ( A ( I ) .E0.A2)  GO  TO  540 
IF  ( A  ( I )  . F. 0 •  A 3 )  GO  TO  SbO 
IF  (A ( I ) .EQ.A4)  GO  TO  580 
IF  ( A ( J ) «F 9. A5)  GO  TO  600 
IF  (A ( I ) .FO.A6)  GO  TO  620 
IF  (  A  (  I  )  # F  U • A 7 )  GO  TO  640 
IF  (A(I).FQ.AB)  GO  TO  660 
IF  ( A ( I ) •  EQ  •  A9)  GO  TO  680 
IF  (All) .FO.ASLSH)  GO  TO  940 
IF  ( A ( I ) .FU.APLUS)  GO  TO  960 
IF  (A (I ) .EQ.AMINUS)  GO  TO  1010 
IF  (  A  ( I )  .E'O.AP)  GO  TO  1040 

IF  (A(I).FO.AF)  GO  TO  1060 

IF  (A(I).tQ.AF)  GO  TO  4500 

IF  ( A ( I) .EO.AS)  GO  TO  1080 

IF  ( A ( I ) .EO.AP)  GO  TO  1140 

IF  (A(I).EO.AL)  GO  TO  1135 

IF  (A(I).EO.AH)  GO  TO  4300 

IF  (A(I).EQ.AN)  GO  TO  4200 

IF  (All)  .EO. AST)  GO  TO  4230 
GO  TO  2000 

500  INU^B  =  INUMd  *  10 
GO  TO  690 

520  I NUMB  =  1NUM6  *  10  ♦  1 
GO  TO  690 

540  INUMH  =  INUMd  *  10  ♦  2 
GO  TO  690 

560  INU6B  =  INUMb  *  10  ♦  3 
GO  TO  690 

56 0  I  NUMB  =  INUMd  *10  +  4 
GO  TO  690 

600  I NUMB  =  INUMH  *  10  ♦  5 
GO  TO  690 

620  I NUMB  =  INUMd  *  10  ♦  6 
GO  TO  690 

640  INUMb  =  INUMd  *  10  ♦  7 
GO  TO  690 


11-82 


660  IMtMB  =  JNUMct  *  10  ♦  8 
GO  TO  690 

660  INUKH  =  INIJM8  *  1  ft  ♦  9 
690  ICT  =  ICT  ♦  1 

IF  <IPPSw  .EG.  1)  GO  TO  2000 
IF  (IS.FO.l)  GO  TO  2000 
IF  (ICT.GT.5)  GO  TO  1210 
IF  (ICT.LT.5)  GO  TO  2000 
GO  TO  ( 730  »  790). IXY 
730  I X  Y  =  2 

XSS  =  FLOAT  ( I NUMB )  *  SCALE  /  1000.0 
IF  USIGN.NF.2)  GO  TO  660 
XSS  =  -  XSS 
GO  TO  861) 

790  I X  Y  =  1 

YSS  =  FLOAT  (ItMUMM  *  S  C  A  L  k  /  1000.0 
IF  (  IMGN.NE  •?)  GO  TO  H30 
YSS  =  -  YSS 

830  IX  (II)  =  XSS  ♦  0.5 
IY  (II)  =  YSS  ♦  O.S 
IF  ( 1 ST  OP  .F  0 •  1)  GO  TO  860 
II  =  II  ♦  1 

IF  ( I  I .Gl  . 20000)  GO  TO  1310 

850  CONTINUE 

IF  (II  .LE.  2)  GO  TO  851 

ITHX  =  JAMS  (IX  (II  -  1)  -  IX  (II  -?>> 

XHl l=FLOAT (IIHX) 

IF  (XMII  .LT.  APS  <  GW  T  0* . 5) )  GO  TO  651 
IX  (II)  =  IX  ( I  I - 1 ) 

IY  (II)  =  TY  < 1 1-1  ) 

IX  ( I  I - 1 )  =  (IX  (II)  ♦  IX  ( II-2) )  /  2 
IY  (IT-1)  =  (IY  (II)  ♦  IY  (11-2) )  /  2 
II  =  II  ♦  1 

851  CONTINUE 

IF  (IP  .LE.  5)  GO  TO  860 
OP TH= IOEPTH ( ILL) /10.0 
wPITE  (11,9980)  XSS.YSS.OPTH, IPEN 
9980  FORMAT  (4A4) 

IF  (1SHP.EQ.1)  GO  TO  4400 
IPEN=? 

860  I S I GN  =  0 
ICT  =  0 
I NU^H  =  0 
GO  TO  2000 

940  IF  (ISLSH.EQ.l)  GO  TO  1330 
IE  (ICT  .LT.  5)  GO  TO  1220 
ISLSH  =  1 
ILL  =  ILL  ♦  1 

IF  (ILL. GT. 100)  GO  TO  1360 

IF  (  1  5  I  CiN  .  EO.  2)  I  Nt  (MB  =  -  I  NUMB 
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I!htJ7H  (ILL)  =  INUW‘- 
Ic>  =  0 
10  =  0 
00  TO  HOO 

IF  (ISIGN.NE.O)  GO  TO  13S0 
IS i GO  =  1 
GO  TO  ?0  0(> 

IF  (  I  S  I  G I :  .  N  F  .0)  On  TO  1  ISO 
IS  I  GO  =  ? 

GO  TO  ?0  0  0 

IF  t  Ili.F  U.  1  )  GO  TO  1370 

I  GI..SO  =  0 
IS  =  1 
I^F.N=3 
I!'  =  1 
IS  10?  =  0 
GO  TO  hG  0 

I  STOP  =  1 
IShF>  =  0 

IF  (1P.OI.S)  GO  TO  4100 
IF  (  IPSw.FQ.O)  GO  TO  10HO 
IF  (lPf.GF.l)  GO  TO  1170 
IXY  =  1 

II  =  II  ♦  1 

ISLSH  =  0 

I  PSw  =  0 

GO  TO  H  GO 

IlL  =  ILL  +  l 

I  OF  PTH  (ILL)  =  I  HOP T  H 

IX  (IT)  =  0.0 

I Y  (IT)  =  0.0 

[PSW  =  1 

jPT  =  IPT  ♦  I 

IP  =  IP  ♦  1 

IPX ( TP)  =  IX ( I  I  ) 

I P  Y  ( IP)  =  I  Y  (  I  1  ) 

GO  TO  «60 

IF  ( IPT .PT . 1 )  GO  TO  3010 


GO  70  3010 


X  -['■  XI  HOT  TOM 


I  WOP  TH  =  IOE  hTh  U) 

NON  =  0 

IXOIS  =  IX  (ID  -  IX  (1) 
T  Y  •  *  f  S  =  IY(II)  -  I Y  < 1 ) 
XI  >  I  S  =  Ft  oAT  (  I  XT'  1  S  > 
fl-IS  -  FI  ;,AT  (  1  YOIS) 

/  : )  I  S  r  AI'G(XOlS) 

Y[.)c.  ;  AHS(YIjTS) 

I  I-  (  .•  1 1  !  N  .  >  ')  .  li  .  0  )  xl'i  I  S 


=  0.(1001 


•msssst. 


.><*x^hiiiiaitoGw. 


ono  nno  noo 


t  m  t  t  a  i  =  atan< ydis/xdis) 

DIb  =  so«t  (X0IS**2  ♦  ynis**?) 
orsi  =  ois 

ALEET  =  0.0 

CALCULATE  LENGTH  OF  INCPFMENTS  ALONG  X  -  BOTTOM 


AXO  =  OIS  /  FLOAT  ( I X AX  -  1) 

1  I  =  T I  -  1 

00  3000  LM  =  1*11 

IXX  =  TAHSdX(LM)  -  I  X  <  LM*  1  )  ) 

IYY  =  IAH5(IY(LM)  -  lY(LNd>> 

XX  =  IXX 
.  YY  =  IYY 

CALCULATE  DISTANCE:  TO  NEXT  POINT 

DO  =  SUHT  (XX  **  2  ♦  YY  »*2> 

IOC=  I  DEPTH  (LM*  1 )  -  IOFVTrt(LM) 

ADC  =  FLOAT  (IDO  /  DO  IV 

0PH  =  FLOAT  (IDFPTH  (Lm))  /  ODIV 

DPH2  =  ELOAT  (IDFPTH  <LM*1)>  /  DOIV 

IF  (ALFFT  .G I .  00)  GO  TO  2990 

IF  (ALEFT.LE.0.0)  GO  TO  ?“00 

NNN  =  NNN  ♦  1 

DEPTH  (NNN)  =  OPH  ♦  ALFFT  *  ADC  /  DO 
OPH  =  DEHTH  (NNN) 

00  =  00  -  ALcFT 
2900  AN  =  PO/AXO 

AIDC  s  (0PH2  -  QPh )  /  AN 
NN  =  AN 

IF  (NN.FO.O)  GO  TO  ?9b0 
CALCOLATE  DEPTH  INCREMENTS 


DO  2950  LN  =  1»NN 
NNN  =  NNN  ♦  I 

DEPTH  (NNN)  =  OPH  ♦  FLOAT  (LN)  *  AIDC 
2950  CONTINUE 
2960  CONTINUE 

ALEE  1  =  DD  -  ELOAT (NN)  *  AXD 
ALEET  =  AXO  -  ALFFT 
GO  TO  3000 

2990  ALEET  =  ALFFT  -  DD 
3000  CONTINUE 

II  =  II  ♦  1 

IF  dY  (1  )  .LT.  IY  (I  I)  )  GO  TO  3003 
I  Y  (  1  )  =  1YIII) 

GO  TO  3004 
3003  IY(II)  =  I Y ( 1 ) 
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jo 04  continue: 

IP  =  I X  AX  -  I 
00  3005  J  =  1  #  I H 
ppo  <jm>  =  ufpth<j) 

3005  CONTINUE 

MMM  =  IH  ♦  1 
IX  (  1  )  =  IX ( I  I ) 

I  V  (  1  )  =  IV(II) 

II  =  1 

I  DEPTH  (1)  =  IOEPTH  (ILL) 

ILL  =  1 

GO  TO  10R0 

3010  IF  ( I p T . GT  •  2 )  GO  TO  3410 
ALEFT  =  0.0 

Y-AXIS  -  FAR  S I OF 

NNN  =  0 

IXUIS  =  IX  (  I  I  )  -  1 X  (  1  ) 

IYDIS  =  IY  (  I  1  )  -  IY  (  1  ) 

XOIS  =  FLOAT  (IXCIS) 

YD  I S  =  FLOAT  (IYOIS) 

XOIS  =  AMS(XOIS) 

YDIS  =  ABS(YDIS) 

IF  (XOIS.FO.O.O)  XDIS  =  0.0001 
THETA?  =  AT  AN  (YD1S/XDIS) 

ms  =  sort  ( x r > i s  **  2  ♦  ydis  **?> 

OIS2  =  01 s 

CALCULATE  lfngth  OF  INCREMENTS  along  y  -  right  side 

AYD  =  DIS  /  FLOAT  (IYAX  -  1) 

II  =  II  -  1 

00  33N0  LM  =  1*11 
I  XX  =  IAHS  (IX(LM)  -  1 X ( L  N ♦ 1 )  ) 

IYY  =  I  APS  (  1 Y  ( L  m )  -  IV (LK* 1)1 
XX  =  IXX 
YY  =  IYY 

DD  =  SUPT  (XX  **  2  ♦  YY  2) 

IOC  =  I  DEPTH  (LM  ♦  1)  -  Il)F.RTh(LM) 

ADC  =  FLOAT  (IOC)  /  00 1 V 

OPH  =  FLOAT  (IOEPTH  (LM))  /  DUIV 

DPh2  =  FLOAT  'IDFPTh  (LM* 1 ) )  /  ODIV 

IF  (ALFFT  . G i .  00)  GO  TO  3350 

IF  (ALFFT. LF. O.o)  GO  TO  3270 

NNN  =  NNN  ♦  1 

DEPTH  (NNN)  =  OPH  ♦  ALEFT  *  ADC  /  DO 
DPH  =  DEPTH  (NNfJ) 

OH  =  DO  -  ALEFT 
ism  an  =  no/  ayo 
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a  I  DC  =  <C>ph2  -  U PH)  /  an 

NN  =  AN 

IF  (NU.H1.0)  GO  TO  333<» 

*  00  3320  LN  =  l.NN 

NNN  =  NNN  ♦  1 

DEPTH  (NNN)  =  DPP  ♦  FLOAT  (LN)  *  Al()C 
3320  CONTINUE 
3330  CONTINUE 

ALEFT  =  DO  -  FLOAT (NN)  *  A YD 
ALEFT  =  A YD  -  ALFFT 
GO  TO  3360 

33S0  ALEFT  =  ALFFT  -  DO 
3360  CONTINUE 

I  1  =  1 1  ♦  I 

IF  (  IX ( 1 ) .GT. IX ( T I ) )  GO  TO  336b 
IX ( 1 )  =  IX  (  I  I ) 

GO  TO  3367 
336b  IX(II)  =  IX<1) 

3367  CONTINUE 

MHMM  =  MHM  ♦  1 Y  AX  -  2 

00  3370  J  = 

JJJ  =  J  -  ♦  1 

PP!)  <J*1)  =  DEPTH  (JJJ) 

3370  CONTINUE 

MMP  =  MHM.H  ♦  1 

IX ( 1  )  =  I  X  (  I  1  ) 

1 Y ( 1 )  =  lY(Ii) 

II  =  1 

IDE.PTH  (1)  =  IOEPTH  (ILL) 

ILL  =  1 
GO  TO  10S0 

3410  IF  (IPT.GT.3)  GO  TO  37«0 
X-AXIS  -  HACK 
NNN  =  0 

IXDIS  =  IAHS  (  I  X ( I  1 )  -  t  X  (  1  ) ) 

IYOIS  =  IAHS  (IY(II)  -  I  Y  (  1  ) ) 

ALEFT  =  0.0 
XDIS  =  FLOAT  (IX0I5) 

YO IS  =  FLOAT  (  IYOIS) 

IF  (XOIS.FU.U.O)  XOIS  s  0.0001 
THETA3  =  A TAN  (YCIS/XUIS) 
ms  =  SOPT  (XOIS  **  ?  ♦  YU  I S  **  2) 
nis3  =  ois 

CALCULATE  LENGTH  OF  INCPF  HE. NTS  ALONG  X 

AXD  =  DIS  /  FLOAT  (  J  X  A  X  -  1) 

II  =  II  -  1 


-  TOP 
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00  3730  LM  =  1  •  1 1 

IXX  =  IflHS  (lX(L^)  -  IX(LKM>> 

IVY  =  IAMS  (1Y(LV)  -  I Y (LM* 1 ) ) 

XX  =  IXX 
YY  =  IYY 

no  =  SOOT  (XX  **  2  *  YY  **  2) 

IOC  =  I0KPTH(LM*1)  -  II)h>TH(LM) 

ADC  =  FLOAT  (IOC)  /  Dl)I V 

DPh  =  FLOAT  (IOEPTH  (Lm) )  /  UDIV 

DPH?  =  FLOAT  ( IOFPTrt  < L M ♦  1 1 )  /  ODIV 

IF  (ALfcFT  .GT.  DO)  GO  TO  3/20 

IF  (ALEFT.LF.C.O)  GO  TO  JGbO 

NNN  =  NNN  ♦  1 

OFPTH  ( NNN )  =  OPH  ♦  ALfcFT  *  ADC  /  DO 
HPH  =  DEPTH  (NNN) 

00  =  00  -  ALfcFT 
36bO  AN  =  00  /  AXO 

A I  DC  =  ( OPH 2  -  OPH)  /  AN 

NN  =  AN 

IF  (NN.EO.O)  GO  TO  3710 
00  3700  LN  =  I « NN 
NMN  =  NNN  ♦  1 

DEPTH  (NNN)  =  OPH  ♦  FLOAT  (LN)  *  AIOC 
3700  CONTINUE 
3710  CONTINUE 

ALfcFT  =  00  -  FLOAT (NN)  *  AXO 
ALfcFT  =  AXO  -  ALFFT 
GO  TO  3730 

3720  ALFFT  =  ALFFT  -  OD 
3730  CONTINUE 

II  =  II  *  1 

IF  ( I Y ( 1 ) .GT.IY(TI) )  GO  TO  373b 
I  Y  (  1 )  =  1  Y  (  I  1  ) 

GO  TO  3737 
373b  IY(II)  =  I Y  <  1 ) 

3737  CONTINUE 

MM  MM  =  MN,M  ♦  IX  AX  -  2 
00  3740  J  =  viMM.mmmm 
JJJ  =  J  -  MMsA  ♦  1 
PPD  <J*1)  =  UfcPTH  (JJJ) 

3740  CONTINUE 

MWM  =  ♦  1 

IX (1)  =  1X(II) 

I  Y  (  1  )  =  I  Y  <  I  I  ) 

IT  =  1 

I  DEPTH  (1)  =  I  OF  P  T  H  (ILL) 

ILL  =  1 
GO  TO  10S0 

3780  IF  (JPT.NF.4)  GO  TO  4020 
C 
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Y-AXIS  DOWN 


NNN  =  0 

IXDIS  =  I AHS  (IX  (  II)  -  IX  (  1  )  ) 

I  YD  IS  =  I AHS (  I  Y  (  II)  -  IY  (1  )  ) 

xrns  =  float  (ixnis) 

YOIS  =  FLOAT  < I  YD  I S ) 

IF  (X f IS. F O.O.O)  XOIS  =  0.000] 

THETA*  =  ATAN  (YDIS/XDIS) 

ms  =  SOFT  (XOIS  **  2  ♦  Yt 1 1 S  **  ?) 

0IS4  =  f  •  1 S 

ayd  =  ois  /  float  uyax  -  n 

11=71-1 

ALEFT  =  0.0 

00  3990  LM  =  1,11 

I  XX  =  IAHS  (IX(LV)  -  I  X ( LM* 1 ) ) 

I Y  Y  =  JAMS  UYtLFO  -  IY(LM*D) 

XX  =  IXX 
YY  =  IYY 

00  =  SORT  (XX  **  2  ♦  YY  **  ?) 

IOC  =  I Ok PTH (LM* 1 )  -  I DEPTH (LM) 

ADC  =  FLOaT  (IOC)  /  DO IV 

OPH  =  FLOAT  (IOEPTH  (L«) )  /  001V 

()Ph2  =  FLOAT  (IL'FPTh  (LM*1>)  /  00 IV 

IF  (ALEFT  ,GT.  DO)  GO  TO  3980 

IF  (ALFF T.LF.0.0)  GO  TO  3910 

NNN  =  NNM  ♦  1 

DEPTH  (NNN )  =  OPH  ♦  ALEFT  *  ADC  /  00 
DPH  =  DEPTH  ( NNN ) 

DO  =  00  -  ALEF  T 
3910  AM  =  00  /  AYD 

A I OC  =  (OHM?  -  DPH)  /  AN 
MM  =  AN 

IF  (NM.fG.O)  GO  TO  3970 
DO  3960  LM  =  1,MN 
MM N  =  NNN  ♦  1 

OfePTH  (IMMN)  =  DPH  ♦  FLOAT  (LN)  *  AlDC 
3960  COMT1MUE 
3970  CONTINUE 

ALEFT  -  OD  -  FLOAT  (NN)  *  AYD 
ALEFT  =  AYD  -  ALFF  T 
GO  10  3990 

39m0  A|  FFT  =  ALEFT  -  00 
3990  CONTINUE 
IP  =  G 

I  X  ( 1 )  =  0 

IX ( II )  =  0 

MPMM  =  MMM  ♦  I YA  X  -  3 
DO  A 000  J  = 

JJJ  -  J  -  MMM  ♦  l 


no  onnnoo 


PPD  ( J ♦ 1 >  =  DF^TH  (JJJ) 

PI  =  3.141693 
4000  CONTINUE 

PPD  (1)  =  OEPTHIPMMM  -  MMM  ♦  ?) 

II  =  0 

THETA1  =  THETA1  *  PI  /  IttC.O 

THETA?  =  THFTA2  •  PI  /  1E0.0 

THETA3  =  THE  T  A3  *  PI  /  1*0.0 

THETA4  =  T  H  E  T  A  4  *  PI  /  1H0.0 

FIND  INCPFMENTS  OM  X-AXIS  -  hOTTOM 

ANC1  =  nisi  /  FLOAT  ( 1 X a X  -  1) 

FINP  IHCPEMEwTS  ON  X-AXIS  -  TOP 

ANC?  =  l»  I  S3  /  FLOAT  (  I X  AX  -  1) 

ANCAV  =  <I)IS1  ♦  OIS3)  /  ?.0 
ancav  =  AnCAV  /  FLOAT  (IXAX  -  1) 

ILL  =  0 
00  TO  1090 
4G?0  WPITE  ( 6i 6030) 

0030  FOPmat  (»1  L  ENTFPEf)  E  AWLY  •  > 

STOP 

4100  II  =  II  -  1 

IF  (II  .LF.  1)  00  TO  1090 
IXD  =  99999999 
DPTh  =  0.0 

I V  H  =  0 
I  X  H  =  0 
IISW  =  0 

MMM  =  0 

XVAL  =  FLOAT  (IX  ( 1 ) )  /  ANCAV 
IXVAL  =  I F I X  (XVAL  ♦  0.6) 

XI.EN  =  FLOAT  (  I  XVAL )  *  ANCAV 
IXH  =  XLF  N 
IYH  =  IY  (1) 

DPTH  =  IDFPTH  (ILL)  /  10.0 

wP I TF  (JQISK,  3000)  IXH,  IYH,  DPTH 

no  4ioo  j  =  ?,n 

4110  JL  =  J  -  1 

IF  (  (  IX  ( J)  ,F'J,  IX  ( JL)  )  ,AN(>.  (  IY  (J)  ,F0.  IY  (JL)  )  )  60  TO  4190 

CALCULATE  THE  MANGE 
XVAL  =  FLOAT  <IX(J))  /  ANCAV 
C 

I XVAL  =  I F I X  (XVAL  ♦  0.6) 

XL EN  =  FLOAT  (IXVAL)  *  ANCAV 
IF  (ITSw.FQ.l)  GO  TO  41S0 
4120  I XH 1 =  I F I X  (XLEN  ♦  0.6  *  ANCAV  ♦  0.6) 
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IXH2  =  [FIX  ( XL t N  -  0.5  *  ANCAV  ♦  0.5) 

IXXAV  =  XL£N 

4150  IF ( ( IX (J> .GT.IXHl) .OP. ( IX ( J) .LT. IXR2) )  GO  TO  4]R0 
I IXD  =  1 AhS  ( I  XX  AV  -  IX (j>) 

IF  (IIXO.GE.1XO)  GO  TO  4190 
IXfj  =  IIXU 

I  XH  =  I  X  X  A V 
I YH  =  I Y ( J ) 

USW  =  1 

OPTH  =  IOFPTH  (ILL)  /  10. 0 
GO  TO  4190 

4 1  HO  IF  (HSW.fcO.O)  GO  TO  4190 
1 1 SW  =  0 

WRITE  (IOTSK.HOOO)  IXH,  1  YM,fjPTH 

I X  0  =  99999999 
GO  TO  4120 
4190  CONTINUE 
ILL  =  0 

II  =  0 

GO  TO  1090 

SET  PEN  SWITCH 

4200  IPPSW=1 

GO  TO  bh 0 
4230  IPFSW=o 

I  PF.N=  I  NUMh 

IF  (  IMIImh.LT.  1)  GO  TO  42H0 
IF  (INUMH.GT.3)  GO  TO  4250 
GO  TO  HHO 

42H0  WRITE  ( 1  OUT  *  9020 )  INIJMS 

9020  FORMAT  Cl  PEN  VALUE  LT  1  OR  GT  3  C*I3*»>*> 

STOP 

4300  ISHP  =  1 

GO  TO  1040 

4400  WRITE  (10*9990)  XSS*YSS.IHtN 
9990  FORMAT  (3A4) 

IPFM=? 

GO  TO  A60 
4500  REWIND  12 
RFWIND  10 
REWIND  17 
CALL  SORT l 

NOw  RUT  PER  I ME  TF  G  OM 

III  =  I X  A  X  *  lYAX 
00  4220  K  =  1 . I  IT 
OFPTH  (K)  =  9999 • 9 

4220  CONT1MUF 


n  o  o  o  r>  n 


no  42SO  K  =  l.IXAX 
DEPTH  (K )  =  PPO  (M 
42S0  CONTINUE 

INI  =  I X AX  ♦  1 
IN  =  I Y AX  -  2 
IP  =  I X AX  ♦  1 

IPl  =2*  IX AX  ♦  2*  IYAX  -  4 
DO  4320  K  =  1  .  I N 
TN2  =  INI  ♦  I X  AX  -  1 
DEPTH  (IM)  =  PPO  (]P1) 

DEPTH  (IN2)  =  PHD  (]P) 

INI  =  INI  ♦  IX  AX 
IP  =  IP  ♦  1 
IPl  =  IPl  -  1 
4320  CONTINUE 

INI  =  (IYAX  -  1)  *  1 X AX  ♦  1 
IP  =  2  *  I X AX  ♦  IYAX  -  2 
00  4300  K  =  l.IXAX 
DEPTH  (INI)  =  PPD  (IP) 

IP  =  TP  -  1 
INI  =  INI  ♦  i 
4300  CONTINUE 

now  STAPT  FILLING  In  THE  POINTS 

IDIS  =  0 
NSE.O  =  0 
IP  =  2 
IPT  =  2 
ISP  =  I  X  A X 
4SS0  CONTINUE 

PF  AD  (13,8000)  IXH,IYM,OPTH 
IF  (IXH  .FQ.  0)  00  TO  4SS0 
4560  CONTINUE 

ALEFT  =  0.0 
ODPTH  =  DEPTH  (IPT) 

AYD1S  =  ( D I S 2  ♦  MS4)  /  2.0 
AYNC  =  AYDIS  /  FLOAT  (IYAX  -  1) 
ISTY  =  0 
I/SW  =  0 
4S80  CONTINUE 

CALCULATE  4  OF  SEGMENTS 

INC  =  TAPSUYH  -  1  STY ) 

00  =  FLOAT  (INC) 

CHGD  =  DPTH  -  ODPTH 
IF  (ALEFT .LE. 0.0)  GO  TO  4bbS 
DD  =  DO  -  ALEFT 
IF  (PD.LT.O.U)  GO  TO  4870 
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NFwPT  =  IPT  ♦  ISP 

DEPTH  (imFwPT)  =  OOPTH  ♦  ALEFT  *  CHGD  /  FLOAT  (INC) 

OOPTH  =  DEPTH  (NFwPT) 

TFT  =  IFT  ♦  IX AX 
4585  CONTINUE 

IF  (DO.FO.O.U)  (30  TO  46TB 
an  =  no  /  aync 

NSEG  =  AN 

A  IOC  =  (FIFTH  -  OOFTH)  /  AN 
IF  (NSEG.LE.O)  00  TO  4601 
00  4600  LMN  s  l.NSEG 
NFwPT  =  1  FT  ♦  ISP 

DEPTH  (NEWPT)  =  OfiPTH  ♦  a  I  DC  *  FLOAT (LMN) 

IPT  =  IFT  ♦  I X AX 

4600  CONTINUE 

4601  CONTINUE 

ALF.FT  =  DO  -  FLOAT (NSEG)*  AYNC 
ALtFT  *  AYNC-  AL^FT 
460ic!  CONTINUE 

OOPTH  =  OPTH 

READ  ( 13.8000 ♦  EN0  =  457  0)  1XA.1YA.DPTH 

IF  ( 1 X A  .GE.  IPX  (2)  .OF.  1 X A  .Of.  IPX  (3))  60  TO  4570 
IF  (IXA.E0.IXH)  GO  TO  4660 

4605  ISTY  =  AYDIS  ♦  0.5 

] PPT  =  I X  AX  *  ( I Y A X  -  1)  ♦  IP 
CHGD  =  DEPTH  (IPPT)  -  OOPTH 
INC  =  IAHS  (IYH  -  ISTY) 

00  =  FLOAT  (INC) 

IF  (ALEFT .LE .0.0)  GO  TO  4606 

no  =  on  -  ALtFT 

IF  (I'D.  IT.  0.0)  GO  TO  4590 
NFwPT  =  IPT  ♦  ISP 

OFPTH  (Nfc'wPT)  =  ODPTH  ♦  ALF.FT  *  CHGD  /  FLOAT  (INC) 
ODPTK=  DEPTH  (NEWPT) 

CHGD  =  DEPTH  (IPPT)  -  OOPTH 
IPT  =  IP!  ♦  1 X  A  X 

4606  CONTINUE 

IF  (OO.FO.O.O)  GO  TO  4590 

an  =  no  /  aync 

A I  DC  =  CHGD  /  AN 
NSEG  =  AN  -  1.0  ♦  0.5 
IF  (NSEG.LE.O)  GO  TU  4590 
DO  4b50  LHN  =  1 .NSEG 
NErfPT  =  IFT  ♦  ISP 

uepth(newpt)  =  nnpTH  ♦  aidc  *  float  (lmn) 

IPT  =  IPT  ♦  i X AX 
<♦650  CONTINUE 
4590  CONTINUE 

IF  (IEND.FO.n  GO  TO  46*0 
IP  =  IP  +  1 
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IPT  =  IP 
I  X  H  =  I  X  A 
IYH  =  I Y  A 
60  TO  4560 
4b  7  0  TEND  =  1 

GO  TO  460b 
4660  I  X  H  =  I  X  A 
ISTY  =  IYH 
IYH  =  I Y  A 
60  TO  4560 
4670  ALEFT  =  AHS<f>U) 

GO  TO  4602 
467b  ALEFT  =  0.0 
GO  TO  4602 
4690  I  I  AX  =  1 

I  I  I  AX  =  I  X  AX 
ITOT=IYAX*IXAX 
DO  100  K  =  1 » I  TOT 

IF  (DEPTH (K)  .FO.  9999.90)  OF PTh < K ) =OE PTH ( «- 1 ) 

100  continue 

DO  4695  K  =  l.IYAX 

WRITE  (6,6040)  K , ( OF P TH ( 1 1 ) . 1  I  =  I  I  AX , 1 1 1  AX ) 

I I  AX  =  I  I  AX  ♦  1  X  AX 

III  AX  =  III  AX  ♦  I X AX 

6040  FORMAT  (*  *»ll<), (12F10.2) ,/) 

469b  CONTINUE 
RETURN 

NOw  wPITE  AXIS  VALUE  TO  DISK  FILE 

1220  WRITE ( TOUT. 6001 ) 

6001  FORMAT (»i  NUMBER  TOO  SHORT*,//) 

1230  WRITE  (IOUT,6000) 

6000  FORMAT  (*1  NUMBER  TOO  LONG*,//) 

1240  J  =  1  -  20 

IF  (J.LT.l)  J  =  1 
J2  =  J  ♦  40 

IF  (J2.LE.4096)  GO  TO  1290 
vJ  =  J?  -  40 

1290  WRITE  '  I  OUT  ,6005)  (A(KK)  *lvK=J»J2) 

6005  FORMAT  ( •  » ,41  Al ) 

STOP 

1310  WRITE  (I OUT, 60 10) 

6010  FORMAT  (»1  EXCEEDED  20000  POINTS  ON  OnF  CONTOUR*) 
STOP 

1330  WRITE  ( I OUT ,6015) 

6015  FORMAT  <*1  NO  IFPTh  BFTWEEN  SLASHES*) 

STOP 

1350  WRITE  ( I  OUT , 6020 ) 

6020  FORMAT  Cl  TWO  SIGNS  WITH  ONE  NUMBER*) 
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GO  TO  1?40 

1370  WHITE  ( I0UT.6()2t>) 

602b  FOPHAT  (»1  flO  SLASH  MFTwEEN  DEPTHS') 

I  OFF  =  0 
STOP 

13B0  WHITE  (IOUT»«G3S) 

603b  FOHPAT  ('1  MOPE  THAN  100  DEPTHS  ALONG  AN  AXIS*) 
STOP 

2000  CONTINUE 
GO  TO  200 
FM> 
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/* 

//  EXEC  ASMGCL  »PARM ,LKED= •LET«LIST,XREF»CALL* 
//ASM.SYSLIN  DO  UN  I T  = » SPACE  = 

//ASM.SYSIN  DD  * 


SORT  1 

CSFCT 

BEGIN 

SAVE 

(14,12) 

RALR 

11,0 

USING 

*,11 

LP 

12,13 

LA 

1 3,SAVFA 

ST 

13,ti(  12) 

ST 

12,4(13) 

B 

HFRE 

SAVE  A 

OS 

18F 

HERE 

LA 

1 »P  AkL 1ST 

LINK 

EP=SORT,MF=(E, ( 1 ) ) 

L 

1 3,SAVEA*4 

RETURN  (14,12) 

CNOP 

0,8 

PARLISC 

OC 

X  »fto  • 

DC 

AL3( ADLTST) 

DC 

X'OOOO* 

AI)L  1ST 

DC 

X  '0024* 

DC 

A ( SORT  CD ) 

DC 

A (STCDEO) 

DC 

A (RCDCD) 

DC 

A  <  R9CUED ) 

DC 

E  'O' 

DC 

F*0« 

DC 

X'OOOO* 

DC 

X • 6d90 • 

DC 

C'BALN' 

DC 

X  *FFOO ' 

DC 

C'AC' 

SORT  CO 

DC 

C*  SORT  FIELDS=(1,4,A,S»4,A) ,EORNAT=FI* 

STCDEO 

DC 

C*  ♦ 

RCDCD 

DC 

C*  RECORD  LFNGTh=12,TYPF.=F  • 

RDCUEC 

OC 

C*  • 

END 

/* 

//L*EO. 

syslhod 

DD  QSN=N31 0013. WORK (REF RAC) ,D1SP= (NEw,CATLG) 

//  SPACE=(TPK,  ( 10,;:,?fl)  »RLSF>  » DCB=LOADMOD , 

//  UNIT=3330V,MSVGP=USCQ 
//iKtD.SYSLlH  DD  DSN=SYS1 .Fl)PTLlH»DISP=SHk 
//  DD  OSN  =  Swl  .LINKLlt:.DISP  =  SHR 

//  D[;  OSN=ACAO.SUBL IH.OSOOS.OlSP=SHk,UNIT=, VOL=SFR= 

// 
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TAPE  TO  DISK  (TTD) 


The  program  TTD  is  an  Assembler  routine  that  copies  the  digitized 
depth  data  from  magnetic  tape  to  disk.  A  listing  of  the  JCL  and  source 
deck  follows  this  brief  description  of  how  to  use  the  program. 

The  //GO. INTAPE  card  specifies  the  tape  volume  to  be  stored.  The 
?????  in  the  VOL=SER=?????  parameter  is  replaced  by  the  number  on  the 
tape  reel.  For  example,  if  the  number  were  CRD26,  the  parameter  will 
be  VOLUME=SER=CRD26 . 

The  //GO.OUTTAPE  specifies  the  output  data  set  name,  under  which 
the  depth  data  is  to  be  stored.  See  JCL  Considerations  for  rules  on 
coding  the  data  set  names  (DSN) . 
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LISTING  OF  TTO  (TAPE  TO  DISK) 


//TTD  JOB  N3100138.TTD.TI ME =(»9) . 
//  USER=N310013»PASSWORD=CRD 
//SIP  1  EXEC  ASMFCLG 


//ASM.SYSIN  DO 

• 

COPY 

START 

0 

R 15 

EOU 

15 

R14 

EOU 

14 

R 1 3 

EQU 

i.3 

R 1 2 

EQU 

12 

R 1 1 

EQU 

11 

R3 

EQU 

3 

R4 

EOU 

4 

BEGIN 

STM 

R14.R12»12(R13) 

LR 

R12.R15 

USING 

BEGIN.R12 

B 

ARD 

SAVEA 

OS 

18F 

ARD 

LA 

Rll. SAVEA 

ST 

R 1 1 1 8 (H 13) 

ST 

R13»4 (R1 1 ) 

LR 

R13.R11 

L 

R11»=A (TAPE  IN- 1024) 

USING 

TAPE1N-1024.RH 

OPEN 

(TAPEIN, INPUT. TAPEOUT. OUTPUT) 

LOOP 

EQU 

• 

GET 

TAPEIN, RECORO 

PUT 

T  APEOUT .RECORD 

L 

R3»=F  *4096* 

LA 

R4.REC0R0 

LOOP! 

EQU 

• 

CLI 

0  (R4) ,C'F» 

BE 

EOF 

LA 

R4, 1 (R4) 

BCT 

R3.LOOP1 

B 

LOOP 

EOF 

EOU 

» 

PUT 

TAPEOUT.RECOHO 

CLOSE 

(TAPEIN, .TAPEOUT) 

L 

R13.4 (R13) 

LM 

R14,R12,12(R13) 

RR 

LTORG 

R14 

RECORD 

OS 

4096C 

tapein 

DCB 

OSORG=PS,MACRF=GM,ODNAMF=INTAPE,LRECL=4096»BLKSIZE=*096.X 

RECFM=F,EOOAD=EOF,BUFNO=l 

TAPE OUT 

OCR 

OSORG=PS,MACRF=PM,ODNAMEsOUTTAPE .LRECL*4096»  X 

11-98 
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} 


BLKSIZE=4096.RECFM=F 
END  BEGIN 

/* 

//GO.SYSUDUMP  DO  SYSOUT=A 

//GO. INTAPE  DO  OSN=SJC,UNIT=TAPEO»OISP=OLO»VOLUME=SER*?????» 

//  DCB=(HLKSIZE=4096»RECFM=F.EROPT=ACC*&UFNO=1) »LAREL=(1»BLP) 

//GO.OUTTAPE  DO  DSN=N31 0013. EXAMPLE »DISP= <NEW*CATLGtDELETE> . 

//  UNI T  =  3330V  « MS VGP=USCP »SPACE=(TKK» (5*2) tHLSU . 

//  DCB=(bLKSIZE=4096»RECFM=F ) 

/* 


— 3 

DEC  78  C  FJCO  ARO-13237* 15-GS  NL 

UNCLASSIFIED  TR-17-CR0  AR0  13237,1 


MICROtOPY  RESOLUTION  TEST  CHART 


NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


LISTING  DIGITIZED  DATA  (DIGCOPY) 

A  listing  of  the  digitized  depth  information  can  be  printed  from 
disk  by  the  Assember  program  DIGCOPY.  This  is  an  easy  means  of  veri¬ 
fying  the  digitized  data  by  comparing  the  printed  X,Y  coordinates  and 
the  corresponding  depths  with  the  depths  on  the  hydrographic  chart. 

The  data  set  name  referred  to  in  the  //GO. TAPE  card  (third  from 

i 

the  last  card)  is  the  same  data  set  name  used  in  the  //GO.OUTTAPE^  card 
in  the  TTD  program  previously  described.  (This  program,  DIGCOPY,  and 
TTD  were  written  by  Jim  Crabtree  to  aid  in  transferring  data  and  in 
error  checking. ) 


LISTING  OF  OIGCOPY  (DIGITIZED  TAPE) 


//DIGCOPY  JOB  (N310013**9) *0IGC0PY.TIM£s(*9)*MSGLEVEL*(l*l)* 
//  USER*N3lOOl3*PASSWORO*CRO 
//STEP1  EXEC  ASMFCLG 

//ASM.SYSIN  DO  * 

PRINT  NOGEN 

DISPLAY  START  0 

RO  EOU  0 

R1  EOU  1 

R2  EOU  2 

R3  EOU  3 

R4  EOU  4 

R5  EOU  5 

R6  EOU  6 

R7  EOU  7 

RB  EOU  6 

R9  EOU  9 

RIO  EOU  10 

Rll  EOU  11 

R12  EOU  12 

R13  EOU  13 

R14  EOU  14 

HIS  EOU  15 

BEGIN  STM  R14,R12*12(R13) 

LR  R12»RIS 

USING  BEGIN*R 12*R1 l 

L  Rll *BASE2 

B  ARD 

ASAVE  DS  18F 

BASE2  DC  A(BEGIN*4096) 

ARD  EOU  * 

LA  P10» ASAVE 

ST  R10.8(R13) 

ST  R13*ASAVE*4 

LR  R13»R10 

LA  R7* 1 

OPEN  (TAPE  IN* INPUT  * PR INTER* OUTPUT) 

SR  R3*R3 

LA  R5*LINE*10 

BAL  R10*HEAOER 

L OOP  EOU  * 

*  GET  TAPE  RECORD 


GET  TAPEIN,TAPE 
LA  R9*TAPE 


* 


ALOOP 


NOTNUM 


SPEC 


SPEC1 


* 


SPEC2 


START  SCAN 

LA  R8.4095 

LA  R8.i<R8) 

EOU  * 

CLI  0 (R9) *C *0 *  *  NUMBER  TEST 

BL  NOTNUM 

CLI  0<R9)*C»9* 

BH  NOTNUM 


TM  EX.X»FF« 

BO  ARD2 

B  AROl 

TEST  FOR  ALPHA 


EOU 

• 

CLI 

0(R9) *C'Z« 

BH 

SPECERR 

CLI 

0(R9)«C»S» 

BL 

SPEC 

B 

ALPHA 

EQU 

• 

CLI 

0<R9>  *C*R* 

BH 

SPECERR 

CLI 

0 (R9) «C*J« 

BL 

SPECl 

B 

ALPHA 

EOU 

• 

CLI 

0<R9) *C*  A* 

BL 

SPEC2 

CLI 

0IR9) » C* I » 

BH 

SPECERR 

B  ALPHA 


TEST  FOR  SPECIAL  CHARACTERS 


EOU 

• 

CLI 

0  <RV) »C V 

• 

SLASH 

BE 

ALPHA 

CLI 

0  <R9) *C»-» 

• 

MINUS 

BE 

MOVEXV 

CLI 

0  <R9J tC»*» 

• 

ASTER1CK 

BE 

ALPHA 

CLI 

0(R9) *C»*» 

• 

PLUS 

BE 

MOVtXY 

CLI 

0(R9»*C».* 

BE 

AR02 

11-102 


MOVEXY 

MOVfcXY 1 


AROl 

ARD2 

Y 

ALPHA 

ALPHA1 


H  SPECERP 


EOU 

• 

MV  I 

EX  *X • 00  * 

TM 

ALP»X«fF* 

B2 

MOVEXYl 

C 

R5*sA(LINE*10> 

BE 

MOVEXYl 

LA 

R3*1(R3) 

BAL 

R6. WRITER 

EOU 

« 

MV  I 

ALP»X*00* 

CLI 

XY*C»X' 

• 

X  VALUE 

BNE 

Y 

* 

NO 

LA 

R3* 1 (R3> 

MV  I 

X Y* C  *  Y • 

ENOUGH  SPACE  LEFT 

LA 

R4.LINE+130 

• 

SR 

R4  »RS 

C 

R4t=sF  *  12* 

BH 

AR02 

BCTR 

R3*  0 

BAL 

R6«WRITER 

LA 

P3»  l  (R3) 

B 

AR02 

EOU  * 

TM  ALP*X»FF» 

82  AR02 

LA  R3« l  (R3) 

MV I  EX«X*Ff» 

EOU  * 

MVC  O(l.Rb) »0(R9) 

LA  R5.1(R5) 

LA  R9*1<R9) 

B  TEST 

EOU  * 

MV  I  XY.C«X* 

B  ARD2 

EOU  * 

MV I  EX  *  X  *  00  * 

MV  I  ALP.X«FF* 

C  R5**A(LINE*10) 

BE  ALPHAi 

BAL  R6f WRITER 

EOU  ♦ 

LA  R3* i <R3) 

MVC  O(lfRS) « 0(R9> 

BAL  R6*wRITER 

CLI  0<R9)»C'F* 

BE  CLOSE 

LA  R9»1<R9) 


*  END  TEST 


TEST 


TEST1 

HEADER 


SPECERR 


WRITER 


EGlU 

• 

C 

R5«*A (LINE* 130) 

*  LINE  EXCEEOEO 

BL 

TEST1 

*  NO 

BAL 

R6* WRITER 

*  YES 

EOU 

• 

BCT 

RBtALOOP 

MV! 

XY.C*X« 

B 

LOOP 

EOU 

« 

MV  I 

PLINE*C»l» 

MV  I 

LINE«C»  » 

MVC 

LINE*l(131).LINfc 

MVC 

LINE*10 (26) **CL26* CHARACTER  DISPLAY  OF  TAPE* 

MVC 

LINE*1Q0<4>  »*CL4*PAGE* 

MVC 

LINE*104(4) » =XL4*40202120» 

ED 

LINE* 104(4) *PAGECT 

AP 

PA6ECT  *  *P» 1 » 

PUT 

PRINTER.PLINE 

MV  I 

PLINE»C*  • 

MVC 

LINE (132) .PLINE 

MVC 

LINE (116) ♦ HE  AO 

PUT 

PRtNTER.PLINE 

ZAP 

LINECT»=P'3* 

MV  I 

PLINE»C»0* 

MV  I 

LINE«C*  » 

MVC 

LIN£*1(131)*LINE 

BR 

RIO 

EQU 

• 

*  UNSPECIFIED  CHARACTER 

CL  I 

0  (R9) «X*00* 

*  BLOCK  FILLER 

BE 

TEST! 

BAL 

ROtWRITER 

MVC 

LINE (5) **CL5 **•**• • 

*  UNIDENTIFIED  CHARACTER 

BAL 

R6* WRITER 

TRT 

0 ( 1 *R9) »T  AB1 

BZ 

ALPHA 

MVC 

8YTES(2) *«XL 2*0000* 

MVZ 

BYTE1»0(R9) 

MVN 

BYTE2*0(R9) 

PACK 

BYTEl*8YTEl 

TR 

BYTES (2) »TAB2 

MVC 

LINE*10(2), BYTES 

MVC 

LINE* 14 (31) « *CL31 * (UNIDENTIFIED  CHARACTER  IN  HEX) 

BAL 

R6«WR1TER 

LA 

R9» 1 (R9) 

B 

TEST1 

EQU 

• 

CLC 

LINE (5) *sCL5'*****» 

BE 

W1 

MVC 

LINE (8) v«XL8* 40202020202021 20* 

CVO 

R7«00U8 

< 


ED 

LINE (6) .D0U8*4 

AH 

R7.R3 

wi 

ECU 

• 

PUT 

PRINTER. PLINE 

MV  I 

PLINE. C*  • 

MVC 

LINE (132) .PLINE 

AP 

LINECT. =P»I* 

CP 

LINECT. sP»60* 

BL 

WRITF.R1 

BAL 

H10.HEA0EH 

WRITER1 

EQU 

« 

LA 

R5.LINEM0 

SH 

R3.H3 

BH 

R6 

CLOSE 

EQU 

« 

CLOSE 

(TAPE  IN. .PRINTER) 

L 

R13.ASAVE+4 

LM 

H14.H12»12(R13) 

BH 

LTORG 

R14 

DOUB 

OS 

0 

TAPE 

OS 

4096C 

XY 

OC 

c»x* 

PLINE 

OS 

C 

LINE 

OS 

CL132 

HEAO 

OC 

CL49»REC0RD  *  X 

OC 

CL48'X  Y  X 

DC 

CL19»X  Y  X 

PAGECT 

OC 

PL2* 1  * 

LINECT 

OC 

PL2«0‘ 

BYTES 

OS 

0CL2 

BYTE  1 

OC 

X»0Q» 

BYTE  2 

OC 

X*00* 

EX 

OC 

X*00» 

ALP 

OC 

X*00» 

TAB2 

OC 

CL16*  0123456789ABC0EF  * 

TAB1 

DC 

256X*FF  * 

ORG 

TABI*74 

OC 

XL 7*00000000000000* 

ORG 

TAB1-.90 

OC 

XL8* 0000000000000000* 

ORG 

TA81*107 

OC 

XLS'OOOOOOOOOO* 

ORG 

TAB1M22 

DC 

XL 6* 000000000000* 

ORG 

TA81M93 

DC 

XL9* 0000 000000 00000000* 

ORG 

TAB1 *209 

OC 

XL9* 000000000000000000* 

ORG 

T  AB1*226 

•  ;%5 


11-105 


DC  XL8*  0000000000000000* 

ORG  TAB l *2*0 

DC  XL10«00000000000000000000* 

OPG  TAB1*2S6 

PRINTER  DCB  DSORG*PS»OONAME*L INES* LRECL*li3»BLKSIZE»133*RECFMs:FA 

MACRFsPM 

TAPEIN  OCR  DSORGsPS*DDNAMfcsTAPE*LR£CLe4096*BLK5IZE«4096*RECFM»F 

EODAD*CLOS£.MACRF*GM 
END  8EGJN 

/* 

//GO. LINES  OD  SYSOUT*A 
//GO.SYSUDUMP  DD  SYSOUTsA 
//GO. TAPE  OD  0SN=N310013. EXAMPLE. DISPsSHR 
/* 


